Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,713 A→T 8.3% intergenic (‑46/‑245) gcvT ← / → yqfI aminomethyltransferase/protein YqfI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,7130AT 8.3% 314.5 / 12.3 106intergenic (‑46/‑245)gcvT/yqfIaminomethyltransferase/protein YqfI
Reads supporting (aligned to +/- strand):  ref base A (61/36);  new base T (7/2);  total (68/38)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

GATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGC  >  NC_000913/3050566‑3050862
                                                                                                                                                   |                                                                                                                                                     
gATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAt                                                                                                                                                     >  2:1068891/1‑150 (MQ=255)
gATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAt                                                                                                                                                     >  1:1315747/1‑150 (MQ=255)
  tttGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTTTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGTGGCTGACAACTTTTTCATGGTGAAATTa                                                                                                                                                   <  2:367681/150‑1 (MQ=255)
    tGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATc                                                                                                                                                 >  1:427732/1‑150 (MQ=255)
     gCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATca                                                                                                                                                >  2:137175/1‑150 (MQ=255)
        aaCCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATcacca                                                                                                                                             >  1:575321/1‑150 (MQ=255)
          ccGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCAGACAACTTTTTCAAGGTGAAATTATCACCACg                                                                                                                                           >  1:63569/1‑150 (MQ=255)
              aaTGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGATAATTTCACCATGaaaa                                                                                                                                       <  1:1027333/150‑2 (MQ=255)
                 gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCa                                                                                                                                    >  2:1618513/1‑150 (MQ=255)
                   aGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGc                                                                                                                                  >  2:1391356/1‑150 (MQ=255)
                       gCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCg                                                                                                                              <  2:684423/150‑1 (MQ=255)
                        catcatCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCgg                                                                                                                             <  1:1633030/150‑1 (MQ=255)
                           catcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGGGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCGCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCCCGAAACCCCGCAGCGgcgc                                                                                                                          >  2:1054221/1‑150 (MQ=255)
                            atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGcc                                                                                                                         <  1:875926/150‑1 (MQ=255)
                            atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGcc                                                                                                                         <  2:1031992/150‑1 (MQ=255)
                             tcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCa                                                                                                                        >  1:1386859/1‑150 (MQ=255)
                                aGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTc                                                                                                                     >  1:730226/1‑150 (MQ=255)
                                           tCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCaa                                                                                                          <  1:877163/150‑1 (MQ=255)
                                           tCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCaa                                                                                                          >  1:656400/1‑150 (MQ=255)
                                             caccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCATCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATAACAACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                                                        >  1:973839/1‑150 (MQ=255)
                                                 aTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGTGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAAt                                                                                                    <  1:214454/150‑1 (MQ=255)
                                                   gCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCTCCAACGCAATCg                                                                                                  <  1:81944/150‑1 (MQ=255)
                                                     gAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTTTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTACCGGTCACCAACGCAATCGtt                                                                                                >  2:47598/1‑150 (MQ=255)
                                                       gcgcCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtct                                                                                              <  1:1068891/150‑1 (MQ=255)
                                                       gcgcCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtct                                                                                              <  1:290974/150‑1 (MQ=255)
                                                          ccgcAAAGCGTGTGTTGTTGGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCCCCACGAAACCCAGCATCGGAGACACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                                                           >  2:529284/1‑150 (MQ=255)
                                                               aaGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                      <  2:1565371/150‑1 (MQ=255)
                                                                 gCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTg                                                                                    <  1:137175/150‑1 (MQ=255)
                                                                     gtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACt                                                                                <  2:48741/150‑1 (MQ=255)
                                                                         tgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacc                                                                            <  2:1321682/150‑1 (MQ=255)
                                                                               tACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGaa                                                                      >  1:900717/1‑150 (MQ=255)
                                                                                 cAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAAc                                                                    >  1:1631530/1‑150 (MQ=255)
                                                                                 cAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAAc                                                                    >  2:626976/1‑150 (MQ=255)
                                                                                 cAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACAACGCAACCCCGCATCGGAGCCCCTCCCGGTCCCCAACGCACTCGTTCACTTTTGCCTGAACTTACACCCGAAAc                                                                    >  2:1210418/1‑150 (MQ=255)
                                                                                   aaGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGt                                                                                                 <  1:22518/119‑1 (MQ=255)
                                                                                   aaGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGt                                                                                                 >  2:22518/1‑119 (MQ=255)
                                                                                     ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                >  1:1027142/1‑150 (MQ=255)
                                                                                     ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                <  2:1562438/150‑1 (MQ=255)
                                                                                     ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                >  1:543410/1‑150 (MQ=255)
                                                                                     ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                <  1:378458/150‑1 (MQ=255)
                                                                                     ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                >  1:1485363/1‑150 (MQ=255)
                                                                                         tCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGtt                                                            <  1:1391356/150‑1 (MQ=255)
                                                                                          cTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTa                                                           <  1:30367/150‑1 (MQ=255)
                                                                                           tgttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCATGAACTTAACACAGAAAAAGACTGTTaa                                                          >  1:693171/1‑150 (MQ=255)
                                                                                            gttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                         <  2:427732/150‑1 (MQ=255)
                                                                                            gttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                         >  1:1547308/1‑150 (MQ=255)
                                                                                               gtgCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCAt                                                      <  1:242058/150‑1 (MQ=255)
                                                                                                tgCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATa                                                     >  1:781689/1‑150 (MQ=255)
                                                                                                 gCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                         >  2:1347275/1‑145 (MQ=255)
                                                                                                 gCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                         <  1:1347275/145‑1 (MQ=255)
                                                                                                 gCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCCAAACAGACTGTTAACCATaa                                                    >  2:978538/1‑150 (MQ=255)
                                                                                                      cTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTaa                                               >  2:1597050/1‑150 (MQ=255)
                                                                                                      cTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTaa                                               >  2:1380395/1‑150 (MQ=255)
                                                                                                      cTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGAATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGACCTTACCACCGAAACCGACTGTTAACCATAAGGTaa                                               >  2:1597051/1‑150 (MQ=255)
                                                                                                       ttGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCCCTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTaaa                                              >  2:574579/1‑150 (MQ=255)
                                                                                                        tGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTTCCACCGAAACAGACTGTTAACCATAAGGTaaaa                                             >  1:1236803/1‑150 (MQ=255)
                                                                                                        tGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTaaaa                                             >  1:824645/1‑150 (MQ=255)
                                                                                                         gTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAt                                            <  2:270381/150‑1 (MQ=255)
                                                                                                           ccTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTg                                          >  1:542415/1‑150 (MQ=255)
                                                                                                           ccTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTg                                          >  1:1123759/1‑150 (MQ=255)
                                                                                                             tCATTTAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGat                                        <  2:255227/150‑1 (MQ=255)
                                                                                                             tCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGat                                        >  1:1002210/1‑150 (MQ=255)
                                                                                                             tCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGat                                        >  1:927698/1‑150 (MQ=255)
                                                                                                             tCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGat                                        >  1:933408/1‑150 (MQ=255)
                                                                                                               aTTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatca                                      >  1:784220/1‑150 (MQ=255)
                                                                                                                ttGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatcat                                     >  1:1339361/1‑150 (MQ=255)
                                                                                                                  gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGAtcatca                                   >  1:276584/1‑150 (MQ=255)
                                                                                                                  gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGAtcatca                                   >  1:1055937/1‑150 (MQ=255)
                                                                                                                   aaTAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATcac                                  <  2:1257592/150‑1 (MQ=255)
                                                                                                                     tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAt                                >  1:31350/1‑150 (MQ=255)
                                                                                                                       aGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTTATCATCACAttt                              <  2:730226/150‑2 (MQ=255)
                                                                                                                       aGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTa                              >  1:1188576/1‑150 (MQ=255)
                                                                                                                       aGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTa                              >  1:1355903/1‑150 (MQ=255)
                                                                                                                        gCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAg                             >  2:1457709/1‑150 (MQ=255)
                                                                                                                         cGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGc                            <  2:266706/150‑1 (MQ=255)
                                                                                                                         cGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGc                            <  2:656400/150‑1 (MQ=255)
                                                                                                                          ggggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCt                           >  1:1431008/1‑150 (MQ=255)
                                                                                                                           gggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCtt                          >  2:156308/1‑150 (MQ=255)
                                                                                                                               tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                      <  2:543410/150‑1 (MQ=255)
                                                                                                                               tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAATAGCTTATgg                      >  2:1423375/1‑150 (MQ=255)
                                                                                                                                gACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGt                     >  2:741074/1‑150 (MQ=255)
                                                                                                                                    tgggtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaa                 >  2:687944/5‑150 (MQ=255)
                                                                                                                                     cTTTTTCATGGTGATATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa                <  1:737888/150‑1 (MQ=255)
                                                                                                                                     cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa                <  2:1547308/150‑1 (MQ=255)
                                                                                                                                     cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa                >  2:970556/1‑150 (MQ=255)
                                                                                                                                      ggtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaa               >  1:948820/3‑150 (MQ=255)
                                                                                                                                       gtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaaa              >  1:1299566/2‑150 (MQ=255)
                                                                                                                                       gtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaaa              >  2:393252/2‑150 (MQ=255)
                                                                                                                                         ttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATg            >  1:851459/1‑150 (MQ=255)
                                                                                                                                         ttCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAGAATg            >  1:265962/1‑150 (MQ=255)
                                                                                                                                           cgtGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCa          >  2:214553/3‑150 (MQ=255)
                                                                                                                                            aTGGTGAAATTATCACCACGAAACCCAGCATCGGATCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaa         <  1:47598/150‑1 (MQ=255)
                                                                                                                                              ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa       >  2:253637/1‑150 (MQ=255)
                                                                                                                                              ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa       >  2:1310151/1‑150 (MQ=255)
                                                                                                                                              ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCAAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa       >  1:1002067/1‑150 (MQ=255)
                                                                                                                                                 gAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATc    >  1:1186270/1‑150 (MQ=255)
                                                                                                                                                 gAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATc    >  1:1280472/1‑150 (MQ=255)
                                                                                                                                                  aaaTTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATcg   <  2:1577748/150‑1 (MQ=255)
                                                                                                                                                   aaTTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATcgc  >  1:134342/1‑150 (MQ=255)
                                                                                                                                                   aaTTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATcgc  >  1:536000/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGC  >  NC_000913/3050566‑3050862

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