Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,577,333 C→A 6.6% C209F (TGC→TTC)  mcrC ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,577,3330CA 6.6% 552.5 / 17.2 153C209F (TGC→TTC) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base C (74/68);  new base A (9/2);  total (83/70)
Fisher's exact test for biased strand distribution p-value = 6.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

CAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGC  >  NC_000913/4577190‑4577482
                                                                                                                                               |                                                                                                                                                     
cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAAtt                                                                                                                                                 <  1:1116115/150‑1 (MQ=255)
cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAAtt                                                                                                                                                 <  1:724888/150‑1 (MQ=255)
   cGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCAt                                                                                                                                              >  2:331109/1‑150 (MQ=255)
    gATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATa                                                                                                                                             >  2:1725616/1‑150 (MQ=255)
    gATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATa                                                                                                                                             >  1:2167990/1‑150 (MQ=255)
     aTAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATaa                                                                                                                                            >  1:1094048/1‑150 (MQ=255)
     aTAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATaa                                                                                                                                            >  1:1932908/1‑150 (MQ=255)
         gTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGaa                                                                                                                                        <  1:1856658/150‑1 (MQ=255)
               tttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACtttt                                                                                                                                  <  2:1214343/150‑1 (MQ=255)
                ttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTg                                                                                                                                 <  1:2325645/150‑1 (MQ=255)
                     atgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATaaa                                                                                                                            <  1:339131/150‑1 (MQ=255)
                      tgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAg                                                                                                                           <  2:2172117/150‑1 (MQ=255)
                       gaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGt                                                                                                                          >  1:140525/1‑150 (MQ=255)
                       gaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGt                                                                                                                          >  1:1310375/1‑150 (MQ=255)
                       gaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGt                                                                                                                          >  1:2109206/1‑150 (MQ=255)
                       gaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGt                                                                                                                          >  2:2588206/1‑150 (MQ=255)
                        aGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTa                                                                                                                         >  2:672193/1‑150 (MQ=255)
                        aGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTa                                                                                                                         >  2:162494/1‑150 (MQ=255)
                        aGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTa                                                                                                                         >  2:162565/1‑150 (MQ=255)
                        aGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTa                                                                                                                         >  2:455830/1‑150 (MQ=255)
                         gCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTaa                                                                                                                        >  1:1226491/1‑150 (MQ=255)
                          cGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAAt                                                                                                                       >  1:122382/1‑150 (MQ=255)
                            aaTGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGa                                                                                                                     <  1:1512883/150‑1 (MQ=255)
                                        aGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTc                                                                                                         >  1:2251960/1‑150 (MQ=255)
                                          tttCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGt                                                                                                       <  2:274083/150‑1 (MQ=255)
                                          tttCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGt                                                                                                       <  2:1176816/150‑1 (MQ=255)
                                                ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                                                                 <  2:1959607/150‑1 (MQ=255)
                                                ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                                                                 <  1:316879/150‑1 (MQ=255)
                                                ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                                                                 <  2:1930272/150‑1 (MQ=255)
                                                ttCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGtttcttt                                                                                                 <  2:2280550/150‑1 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  1:790363/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  1:2451684/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  1:1678854/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  1:2106445/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGAAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  2:2496959/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGAAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  2:1762291/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGAAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  2:20874/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGAAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  2:667948/1‑150 (MQ=255)
                                                      gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGAAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  2:1002537/1‑150 (MQ=255)
                                                      gTAACAAATTAAGAGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGAAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAAt                                                                                           >  2:2496907/1‑150 (MQ=255)
                                                         aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGTGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                                                        <  1:1428892/150‑1 (MQ=255)
                                                         aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                                                        <  1:34380/150‑1 (MQ=255)
                                                         aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                                                        <  2:622809/150‑1 (MQ=255)
                                                         aCAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCAt                                                                                        <  2:239430/150‑1 (MQ=255)
                                                          cAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATa                                                                                       >  1:1660322/1‑150 (MQ=255)
                                                                aaGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACg                                                                                 <  1:1780035/150‑1 (MQ=255)
                                                                   tGACTGATCCGATATACTCGATGGATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                                                              <  2:492263/150‑1 (MQ=255)
                                                                   tGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                                                              <  1:185152/150‑1 (MQ=255)
                                                                   tGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                                                              <  1:2695185/150‑1 (MQ=255)
                                                                       tGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgtgt                                                                          >  1:545386/1‑150 (MQ=255)
                                                                       tGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgtgt                                                                          >  1:792598/1‑150 (MQ=255)
                                                                       tGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgtgt                                                                          >  1:738746/1‑150 (MQ=255)
                                                                        gATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTc                                                                         <  2:400078/150‑1 (MQ=255)
                                                                              atatACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGt                                                                   >  2:249021/1‑150 (MQ=255)
                                                                              atatACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGt                                                                   >  2:1278621/1‑150 (MQ=255)
                                                                                ataCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGttt                                                                 <  2:2701925/150‑1 (MQ=255)
                                                                                ataCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGttt                                                                 <  2:2701933/150‑1 (MQ=255)
                                                                                ataCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGttt                                                                 <  2:781702/150‑1 (MQ=255)
                                                                                 taCTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGtttt                                                                <  1:575578/150‑1 (MQ=255)
                                                                                   cTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGa                                                              <  2:608212/150‑1 (MQ=255)
                                                                                       atccctCCCATATTAAATACGACCGGGTGGTGCTCGCAGACGTTAATTCACGACGGCAAAATCCATAAAGAAACATTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTg                                                          <  2:2170014/145‑1 (MQ=255)
                                                                                              ccATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGa                                                   >  2:2546734/1‑150 (MQ=255)
                                                                                                aTTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAat                                                 <  2:1355360/150‑1 (MQ=255)
                                                                                                aTTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAat                                                 >  2:2405913/1‑150 (MQ=255)
                                                                                                 ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAatt                                                >  2:2302789/1‑150 (MQ=255)
                                                                                                 ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAatt                                                >  1:190463/1‑150 (MQ=255)
                                                                                                 ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAatt                                                >  1:2217393/1‑150 (MQ=255)
                                                                                                 ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAatt                                                >  1:2217398/1‑150 (MQ=255)
                                                                                                 ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAatt                                                <  1:331109/150‑1 (MQ=255)
                                                                                                 ttttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATAATAGAAACGGTAGTGTCCTTTGTTTTGTCCTGGAATAGAatt                                                >  1:2191437/1‑150 (MQ=255)
                                                                                                   ttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAattat                                              <  1:1083589/150‑1 (MQ=255)
                                                                                                    tAAATAAGAGCGTGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAattatt                                             <  2:2567797/150‑1 (MQ=255)
                                                                                                    tAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTGTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAattatt                                             <  2:1752896/150‑1 (MQ=255)
                                                                                                       aTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGAc                                          >  2:1935521/1‑150 (MQ=255)
                                                                                                       aTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGAc                                          <  1:1088990/150‑1 (MQ=255)
                                                                                                         aaGAGCGGGTTGTGTTGGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGa                                        <  1:634951/150‑1 (MQ=255)
                                                                                                          agagCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGAt                                       <  2:1249527/150‑1 (MQ=255)
                                                                                                            cgCGTGTTGTGTTTGGAGACTTTAATTCAAGTCGTAAAAATTCATAAAGTATGTTTTGATAATGTAATGACTTCTCTTTTTCGTTTCTTTCAATATCTTAGAAAAGGTATTGTCCTTTGTTTTGACCTGGAATTGAATTATTGACGATgt                                     <  1:1694527/149‑2 (MQ=255)
                                                                                                            agCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGa                                     <  2:546515/150‑1 (MQ=255)
                                                                                                            agCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGa                                     <  2:834152/150‑1 (MQ=255)
                                                                                                                       tttGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                          >  2:912143/1‑150 (MQ=255)
                                                                                                                       tttGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                          >  2:320428/1‑150 (MQ=255)
                                                                                                                       tttGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                          >  2:226756/1‑150 (MQ=255)
                                                                                                                       tttGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                          >  2:1475106/1‑150 (MQ=255)
                                                                                                                       tttGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                          >  2:1032961/1‑150 (MQ=255)
                                                                                                                       tttGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGAAGaaa                          >  1:2571770/1‑148 (MQ=255)
                                                                                                                       tttGCAGACGTTAATACACGACGCAAAAATTCAAAAACAAACTATTGATAAACTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCATTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                          <  1:2208671/150‑1 (MQ=255)
                                                                                                                           cAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGa                      <  2:212981/150‑1 (MQ=255)
                                                                                                                           cAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGa                      <  1:2624029/150‑1 (MQ=255)
                                                                                                                           cAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGa                      <  1:2406663/150‑1 (MQ=255)
                                                                                                                                gTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTCGTGTCCTTTGGTGTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACCCTGATAACg                 >  2:344350/1‑150 (MQ=255)
                                                                                                                                 ttAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGa                >  2:118130/1‑150 (MQ=255)
                                                                                                                                 ttAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGa                >  1:928948/1‑150 (MQ=255)
                                                                                                                                 ttAAGTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGa                <  1:530284/150‑1 (MQ=255)
                                                                                                                                        acgacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTtata         >  2:2585129/1‑150 (MQ=255)
                                                                                                                                        acgacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTtata         <  1:403848/150‑1 (MQ=255)
                                                                                                                                          gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTataat       >  2:1636756/1‑150 (MQ=255)
                                                                                                                                               cAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAAcgc  >  1:1485723/1‑150 (MQ=255)
                                                                                                                                               cAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAAcgc  >  1:1485736/1‑150 (MQ=255)
                                                                                                                                               cAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAAcgc  >  1:531836/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGC  >  NC_000913/4577190‑4577482

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