Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 604,752 G→T 6.2% intergenic (‑89/+19) ybdG ← / ← nfsB miniconductance mechanosensitive channel YbdG/NAD(P)H nitroreductase NfsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913604,7520GT 6.2% 275.9 / 4.9 81intergenic (‑89/+19)ybdG/nfsBminiconductance mechanosensitive channel YbdG/NAD(P)H nitroreductase NfsB
Reads supporting (aligned to +/- strand):  ref base G (47/29);  new base T (2/3);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATT  >  NC_000913/604603‑604901
                                                                                                                                                     |                                                                                                                                                     
aGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAg                                                                                                                                                       <  1:335743/150‑1 (MQ=255)
 ggtggtGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGa                                                                                                                                                      <  1:1003094/150‑1 (MQ=255)
    ggtgTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGc                                                                                                                                                   >  2:638755/1‑150 (MQ=255)
       gtgATCGATCTCAATACGCTCTAAATCTTCAACGTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCg                                                                                                                                                <  1:378239/150‑1 (MQ=255)
       gtgATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCg                                                                                                                                                <  2:1780026/150‑1 (MQ=255)
         gatcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGc                                                                                                                                              >  1:1505731/1‑150 (MQ=255)
           tcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGcta                                                                                                                                            <  2:405822/150‑3 (MQ=255)
                ctcAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGctatttcc                                                                                                                                       <  2:570550/150‑8 (MQ=255)
                  cAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGaa                                                                                                                                     <  1:638755/150‑1 (MQ=255)
                    aTACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAAtt                                                                                                                                   <  1:1514621/150‑1 (MQ=255)
                       cccGCTTAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                <  1:328301/150‑1 (MQ=255)
                       cccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                >  1:457138/1‑150 (MQ=255)
                       cccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                >  1:564869/1‑150 (MQ=255)
                       cccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                >  1:1727116/1‑150 (MQ=255)
                       cccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                <  1:1648124/150‑1 (MQ=255)
                                 cTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                      >  2:401404/1‑150 (MQ=255)
                                  ttCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAg                                                                                                                     >  1:1990656/1‑150 (MQ=255)
                                              tatTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGc                                                                                                         >  1:1264862/1‑150 (MQ=255)
                                                ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCgg                                                                                                       <  1:807548/150‑1 (MQ=255)
                                                 tAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAATGTGATGTTTTGCGGc                                                                                                      <  2:564869/150‑1 (MQ=255)
                                                    aTCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGa                                                                                                   >  2:1025598/1‑150 (MQ=255)
                                                      ccTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACCCTTCGGTTAAGGTGATGTTTTGCGGCAGACg                                                                                                 >  2:28787/1‑150 (MQ=255)
                                                      ccTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACg                                                                                                 >  1:241327/1‑150 (MQ=255)
                                                       cTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACga                                                                                                >  2:1137500/1‑150 (MQ=255)
                                                        tGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACgag                                                                                               >  1:1135809/1‑150 (MQ=255)
                                                          cATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGAt                                                                                             >  1:1085479/1‑150 (MQ=255)
                                                            tAACGTCTCCTTTAAACAACAGGAGCCTATGATGACGGCTGAAACAGGGTTATGCAACTCAGGAGAATCTGAGAGGAAATACCCTGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGAttt                                                                                           <  2:1409360/150‑1 (MQ=255)
                                                              aCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCg                                                                                         <  2:959154/150‑1 (MQ=255)
                                                                  ctcCTTTATAGAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCag                                                                                     >  2:678310/1‑150 (MQ=255)
                                                                    ccTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcg                                                                                   >  2:1331692/1‑150 (MQ=255)
                                                                    ccTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcg                                                                                   <  2:1505731/150‑1 (MQ=255)
                                                                        tataCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCag                                                                                     >  2:1366181/1‑144 (MQ=255)
                                                                        tataCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCag                                                                                     <  1:1366181/144‑1 (MQ=255)
                                                                         ataCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgtagcg                                                                              >  1:1935192/1‑150 (MQ=255)
                                                                         ataCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgtagcg                                                                              >  2:1446548/1‑150 (MQ=255)
                                                                          taCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgtagcgt                                                                             >  1:441561/1‑150 (MQ=255)
                                                                           acaacaGCAGCCTATGATGACGGCGGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGtt                                                                            <  1:1137500/150‑1 (MQ=255)
                                                                            caacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                           >  1:99224/1‑150 (MQ=255)
                                                                               cagcagCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTaaaa                                                                        >  1:175716/1‑150 (MQ=255)
                                                                                    cccTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTc                                                                   <  2:513140/149‑1 (MQ=255)
                                                                                     ccTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCa                                                                  >  1:922019/1‑150 (MQ=255)
                                                                                         tgatgaCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGc                                                              >  2:864879/1‑150 (MQ=255)
                                                                                          gatgaCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCt                                                             >  2:1253524/1‑150 (MQ=255)
                                                                                              aCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtga                                                         <  2:508335/150‑1 (MQ=255)
                                                                                                ggCTGAAACAGGGTTATGCAATTGAGGAGAATCTGATAGGAAATTGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGTCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgatg                                                       <  1:1835497/150‑1 (MQ=255)
                                                                                                  cTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgatgag                                                     >  2:1796304/1‑149 (MQ=255)
                                                                                                     aaaCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTa                                                  >  2:1229543/1‑150 (MQ=255)
                                                                                                     aaaCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTa                                                  >  2:700045/1‑150 (MQ=255)
                                                                                                       aCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTAcc                                                <  1:1025598/150‑1 (MQ=255)
                                                                                                       aCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTAcc                                                >  1:650765/1‑150 (MQ=255)
                                                                                                            gTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGaac                                           >  2:1494494/1‑150 (MQ=255)
                                                                                                              tATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGaacaa                                         >  1:1703778/1‑150 (MQ=255)
                                                                                                              tATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGaacaa                                         >  1:1909735/1‑150 (MQ=255)
                                                                                                                tGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAacc                                       >  2:1959356/1‑150 (MQ=255)
                                                                                                                tGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAacc                                       >  1:1228802/1‑150 (MQ=255)
                                                                                                                  cAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATTTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                     <  1:678310/150‑1 (MQ=255)
                                                                                                                   aaaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccacc                                    <  2:1135809/150‑1 (MQ=255)
                                                                                                                    aaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAAccacca                                   >  1:1476366/1‑150 (MQ=255)
                                                                                                                     aTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAg                                  <  2:441561/150‑1 (MQ=255)
                                                                                                                          gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTGAACGATGTGATGACCTACCGGAACAACCACCAGACTgg                             >  2:872902/1‑150 (MQ=255)
                                                                                                                          gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTAAACGCTGTGATGACCTACCGGAACAACCAGCAGACTgg                             >  2:866703/1‑150 (MQ=255)
                                                                                                                              aTCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTa                         >  2:287224/1‑150 (MQ=255)
                                                                                                                              aTCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTa                         >  1:240513/1‑150 (MQ=255)
                                                                                                                                cTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGc                       <  1:1446548/150‑1 (MQ=255)
                                                                                                                                  gagaGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCt                     >  2:1966562/1‑150 (MQ=255)
                                                                                                                                  gagaGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCt                     >  2:281411/1‑150 (MQ=255)
                                                                                                                                   agagGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCtt                    >  1:611422/1‑150 (MQ=255)
                                                                                                                                    gagGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCttt                   <  2:1228802/150‑1 (MQ=255)
                                                                                                                                     agGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtc                  <  1:815465/150‑1 (MQ=255)
                                                                                                                                      ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtct                 >  2:1488741/1‑150 (MQ=255)
                                                                                                                                      ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtct                 >  2:113474/1‑150 (MQ=255)
                                                                                                                                      ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtct                 >  2:457510/1‑150 (MQ=255)
                                                                                                                                      ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtct                 >  2:1435103/1‑150 (MQ=255)
                                                                                                                                          cttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCttt             >  2:1594138/4‑150 (MQ=255)
                                                                                                                                           tAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTc            <  2:457138/150‑1 (MQ=255)
                                                                                                                                           tAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTc            <  2:241327/150‑1 (MQ=255)
                                                                                                                                           tAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTc            <  2:922019/150‑1 (MQ=255)
                                                                                                                                            tgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCa           >  2:422475/2‑150 (MQ=255)
                                                                                                                                            aGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCa           >  2:239382/1‑150 (MQ=255)
                                                                                                                                             gCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAg          <  2:240513/150‑1 (MQ=255)
                                                                                                                                             gCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAg          <  2:1990656/150‑1 (MQ=255)
                                                                                                                                               cGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGAc        <  1:422475/150‑1 (MQ=255)
                                                                                                                                                     gATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAAtt  >  1:406241/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATT  >  NC_000913/604603‑604901

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: