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

CAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCG  >  NC_000913/4578311‑4578781
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
caacaaCGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCGTGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTaa                                                                                                                                                                                                                              <  1:508375/251‑1 (MQ=255)
caacaaCGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTaa                                                                                                                                                                                                                              <  1:493834/251‑1 (MQ=255)
    aacGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGACTCGTTCTTCATCGTTTTAGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACCAAAATTGTAAAat                                                                                                                                                                                                                           >  1:32925/1‑250 (MQ=255)
     acGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAatat                                                                                                                                                                                                                         >  1:297857/1‑251 (MQ=255)
           agagatCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTc                                                                                                                                                                                                                   <  2:120501/251‑1 (MQ=255)
                  gATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAACTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGt                                                                                                                                                                                                            >  1:762281/1‑251 (MQ=255)
                                          tgatATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACg                                                                                                                                                                                     >  1:342083/1‑250 (MQ=255)
                                            atatAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGAt                                                                                                                                                                                   >  1:807974/1‑250 (MQ=255)
                                             tataAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGAt                                                                                                                                                                                   >  1:68966/1‑249 (MQ=255)
                                                    aTTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGt                                                                                                                                                                           >  1:111318/1‑250 (MQ=255)
                                                       ctcCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCAGCCGACTCCATTCGGACGTTACCCCTGtata                                                                                                                                                                        >  2:114114/1‑250 (MQ=255)
                                                       ctcCGGGACATAGAATCGTTCTTCAACGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGATCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTAGGCCGATAGCCCTGtata                                                                                                                                                                        >  1:480622/1‑250 (MQ=255)
                                                           gggACATAGAATCGTTCTTCATCGTTTTCGGCGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAATCCGACGCCATTCGGACGATAGCCCTGTATAAAATc                                                                                                                                                                   >  1:178961/1‑251 (MQ=255)
                                                           gggACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATc                                                                                                                                                                   <  1:590922/251‑1 (MQ=255)
                                                               cATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTATAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCa                                                                                                                                                               <  2:155750/251‑1 (MQ=255)
                                                                aTAGAATCGATCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGCTT‑TCATCTATACTAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGAACGATAGCCCTGTATAAAATCCTCa                                                                                                                                                               >  2:429412/1‑250 (MQ=255)
                                                                             tcATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTATCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCAGCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAATCCGTAGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGa                                                                                                                                                    <  2:631830/248‑1 (MQ=255)
                                                                                tCGTTTTCGGAGGTGGTTAGGGTAACAGACCAGTTTTCACCTCGTGTATCATGTTCCATAAACATCATCACTTCGCCAAATACTTTACTGAGATTGTCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGACATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt                                                                                                                                               <  2:762281/250‑1 (MQ=255)
                                                                                    tttCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCCTTCGGCCGATAGCCCTGTATAAAATCCTCATCGCTATAAGATTGATTGAAc                                                                                                                                          >  1:811452/1‑251 (MQ=255)
                                                                                    tttCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                                                                                                                           <  2:111318/250‑1 (MQ=255)
                                                                                           gTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTAGAAAATCCTCAGATCTATAATATTGATGGAACTGAACCa                                                                                                                                   <  2:801233/251‑1 (MQ=255)
                                                                                            tAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCat                                                                                                                                  >  1:402749/1‑251 (MQ=255)
                                                                                                         cagtccAGTTTTCACCTCGTTTATCATGTTCCAATAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACTCGTTGc                                                                                                                     <  2:728256/247‑1 (MQ=255)
                                                                                                         cAGACCAGTTTTCAACTCGTTTATCATGTTCCATTTACTACAGCCCTTCGCCAAATTCTTTACTGTGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGATGCCATTCGGACGATAGCCCTTTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATGGACGCGTTGc                                                                                                                     <  2:83587/251‑1 (MQ=255)
                                                                                                          agtccAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                                                                                                                    <  1:278486/248‑1 (MQ=255)
                                                                                                          aGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                                                                                                                    <  1:319789/251‑1 (MQ=255)
                                                                                                                                            aaCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAATCCGACGCCATTCGGACGATAGCCCTGTATAAAATCATCATAGCTATAAGATTGATGGAAATGACACATATTTACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACg                                                                                  >  2:464517/1‑251 (MQ=255)
                                                                                                                                             aCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACgg                                                                                 <  1:12759/251‑1 (MQ=255)
                                                                                                                                                  atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGc                                                                            >  2:206029/1‑251 (MQ=255)
                                                                                                                                                    cccTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGcaa                                                                          <  2:63024/249‑1 (MQ=255)
                                                                                                                                                         cGCCAAATACTTTACTGAGATTGGCACGATTGATG‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAg                                                                     <  2:297857/251‑1 (MQ=255)
                                                                                                                                                            cAAATACTTTACTGAGATTGGCACGATTAATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTATTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGtt                                                                  >  2:442693/1‑251 (MQ=255)
                                                                                                                                                                             gtGGCACGATTTATTGTCA‑CTATAATAAAAAAATACGTTTTCGCTGTCTGCTCCTTAGCTTGCTAACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGAAGCCATTCGGACGATAGCCCTGTATAAAAGCCTCATAGCTATAAGATTGTTGGAACTGAACCATATGGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGggcgg                                                 <  2:306182/250‑1 (MQ=255)
                                                                                                                                                                              tGGCCCGATTGATT‑TCATCTATAATAAAAATATTCTTTGTCTCTGGCTGCTCTTTAGCGTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGTCGTGCAACAAAGGTTTTTCCAACGCCGggcgg                                                 <  2:414129/250‑1 (MQ=255)
                                                                                                                                                                                gCACGATTGGTT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAACTTGCTAACAAAAATTGTAAATTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAACCTCATAGCTATAAGATTGATGGACCTGAACCATATTGACGCGTTGCGGCGCCTTTTCTCCTGTCAGCACGTAAGCCAGACGGCGTGCAACAAAGGCTTTACCAACGCCGGTCGGccc                                              >  2:379230/1‑251 (MQ=255)
                                                                                                                                                                                       tgatt‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTggag                                        <  1:723763/250‑1 (MQ=255)
                                                                                                                                                                                        gatt‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTggagg                                       <  1:193828/250‑1 (MQ=255)
                                                                                                                                                                                          tt‑tCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGAGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGat                                     >  2:755087/1‑250 (MQ=255)
                                                                                                                                                                                           t‑tCATCTATAATAAAAATATACTTTTTATCTGACTGCTCTTTAGCTTGCTAAAAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGata                                    >  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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
CAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATT‑TCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCG  >  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: ‑