Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,579,468 G→T intergenic (‑131/+31) mcrB ← / ← symE 5‑methylcytosine‑specific restriction enzyme subunit McrB/toxic protein SymE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,4680GT100.0% 602.1 / NA 167intergenic (‑131/+31)mcrB/symE5‑methylcytosine‑specific restriction enzyme subunit McrB/toxic protein SymE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (103/64);  total (103/64)

Only 100 of 170 total aligned reads displayed.

TTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCG  >  NC_000913/4579332‑4579604
                                                                                                                                        |                                                                                                                                        
ttCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGtctc                                                                                                                                       >  1:1021801/1‑140 (MQ=255)
 tCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTcc                                                                                                                                      <  2:866521/140‑1 (MQ=255)
  ccATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTccc                                                                                                                                     <  1:100155/140‑1 (MQ=255)
    aTAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCgg                                                                                                                                   >  1:2369695/1‑140 (MQ=255)
    aTAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCgg                                                                                                                                   >  2:634681/1‑140 (MQ=255)
         atatGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGtt                                                                                                                              >  1:2781291/1‑140 (MQ=255)
         atatGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGtt                                                                                                                              >  1:2173674/1‑140 (MQ=255)
             gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAAcc                                                                                                                          >  1:1984665/1‑140 (MQ=255)
             gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAAcc                                                                                                                          >  1:1201735/1‑140 (MQ=255)
             gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAAcc                                                                                                                          >  1:304327/1‑140 (MQ=255)
             gCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAAcc                                                                                                                          <  1:1645121/140‑1 (MQ=255)
                 tCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGa                                                                                                                      >  1:2137133/1‑140 (MQ=255)
                  ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCggg                                                                                                                       <  1:2493434/138‑1 (MQ=255)
                  ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCggg                                                                                                                       >  2:2493434/1‑138 (MQ=255)
                  ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAg                                                                                                                     >  1:1785682/1‑140 (MQ=255)
                  ccTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAg                                                                                                                     >  2:1419806/1‑140 (MQ=255)
                   cTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGt                                                                                                                    <  2:2173674/140‑1 (MQ=255)
                      aTTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTaaa                                                                                                                 <  1:1419170/140‑1 (MQ=255)
                        tGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACg                                                                                                               >  1:1046357/1‑140 (MQ=255)
                        tGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACg                                                                                                               >  1:1319771/1‑140 (MQ=255)
                         gCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGc                                                                                                              >  1:841550/1‑140 (MQ=255)
                             aGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGtt                                                                                                          >  1:2463757/1‑140 (MQ=255)
                               ccTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAc                                                                                                        <  1:821793/140‑1 (MQ=255)
                                  ctATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAcgcg                                                                                                     >  2:2031841/1‑140 (MQ=255)
                                  ctATTATCGCTCTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAcgcg                                                                                                     >  2:2133289/1‑140 (MQ=255)
                                    attatCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGAc                                                                                                   >  2:802604/1‑140 (MQ=255)
                                            tttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGtt                                                                                           <  2:75163/140‑1 (MQ=255)
                                             ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTAcg                                                                                                       >  1:1928776/1‑127 (MQ=255)
                                             ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                                                                          >  1:426506/1‑140 (MQ=255)
                                             ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                                                                          >  1:430234/1‑140 (MQ=255)
                                             ttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                                                                          >  2:1355573/1‑140 (MQ=255)
                                               cGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTAc                                                                                        >  2:700147/1‑140 (MQ=255)
                                                gCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTAcc                                                                                       <  1:674804/140‑1 (MQ=255)
                                                    cGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCa                                                                                   <  1:700147/140‑1 (MQ=255)
                                                     gTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCaa                                                                                  <  1:2133289/140‑1 (MQ=255)
                                                      tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                                                                                 <  1:286510/140‑1 (MQ=255)
                                                      tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                                                                                 <  2:1021801/140‑1 (MQ=255)
                                                      tACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                                                                                 <  1:782335/140‑1 (MQ=255)
                                                       aCTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                                                                                <  2:1984665/140‑1 (MQ=255)
                                                        cTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCa                                                                               <  2:1785682/140‑1 (MQ=255)
                                                        cTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCa                                                                               <  2:229827/140‑1 (MQ=255)
                                                                aTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa                                                                       >  1:2593786/1‑140 (MQ=255)
                                                                aTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa                                                                       >  2:578518/1‑140 (MQ=255)
                                                                aTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATa                                                                       >  1:1293440/1‑140 (MQ=255)
                                                                 tAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCa                                                                               <  2:608083/131‑1 (MQ=255)
                                                                 tAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTc                                                                            >  1:1383049/1‑134 (MQ=255)
                                                                 tAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                                      >  1:1088600/1‑140 (MQ=255)
                                                                 tAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                                      >  1:1557777/1‑140 (MQ=255)
                                                                         tACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                              <  1:609916/140‑1 (MQ=255)
                                                                         tACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                              >  2:311170/1‑140 (MQ=255)
                                                                             gCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCt                                                          >  2:2043917/1‑140 (MQ=255)
                                                                             gCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCt                                                          >  2:1447085/1‑140 (MQ=255)
                                                                              cAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTg                                                         >  1:2582217/1‑140 (MQ=255)
                                                                                      gACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGtt                                                 >  2:227424/1‑140 (MQ=255)
                                                                                        cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                               <  1:1419806/140‑1 (MQ=255)
                                                                                         tGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAc                                              <  1:2263374/140‑1 (MQ=255)
                                                                                              aaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGcc                                         <  1:802604/140‑1 (MQ=255)
                                                                                               atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                                        <  1:1355573/140‑1 (MQ=255)
                                                                                               atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                                        <  2:1372174/140‑1 (MQ=255)
                                                                                                 atTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGAc                                      >  1:349616/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  2:488844/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  1:647565/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  2:2337086/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  2:2375398/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  2:1769795/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  2:649170/1‑140 (MQ=255)
                                                                                                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                                  >  2:1238623/1‑140 (MQ=255)
                                                                                                      cTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGttt                                 >  1:922307/1‑140 (MQ=255)
                                                                                                          gggggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcac                             <  1:2375398/140‑1 (MQ=255)
                                                                                                          gggggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcac                             >  1:565702/1‑140 (MQ=255)
                                                                                                           ggggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca                            <  2:2502529/140‑1 (MQ=255)
                                                                                                           ggggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca                            >  1:57883/1‑140 (MQ=255)
                                                                                                           ggggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca                            >  1:1137100/1‑140 (MQ=255)
                                                                                                           ggggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGCTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca                            >  2:1459295/1‑140 (MQ=255)
                                                                                                            ggtaaaGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCTCTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCacac                           <  2:2781291/137‑1 (MQ=255)
                                                                                                            gggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCacac                           >  2:2366894/1‑140 (MQ=255)
                                                                                                             ggAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACAcc                          <  1:2031841/140‑1 (MQ=255)
                                                                                                              gAAAGTTTCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCt                         <  1:578518/140‑1 (MQ=255)
                                                                                                                   tttCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGt                    <  2:565702/140‑1 (MQ=255)
                                                                                                                   tttCTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGt                    >  2:2738777/1‑140 (MQ=255)
                                                                                                                      cTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGc                 >  2:541282/1‑140 (MQ=255)
                                                                                                                      cTATTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGc                 >  2:1061586/1‑140 (MQ=255)
                                                                                                                        aTTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGa               <  1:2729985/140‑1 (MQ=255)
                                                                                                                        aTTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGa               >  1:2388083/1‑140 (MQ=255)
                                                                                                                        aTTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGa               <  1:2331234/140‑1 (MQ=255)
                                                                                                                        aTTGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGa               >  2:2331234/1‑140 (MQ=255)
                                                                                                                         ttGAGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGAc              <  1:15407/140‑1 (MQ=255)
                                                                                                                            aGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGc           >  1:2287164/1‑140 (MQ=255)
                                                                                                                            aGTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCACCGACAGTTTGCACACCTGGCGTAGCGACTGc           >  2:743842/1‑140 (MQ=255)
                                                                                                                             gTCAGTGGAAGTGTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCa          >  2:157140/1‑140 (MQ=255)
                                                                                                                             gTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCa          <  1:2302979/140‑1 (MQ=255)
                                                                                                                             gTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCa          >  1:742643/1‑140 (MQ=255)
                                                                                                                             gTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCa          <  2:426506/140‑1 (MQ=255)
                                                                                                                             gTCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCa          >  2:2693974/1‑140 (MQ=255)
                                                                                                                             gTCAGTGGAAGTCTCCCGGTGGCTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCa          <  1:1459295/140‑1 (MQ=255)
                                                                                                                              tCAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCAt         >  2:1116378/1‑140 (MQ=255)
                                                                                                                               cAGTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATc        <  2:293883/140‑1 (MQ=255)
                                                                                                                                 gTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAg      >  1:2257877/1‑140 (MQ=255)
                                                                                                                                 gTGGAAGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAg      >  1:900327/1‑140 (MQ=255)
                                                                                                                                     aaGTCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCg  <  2:2388083/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCG  >  NC_000913/4579332‑4579604

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