Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,192,685 A→C 5.7% R363R (CGT→CGG mnmA ← tRNA‑specific 2‑thiouridylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,192,6850AC 5.7% 308.2 / 7.4 105R363R (CGT→CGGmnmAtRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base A (49/50);  new base C (3/3);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 107 total aligned reads displayed.

ATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCGG  >  NC_000913/1192536‑1192830
                                                                                                                                                     |                                                                                                                                                 
aTCACAATGCCCCTGGTGAGCGAGTTGTTGCCCCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTATGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGa                                                                                                                                                   >  2:1069766/1‑150 (MQ=255)
aTCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGa                                                                                                                                                   >  2:712867/1‑150 (MQ=255)
 tCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGAc                                                                                                                                                  >  2:1110443/1‑150 (MQ=255)
    caATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgct                                                                                                                                               <  1:985067/150‑1 (MQ=255)
    caATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACAGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGATTACTTGTTGTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGCCgct                                                                                                                                               >  2:770885/1‑150 (MQ=255)
     aaTGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgctg                                                                                                                                              <  2:1483874/150‑1 (MQ=255)
      aTGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgctgc                                                                                                                                             <  1:1197886/150‑1 (MQ=255)
        gCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTc                                                                                                                                           <  2:249449/150‑1 (MQ=255)
         ccccTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCa                                                                                                                                          <  2:1208739/150‑1 (MQ=255)
             tGGTGAGCGAGTTGTGGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGTCTGCCgctgccggtctg                                                                                                                                      <  2:1286410/150‑8 (MQ=255)
              ggTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCaataat                                                                                                                                     <  2:130466/150‑1 (MQ=255)
                 gAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgctgc                                                                                                                                             >  2:953201/1‑139 (MQ=255)
                 gAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACgctgc                                                                                                                                             <  1:953201/139‑1 (MQ=255)
                  aGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCTCGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCg                                                                                                                                 <  1:147069/150‑1 (MQ=255)
                  aGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCg                                                                                                                                 <  1:880117/150‑1 (MQ=255)
                    ggAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCgtt                                                                                                                               <  1:299222/149‑3 (MQ=255)
                     gAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGTCTGCCgctgccgg                                                                                                                                          <  1:1289763/138‑4 (MQ=255)
                     gAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGTCTGCCgctgccgg                                                                                                                                          >  2:1289763/1‑135 (MQ=255)
                       gttgttGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCAcc                                                                                                                            <  2:245393/150‑1 (MQ=255)
                        ttgttgCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCg                                                                                                                           <  1:1110443/150‑1 (MQ=255)
                          gttgCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCTACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAg                                                                                                                         >  1:761689/1‑150 (MQ=255)
                           ttgCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAgg                                                                                                                        >  2:865936/1‑150 (MQ=255)
                            tgCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGc                                                                                                                       <  1:1315889/150‑1 (MQ=255)
                              cACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGcac                                                                                                                     >  2:1320239/1‑150 (MQ=255)
                              cACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGcac                                                                                                                     >  2:1424544/1‑150 (MQ=255)
                                   ggCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAgg                                                                                                                        >  2:476007/1‑142 (MQ=255)
                                   ggCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAgg                                                                                                                        <  1:476007/142‑1 (MQ=255)
                                      gTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAATATATTAATCAGACAGGAAGCTGAAGAATATGATAAATAATTCTTCCACCTATGGACTCTTCAca                                                                                                             <  2:35408/150‑2 (MQ=255)
                                      gTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCAcc                                                                                                             <  1:730469/150‑1 (MQ=255)
                                           gACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTtat                                                                                                        <  1:1194302/150‑1 (MQ=255)
                                             cTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATag                                                                                                      <  2:381196/150‑1 (MQ=255)
                                             cTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATag                                                                                                      >  2:686459/1‑150 (MQ=255)
                                             cTGACAAATACCGGCCAGGGAGAGGGTTATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATag                                                                                                      <  1:1385159/150‑1 (MQ=255)
                                             cTGACAAATACCCGCCAGGGCGAGGGTGATGTCATAGAAATTCTTTGCCACCTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGCCAGACCCTGCTCAATAATACCGCCACCCAGGCACACTTCACCGTTATag                                                                                                      >  2:334754/1‑150 (MQ=255)
                                               gACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaa                                                                                                    >  2:734346/1‑150 (MQ=255)
                                                 cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagatga                                                                                                  >  1:239622/1‑150 (MQ=255)
                                                 cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaaga                                                                                                  >  1:412885/1‑150 (MQ=255)
                                                 cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaaga                                                                                                  >  1:176500/1‑150 (MQ=255)
                                                 cAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATagaaga                                                                                                  >  1:834903/1‑150 (MQ=255)
                                                  aaaTACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGAAAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGAc                                                                                                 <  2:1148886/150‑1 (MQ=255)
                                                   aaTACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACg                                                                                                >  1:99536/1‑150 (MQ=255)
                                                    aTACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCAGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACgg                                                                                               >  2:761686/1‑150 (MQ=255)
                                                     tACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGc                                                                                              >  1:1079704/1‑150 (MQ=255)
                                                      aCCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACTGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCa                                                                                             >  2:43093/1‑150 (MQ=255)
                                                       ccGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCCCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAg                                                                                            <  1:120096/150‑1 (MQ=255)
                                                         ggCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGAc                                                                                          <  1:1424544/150‑1 (MQ=255)
                                                          gCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACt                                                                                         <  2:763541/150‑1 (MQ=255)
                                                            cAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTgg                                                                                       <  2:1242915/150‑1 (MQ=255)
                                                                   aGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCg                                                                                >  2:877606/1‑150 (MQ=255)
                                                                   aGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCg                                                                                <  1:712867/150‑1 (MQ=255)
                                                                   aGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCg                                                                                >  2:1531812/1‑150 (MQ=255)
                                                                    gggTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGt                                                                               >  1:1581376/1‑150 (MQ=255)
                                                                    gggTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGt                                                                               >  2:1627809/1‑150 (MQ=255)
                                                                          tGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCAcggc                                                                         <  1:1531812/150‑1 (MQ=255)
                                                                          tGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCAcggc                                                                         >  1:1210906/1‑150 (MQ=255)
                                                                           gTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGc                                                                                              >  2:1097003/1‑128 (MQ=255)
                                                                           gTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGc                                                                                              <  1:1097003/128‑1 (MQ=255)
                                                                            tCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCAcggcgg                                                                       <  2:365173/150‑1 (MQ=255)
                                                                            tCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCAcggcgg                                                                       >  1:971938/1‑150 (MQ=255)
                                                                             cATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACggcggc                                                                      >  2:535002/1‑150 (MQ=255)
                                                                                aGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCCCCGAGGCACACTTCACCGTTATTGAAGACGGCAGCCTGGCCCGGCGTCACCGCGGCAAc                                                                   >  2:1420428/1‑150 (MQ=255)
                                                                                  tAATTCTTTGCCACGGTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCg                                                                 <  1:770885/150‑1 (MQ=255)
                                                                                    aTTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGt                                                               <  2:102172/150‑1 (MQ=255)
                                                                                     ttCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGtt                                                              >  1:730291/1‑150 (MQ=255)
                                                                                         ttGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATc                                                          <  2:1210906/150‑1 (MQ=255)
                                                                                           gCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGa                                                        >  2:1062431/1‑150 (MQ=255)
                                                                                            ccACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGAAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGaa                                                       >  2:556488/1‑150 (MQ=255)
                                                                                             cACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGaaa                                                      >  1:1183059/1‑150 (MQ=255)
                                                                                             cACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGaaa                                                      <  2:1414876/150‑1 (MQ=255)
                                                                                              aCGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGaaaa                                                     <  2:242477/150‑1 (MQ=255)
                                                                                               cGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAAt                                                    <  1:1410849/150‑1 (MQ=255)
                                                                                                gTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAAta                                                   <  1:1627809/150‑2 (MQ=255)
                                                                                                    aCTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACtt                                               <  1:968770/150‑1 (MQ=255)
                                                                                                     cTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTc                                              <  2:163638/150‑1 (MQ=255)
                                                                                                        cTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAAt                                           >  1:522786/1‑150 (MQ=255)
                                                                                                        cTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAAt                                           <  2:1079704/150‑1 (MQ=255)
                                                                                                           ccTTGTTTTAAGTAAAGATAATAATCAGACCCGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGCCACACTTCACCGTTATAGAAGACCGCAGACTCGCCCGGCGTCACGGCGGCAACCGGATCAACGAAAATCACTTAAATgag                                        >  1:163728/1‑148 (MQ=255)
                                                                                                               gTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGatca                                    >  1:129599/1‑150 (MQ=255)
                                                                                                               gTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATAACGCCACCGAAGCACACTTCCCCGTTATAGACCACGGCAGACTGGCCCGTCGTAACGGCGGGAACCGGTTCAGCGAAAATCGCGTCAATGCGatcc                                    >  2:363488/1‑149 (MQ=255)
                                                                                                                 tttAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGatcatc                                  >  1:1026981/1‑150 (MQ=255)
                                                                                                                 tttAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGatcatc                                  >  2:261678/1‑150 (MQ=255)
                                                                                                                     aGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTcc                              >  2:1293704/1‑150 (MQ=255)
                                                                                                                     aGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTcc                              <  2:1581376/150‑1 (MQ=255)
                                                                                                                      gTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCa                             >  1:718205/1‑150 (MQ=255)
                                                                                                                       tAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAg                            <  1:1047572/150‑1 (MQ=255)
                                                                                                                       tAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAg                            <  1:1552048/150‑1 (MQ=255)
                                                                                                                        aaaGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGTCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAgc                           >  1:588641/1‑150 (MQ=255)
                                                                                                                         aaGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAgcg                          <  1:686459/150‑1 (MQ=255)
                                                                                                                             taataaTCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCtt                      >  1:1532953/1‑150 (MQ=255)
                                                                                                                             taataaTCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCtt                      >  1:982377/1‑150 (MQ=255)
                                                                                                                                 aatCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACg                  <  1:43666/150‑1 (MQ=255)
                                                                                                                                  atCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACgg                 >  1:912282/1‑150 (MQ=255)
                                                                                                                                     aGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGc              <  2:176500/150‑1 (MQ=255)
                                                                                                                                        ccGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTCTAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAg           <  1:1069766/150‑1 (MQ=255)
                                                                                                                                         gagcagcgtctgCCGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAgg          >  1:1614132/3‑150 (MQ=255)
                                                                                                                                          ggcagcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAggg         >  2:607662/1‑150 (MQ=255)
                                                                                                                                           gcagcggcagACGCTGCTCAATAATACCGCCACCGAGGCAGACTTCACCGTTATAGAAGAGGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGAGGGTGCAAGGGa        <  1:334754/150‑1 (MQ=255)
                                                                                                                                              gcggcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGatat     <  1:837082/150‑1 (MQ=255)
                                                                                                                                                ggcagACGCTGCTCAATAATACCGCCACCGAGGCGCACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTGATCGAAAATCACTTCAATGCGATCGTCGTCCAGCGCCTTGCCGGTGGAAGGGATATCg   >  2:819948/1‑150 (MQ=255)
                                                                                                                                                 gcagACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCgg  <  2:412885/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAATACCGGCCAGGGCGAGGGTGATGTCATAGTAATTCTTTGCCACGTTCACTGCTTCCTTGTTTTAAGTAAAGATAATAATCAGACCGGCAGCGGCAGACGCTGCTCAATAATACCGCCACCGAGGCACACTTCACCGTTATAGAAGACGGCAGACTGGCCCGGCGTCACGGCGGCAACCGGTTCATCGAAAATCACTTCAATGCGATCATCGTCCAGCGCCTTGACGGTGCAAGGGATATCGG  >  NC_000913/1192536‑1192830

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: