Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,578,548 G→A Q264* (CAG→TAG) 
Q264* (CAG→TAG) 
mcrB ←
mcrB ←
5‑methylcytosine‑specific restriction enzyme subunit McrB
5‑methylcytosine‑specific restriction enzyme subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,578,5480GA98.1% 178.3 / ‑3.6 53Q264* (CAG→TAG) 
Q264* (CAG→TAG) 
mcrB
mcrB
5‑methylcytosine‑specific restriction enzyme subunit McrB
5‑methylcytosine‑specific restriction enzyme subunit McrB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/24);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

CAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCG  >  NC_000913/4578311‑4578781
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
caacaaCGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCGTGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTaa                                                                                                                                                                                                                              <  1:508375/251‑1 (MQ=255)
caacaaCGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTaa                                                                                                                                                                                                                              <  1:493834/251‑1 (MQ=255)
    aacGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGACTCGTTCTTCATCGTTTTAGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACCAAAATTGTAAAat                                                                                                                                                                                                                           >  1:32925/1‑250 (MQ=255)
     acGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAatat                                                                                                                                                                                                                         >  1:297857/1‑251 (MQ=255)
           agagatCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTc                                                                                                                                                                                                                   <  2:120501/251‑1 (MQ=255)
                  gATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAACTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGt                                                                                                                                                                                                            >  1:762281/1‑251 (MQ=255)
                                          tgatATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACg                                                                                                                                                                                     >  1:342083/1‑250 (MQ=255)
                                            atatAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGAt                                                                                                                                                                                   >  1:807974/1‑250 (MQ=255)
                                             tataAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGAt                                                                                                                                                                                   >  1:68966/1‑249 (MQ=255)
                                                    aTTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGt                                                                                                                                                                           >  1:111318/1‑250 (MQ=255)
                                                       ctcCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCAGCCGACTCCATTCGGACGTTACCCCTGtata                                                                                                                                                                        >  2:114114/1‑250 (MQ=255)
                                                       ctcCGGGACATAGAATCGTTCTTCAACGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGATCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTAGGCCGATAGCCCTGtata                                                                                                                                                                        >  1:480622/1‑250 (MQ=255)
                                                           gggACATAGAATCGTTCTTCATCGTTTTCGGCGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAATCCGACGCCATTCGGACGATAGCCCTGTATAAAATc                                                                                                                                                                   >  1:178961/1‑251 (MQ=255)
                                                           gggACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATc                                                                                                                                                                   <  1:590922/251‑1 (MQ=255)
                                                               cATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTATAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCa                                                                                                                                                               <  2:155750/251‑1 (MQ=255)
                                                                aTAGAATCGATCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGCTTTCATCTATACTAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGAACGATAGCCCTGTATAAAATCCTCa                                                                                                                                                               >  2:429412/1‑250 (MQ=255)
                                                                             tcATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTATCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCAGCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAATCCGTAGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGa                                                                                                                                                    <  2:631830/248‑1 (MQ=255)
                                                                                tCGTTTTCGGAGGTGGTTAGGGTAACAGACCAGTTTTCACCTCGTGTATCATGTTCCATAAACATCATCACTTCGCCAAATACTTTACTGAGATTGTCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGACATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt                                                                                                                                               <  2:762281/250‑1 (MQ=255)
                                                                                    tttCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCCTTCGGCCGATAGCCCTGTATAAAATCCTCATCGCTATAAGATTGATTGAAc                                                                                                                                          >  1:811452/1‑251 (MQ=255)
                                                                                    tttCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                                                                                                                           <  2:111318/250‑1 (MQ=255)
                                                                                           gTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTAGAAAATCCTCAGATCTATAATATTGATGGAACTGAACCa                                                                                                                                   <  2:801233/251‑1 (MQ=255)
                                                                                            tAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCat                                                                                                                                  >  1:402749/1‑251 (MQ=255)
                                                                                                         cagtccAGTTTTCACCTCGTTTATCATGTTCCAATAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACTCGTTGc                                                                                                                     <  2:728256/247‑1 (MQ=255)
                                                                                                         cAGACCAGTTTTCAACTCGTTTATCATGTTCCATTTACTACAGCCCTTCGCCAAATTCTTTACTGTGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGATGCCATTCGGACGATAGCCCTTTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATGGACGCGTTGc                                                                                                                     <  2:83587/251‑1 (MQ=255)
                                                                                                          agtccAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                                                                                                                    <  1:278486/248‑1 (MQ=255)
                                                                                                          aGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                                                                                                                    <  1:319789/251‑1 (MQ=255)
                                                                                                                                            aaCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAATCCGACGCCATTCGGACGATAGCCCTGTATAAAATCATCATAGCTATAAGATTGATGGAAATGACACATATTTACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACg                                                                                  >  2:464517/1‑251 (MQ=255)
                                                                                                                                             aCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACgg                                                                                 <  1:12759/251‑1 (MQ=255)
                                                                                                                                                  atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGc                                                                            >  2:206029/1‑251 (MQ=255)
                                                                                                                                                    cccTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGcaa                                                                          <  2:63024/249‑1 (MQ=255)
                                                                                                                                                         cGCCAAATACTTTACTGAGATTGGCACGATTGATGTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAg                                                                     <  2:297857/251‑1 (MQ=255)
                                                                                                                                                            cAAATACTTTACTGAGATTGGCACGATTAATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTATTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGtt                                                                  >  2:442693/1‑251 (MQ=255)
                                                                                                                                                                             gtGGCACGATTTATTGTCACTATAATAAAAAAATACGTTTTCGCTGTCTGCTCCTTAGCTTGCTAACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGAAGCCATTCGGACGATAGCCCTGTATAAAAGCCTCATAGCTATAAGATTGTTGGAACTGAACCATATGGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGggcgg                                                 <  2:306182/250‑1 (MQ=255)
                                                                                                                                                                              tGGCCCGATTGATTTCATCTATAATAAAAATATTCTTTGTCTCTGGCTGCTCTTTAGCGTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGTCGTGCAACAAAGGTTTTTCCAACGCCGggcgg                                                 <  2:414129/250‑1 (MQ=255)
                                                                                                                                                                                gCACGATTGGTTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAACTTGCTAACAAAAATTGTAAATTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAACCTCATAGCTATAAGATTGATGGACCTGAACCATATTGACGCGTTGCGGCGCCTTTTCTCCTGTCAGCACGTAAGCCAGACGGCGTGCAACAAAGGCTTTACCAACGCCGGTCGGccc                                              >  2:379230/1‑251 (MQ=255)
                                                                                                                                                                                       tgattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTggag                                        <  1:723763/250‑1 (MQ=255)
                                                                                                                                                                                        gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTggagg                                       <  1:193828/250‑1 (MQ=255)
                                                                                                                                                                                          tttCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGAGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGat                                     >  2:755087/1‑250 (MQ=255)
                                                                                                                                                                                           ttCATCTATAATAAAAATATACTTTTTATCTGACTGCTCTTTAGCTTGCTAAAAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGata                                    >  1:632725/1‑250 (MQ=255)
                                                                                                                                                                                             tCATCTATAATACAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGAAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATACGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCACCAAAGGTTTTTCAAACGCCGGGCGGCCCCTGGGGGataat                                  >  2:602859/1‑251 (MQ=255)
                                                                                                                                                                                                tCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATAtt                               >  1:548140/1‑251 (MQ=255)
                                                                                                                                                                                                 cTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATAttt                              >  1:554318/1‑251 (MQ=255)
                                                                                                                                                                                                          aaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACTTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGAGCAACAAAGGTTTTTCCACAGCCGGGCGGCCCCTGGAGGATAATATTTATTTTgttg                     >  2:24638/1‑248 (MQ=255)
                                                                                                                                                                                                                aCTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTaa                <  1:434999/250‑1 (MQ=255)
                                                                                                                                                                                                                  tttttCTACGGCTGCTCTTTATCTTGCTAACAAAAATAGGAAAATATTCCGTCTTTACTTCGGAAGCCGACGGCATGCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCGTTTCTCCTGTCAGCAAGTAAGCCAGACGGCTTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTGTGATTGTTAATc              <  2:639744/250‑1 (MQ=255)
                                                                                                                                                                                                                   ttttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGt            >  1:873964/1‑251 (MQ=255)
                                                                                                                                                                                                                     ttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGttt          >  1:807783/1‑251 (MQ=255)
                                                                                                                                                                                                                        tctGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGGGCAAAAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTg         >  2:517782/1‑249 (MQ=255)
                                                                                                                                                                                                                         ctGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGt      >  2:350590/1‑251 (MQ=255)
                                                                                                                                                                                                                             ctgctCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGCGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCg  >  2:458907/1‑251 (MQ=255)
                                                                                                                                                                                                                               gctCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCg  <  1:458907/249‑1 (MQ=255)
                                                                                                                                                                                                                                      aGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTACAACGCCGggcggc                                                >  2:327692/1‑196 (MQ=255)
                                                                                                                                                                                                                                      aGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTACAACGCCGggcggc                                                <  1:327692/196‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
CAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCG  >  NC_000913/4578311‑4578781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: