Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,126 G→T 7.3% intergenic (+109/‑89) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1260GT 7.3% 266.1 / 8.4 79intergenic (+109/‑89)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base T (2/3);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTC  >  NC_000913/3780978‑3781275
                                                                                                                                                     |                                                                                                                                                     
tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                                      <  2:402013/150‑1 (MQ=255)
tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                                      <  1:903035/150‑1 (MQ=255)
tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                                      <  1:2357298/150‑1 (MQ=255)
tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                                      <  1:1858514/150‑1 (MQ=255)
tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                                                      <  2:2272917/150‑1 (MQ=255)
 gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGTCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                                                                                                                                     <  2:1800133/150‑1 (MQ=255)
     cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTc                                                                                                                                                 <  1:2175108/150‑1 (MQ=255)
     cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTc                                                                                                                                                 <  1:1651274/150‑1 (MQ=255)
       cTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCAAACCCACGCCCTCCGGGTG‑AGGGGGCCGTTCGGCGCGGTATGTAATCCCGCACTCTGAAAACACACCCCACACATTTTTTCTCCCACGCCCCAccc                                                                                                                                               >  2:1392696/1‑149 (MQ=255)
       cTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCGCTCCCTAAGGGCG‑AGGGGACCGTTCGGCGCGGTTTGTACTCCCTCACTCTGAAACGACACAGCACTATTTTTTTCTCCCTCGCCCCTCCg                                                                                                                                               >  2:1392611/1‑150 (MQ=255)
        tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                                              <  2:1548313/150‑1 (MQ=255)
        tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                                              <  1:2388666/150‑1 (MQ=255)
        tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                                              <  2:1548354/150‑1 (MQ=255)
        tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                                              <  2:1548252/150‑1 (MQ=255)
         ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                                             <  2:1628400/150‑1 (MQ=255)
         ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                                             <  2:698266/150‑1 (MQ=255)
         ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                                             <  1:1722057/150‑1 (MQ=255)
         ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                                             <  2:698313/150‑1 (MQ=255)
         ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                                                             <  1:1722028/150‑1 (MQ=255)
           cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                                                           <  1:1572573/150‑1 (MQ=255)
           cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                                                           <  1:1109075/150‑1 (MQ=255)
           cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                                                           <  1:1454774/150‑1 (MQ=255)
           cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                                                           <  1:690327/150‑1 (MQ=255)
            tCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAAACCTCAACCACAAGGAGCGGGGGACCG‑TCGGCGCGGTTAGTAATCCCCCACTCTGAAACCACACCGCACTCATTTTTTCTCCCCCGCCCCACCcggggg                                                                                                                                          >  2:1778268/1‑148 (MQ=255)
              ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGGCCCCGCACCCATAACCACGCCCTCAGGGCG‑AGGGGGCCGTTCGGCGCGGTGAGGACGCCCACACGCGGAAACGACACCACCCACGTTTTTTCGCCCGCGCCCCACCGGGGtgag                                                                                                                                        >  2:1609191/1‑150 (MQ=255)
              ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACACCCACAAGGAG‑AGGGGGCCGTTCGGCGCTGTATTTACTCCCACACACAGTAACGACACCCCACAAGTTTTTTCTACCTCCCCCCACCCGGGagag                                                                                                                                        >  2:1609360/1‑150 (MQ=255)
              ccATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCTAACCCACAACCGCAGGGAG‑AGGGGACCGTTAGGCCCGGTATGTACTCCCCCACACTGAAAACACACCGCACTCATTTTTTCTCCCAAGCCCCACCCGGGagag                                                                                                                                        >  2:1609327/1‑150 (MQ=255)
                aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAggg                                                                                                                                      <  2:152682/150‑1 (MQ=255)
                aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAggg                                                                                                                                      <  2:110687/150‑1 (MQ=255)
                        aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑TGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                              <  1:1915020/150‑1 (MQ=255)
                        aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                              <  1:1788969/150‑1 (MQ=255)
                        aGGGAAGGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCGCGCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTg                                                                                                                              <  1:602602/150‑1 (MQ=255)
                         gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                                                             <  1:1605283/150‑1 (MQ=255)
                          ggAATGCGGCATAGTCGTTTTCGCCCCCAGACCCAAAACCACTACCCCAGCGAG‑AGGGGACCGTTAGGCGCTGTATGTACTCCCACACACAGAAAAGACAACCCACAAATTTTTTATCACGCGCACCAACCGGGAGAGGGGCGGGGTgtg                                                                                                                            >  2:1482217/1‑148 (MQ=255)
                          ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                                            <  1:2402130/150‑1 (MQ=255)
                          ggAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACTCCCTCAGGGAG‑AGGGGGCCCGTAGGCGCTGTATGTAAACCCGCACTCAGAAACGACACCGCACACTTTTTTTATCCCACGCCCCAACGGGGAGGGGGCCGGGGTGAg                                                                                                                            >  2:325850/1‑150 (MQ=255)
                           gattGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGTG‑AGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                           <  1:1115971/147‑1 (MQ=255)
                           gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                           <  1:1373493/150‑1 (MQ=255)
                           gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                           <  1:1080969/150‑1 (MQ=255)
                           gAATGCGGCATAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                           <  2:1667751/150‑1 (MQ=255)
                           gAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCTCGCCCGCCGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGACACGACCCCGCACTCTTTTTTTCTCCCGCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                           <  1:1115882/150‑1 (MQ=255)
                           gAATGCGGAATAGTGGTTTGCGCCCCATCACCCTAACCCTCTCCCTCCGGTTG‑TGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                           <  1:1036711/150‑1 (MQ=255)
                                 ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                                                     <  2:1807778/150‑1 (MQ=255)
                                               cgcCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCa                                                                                                       >  2:2374897/1‑150 (MQ=255)
                                               cgcCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCa                                                                                                       >  2:179137/1‑150 (MQ=255)
                                                            tAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                          >  1:870248/1‑150 (MQ=255)
                                                            tAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                          >  2:1129835/1‑150 (MQ=255)
                                                            tAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                          >  1:1179647/1‑150 (MQ=255)
                                                             aaCCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCtc                                                                                         >  2:1258529/1‑149 (MQ=255)
                                                               cccTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                       >  1:2152073/1‑150 (MQ=255)
                                                                                   gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                                    >  1:201589/1‑150 (MQ=255)
                                                                                   gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                                    >  2:1782751/1‑150 (MQ=255)
                                                                                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                                                  >  1:1877243/1‑150 (MQ=255)
                                                                                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                                                              >  1:924485/1‑150 (MQ=255)
                                                                                          ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAACAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                             >  2:2348833/1‑150 (MQ=255)
                                                                                          ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                             >  2:2304481/1‑150 (MQ=255)
                                                                                          ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                             >  1:804051/1‑150 (MQ=255)
                                                                                          ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                             >  2:1223257/1‑150 (MQ=255)
                                                                                          ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                             >  1:804063/1‑150 (MQ=255)
                                                                                                 cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                                                      >  2:198534/1‑150 (MQ=255)
                                                                                                 cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                                                      >  1:1015138/1‑150 (MQ=255)
                                                                                                           tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                                            >  1:1604044/1‑150 (MQ=255)
                                                                                                                         aCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAtt                              >  2:468853/1‑150 (MQ=255)
                                                                                                                         aCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAtt                              >  2:1823898/1‑150 (MQ=255)
                                                                                                                               ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCa                        >  2:260541/1‑150 (MQ=255)
                                                                                                                                cGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaa                       >  1:125918/1‑150 (MQ=255)
                                                                                                                                 ccACTTTTTCTTCCACCCCCCCCCCCCGGGGGGGGGGGGGGGGGGGCGGGCCCAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaaa                      <  1:2014471/149‑1 (MQ=255)
                                                                                                                                     tctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTCCTTTACGAACCAGAAATTCCGCCAAATACt                  >  2:1688937/1‑150 (MQ=255)
                                                                                                                                       tttttttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTgg                >  2:464741/1‑150 (MQ=255)
                                                                                                                                           tttCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAc            >  2:655670/1‑150 (MQ=255)
                                                                                                                                           cccggcccTCTCCCCGGAGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                                                   >  2:1373070/6‑111 (MQ=38)
                                                                                                                                           cccggcccTCTCCCCGGAGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                                                   <  1:1373070/106‑1 (MQ=38)
                                                                                                                                            ccggcccTCTCCCCGGAGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAca           >  1:1335069/5‑150 (MQ=255)
                                                                                                                                            ccggcccTCACCCCGGAGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAATTTATTCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAAAACGTACTTTACGAACCAGAAATTCCGCCAAAAAATGGCAAca           >  1:1335035/5‑150 (MQ=255)
                                                                                                                                               tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcatca        >  1:51412/1‑150 (MQ=255)
                                                                                                                                                  cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc     >  2:947411/1‑150 (MQ=255)
                                                                                                                                                  cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc     >  1:612225/1‑150 (MQ=255)
                                                                                                                                                  cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc     >  1:1816968/1‑150 (MQ=255)
                                                                                                                                                  cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc     >  1:1031557/1‑150 (MQ=255)
                                                                                                                                                     gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc  >  1:1660038/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAG‑AGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTC  >  NC_000913/3780978‑3781275

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑