Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 557,009 C→T 6.7% V245M (GTG→ATG)  folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913557,0090CT 6.7% 186.0 / 2.7 60V245M (GTG→ATG) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

AAAATGTCGCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACG  >  NC_000913/556861‑557153
                                                                                                                                                    |                                                                                                                                                
aaaaTGTCGCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGGCTTCAAacac                                                                                                                                                  >  6:268410/1‑149 (MQ=255)
 aattGTCGCCATGTTACTCATCCTGTGTATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACacg                                                                                                                                                 <  8:325037/146‑1 (MQ=255)
 aaaTGTCGCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAacac                                                                                                                                                  <  8:45057/148‑1 (MQ=255)
 aaaTGTCGCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACacg                                                                                                                                                 <  5:270493/149‑1 (MQ=255)
   aTGTCGCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACacgac                                                                                                                                               <  1:232435/149‑1 (MQ=255)
        gCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGc                                                                                                                                          >  4:41681/1‑149 (MQ=255)
         ccATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGcc                                                                                                                                         >  2:174111/1‑149 (MQ=255)
             gTTACTCATCCTGTGGATCATGATATTCAACGCACGACTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCACAGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGccccc                                                                                                                                      >  6:213994/1‑146 (MQ=255)
              ttACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACaa                                                                                                                                    >  4:124833/1‑149 (MQ=255)
                  tCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACATGACGTCGCCCACAACttt                                                                                                                                >  1:728/1‑149 (MQ=255)
                  tCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACttt                                                                                                                                >  7:152116/1‑149 (MQ=255)
                    aTCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACttt                                                                                                                                <  6:150125/147‑1 (MQ=255)
                    aTCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACttt                                                                                                                                <  7:3173/147‑1 (MQ=255)
                          tgGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCAtttt                                                                                                                        >  6:203051/1‑149 (MQ=255)
                            gATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTc                                                                                                                       <  3:7240/148‑1 (MQ=255)
                             aTCATGATATTCAACGACCGGCTGTCGCGTGTTTTCAATCAGCGTGGCAAACGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGGCGCCCACAACTTTGCCATTTTCCa                                                                                                                     <  5:97828/149‑1 (MQ=255)
                                 tGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGCCg                                                                                                                 <  7:210124/149‑1 (MQ=255)
                                 tGATATTCAACGCACGCCTGTAGAGTGTTTTCAATCAGCGTGGCAACCGTCATCGGTCCAACGCAGACTTGAACAGGAGTAATGTATGATGAGCTTTTATCCTCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACg                                                                                                                 >  3:107235/1‑149 (MQ=255)
                                  gATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACgg                                                                                                                >  4:175082/1‑149 (MQ=255)
                                   atatTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACgg                                                                                                                <  5:150285/148‑1 (MQ=255)
                                    tatTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACATGACGTCGCCCACAACTTTGCCATTTTCCAGACGGtt                                                                                                              <  2:728/149‑1 (MQ=255)
                                      ttCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGa                                                                                                            <  2:170894/149‑1 (MQ=255)
                                       tCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGAt                                                                                                           >  7:169356/1‑149 (MQ=255)
                                           cgcacgcCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCg                                                                                                       <  2:185696/149‑1 (MQ=255)
                                               cgcCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACCGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACat                                                                                                   <  1:217762/149‑1 (MQ=255)
                                               cgcCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACat                                                                                                   <  2:305220/149‑1 (MQ=255)
                                               cgcCTGTAGCGTGTTTTCAATCAGCGAGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGGCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACat                                                                                                   <  8:290648/149‑1 (MQ=255)
                                                   tGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACatcaat                                                                                               >  7:132302/1‑149 (MQ=255)
                                                     tAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCTGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACatcaatc                                                                                              <  1:258259/148‑1 (MQ=255)
                                                     tAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACATGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACatcaatca                                                                                             >  4:243966/1‑149 (MQ=255)
                                                       gCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGCCGGTTGATGCCGACATCAATcac                                                                                            >  7:311957/1‑148 (MQ=255)
                                                        cGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACaa                                                                                          <  2:24331/149‑1 (MQ=255)
                                                            ttttCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTgc                                                                                      <  5:183121/149‑1 (MQ=255)
                                                             tttCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTgcg                                                                                     >  7:138067/1‑149 (MQ=255)
                                                             tttCAATCAGCGTGGCAACCGTCATCGGGACAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTgcg                                                                                     >  4:174509/1‑149 (MQ=255)
                                                                 aaTCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCct                                                                                  <  2:307652/148‑1 (MQ=255)
                                                                     aGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCttt                                                                             >  7:303487/1‑149 (MQ=255)
                                                                     aGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCttt                                                                             >  8:14414/1‑149 (MQ=255)
                                                                       cGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGa                                                                           <  5:151402/149‑1 (MQ=255)
                                                                         tGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATc                                                                         <  3:124833/149‑1 (MQ=255)
                                                                            cAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAg                                                                      >  8:173577/1‑149 (MQ=255)
                                                                                  tCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTATTTGATCCAGTCACCg                                                                <  2:247894/149‑1 (MQ=255)
                                                                                    aTCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGATGCGCGTTTAGCCGCGTCTTCAAACACGACTTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCgg                                                               <  3:298368/148‑1 (MQ=255)
                                                                                     tCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCggg                                                              >  3:236330/1‑148 (MQ=255)
                                                                                          ccAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATaaa                                                        <  8:303487/149‑1 (MQ=255)
                                                                                            aaCGCCGCCGGGTACAGGCGTAATGTATGAGGCGCGGTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGc                                                      <  4:5134/149‑1 (MQ=255)
                                                                                            aaCGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGc                                                      <  5:268410/149‑1 (MQ=255)
                                                                                              cgccgcCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCt                                                    <  3:41681/149‑1 (MQ=255)
                                                                                                 cgccgGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGc                                                 >  7:168441/1‑149 (MQ=255)
                                                                                                 cgccgGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGc                                                 <  8:8565/149‑1 (MQ=255)
                                                                                                         aCAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAAc                                         <  8:169356/149‑1 (MQ=255)
                                                                                                           aGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACgg                                       <  2:214225/149‑1 (MQ=255)
                                                                                                                    gTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTGCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCAGACGGCAACGATc                               >  3:1157/1‑148 (MQ=255)
                                                                                                                      aTGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCaa                             >  7:165487/1‑148 (MQ=255)
                                                                                                                           gcgcGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACCATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCaa                             >  5:165842/1‑143 (MQ=255)
                                                                                                                             gcgTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGCTGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCCCCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCgg                     >  5:218620/1‑149 (MQ=255)
                                                                                                                                tttAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCAt                  <  5:203051/149‑1 (MQ=255)
                                                                                                                                 ttAGCCGCGTCTTCAAACAAGACGTCGACCACAACTATGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACCATCAAGAGCTCGGCagt                 >  3:25292/1‑147 (MQ=255)
                                                                                                                                        cgTCTTCAAACATGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACg          <  3:243966/149‑1 (MQ=255)
                                                                                                                                         gTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACAGCAACGATCAATAGATCGGCATTTTCTACGt         >  1:98410/1‑149 (MQ=255)
                                                                                                                                          tCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACAGCAACGATCAATAGATCGGCATTTTCTACGtg        <  2:98410/149‑1 (MQ=255)
                                                                                                                                             tCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATaaa                                                        <  7:271425/98‑1 (MQ=255)
                                                                                                                                             tCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATaaa                                                        >  8:271425/1‑98 (MQ=255)
                                                                                                                                              cAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGtgatga    >  3:277220/1‑149 (MQ=255)
                                                                                                                                               aaaCACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGAc   >  5:157376/1‑149 (MQ=255)
                                                                                                                                                   acgacgTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACg  >  3:307007/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AAAATGTCGCCATGTTACTCATCCTGTGGATCATGATATTCAACGCACGCCTGTAGCGTGTTTTCAATCAGCGTGGCAACCGTCATCGGGCCAACGCCGCCGGGAACAGGCGTAATGTATGAGGCGCGTTTAGCCGCGTCTTCAAACACGACGTCGCCCACAACTTTGCCATTTTCCAGACGGTTGATGCCGACATCAATCACAATTGCGCCTTCTTTGATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACG  >  NC_000913/556861‑557153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑