Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,718,318 T→G 7.8% T71P (ACG→CCG)  yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,718,3180TG 7.8% 1714.4 / 25.8 582T71P (ACG→CCG) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (348/176);  new base G (2/56);  total (350/232)
Fisher's exact test for biased strand distribution p-value = 3.52e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 598 total aligned reads displayed.

GAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCAT  >  NC_000913/2718169‑2718467
                                                                                                                                                     |                                                                                                                                                     
gAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGACCCCGt                                                                                                                                                       >  1:3457672/1‑150 (MQ=255)
gAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGt                                                                                                                                                       <  1:579969/150‑1 (MQ=255)
 aaaaGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTg                                                                                                                                                      >  2:3577346/1‑150 (MQ=255)
        ccTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACg                                                                                                                                               >  2:10270230/1‑150 (MQ=255)
              aTGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAAt                                                                                                                                         >  1:5700205/1‑150 (MQ=255)
               tGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGGTGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATc                                                                                                                                        >  2:8506630/1‑150 (MQ=255)
               tGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATc                                                                                                                                        >  1:4525099/1‑150 (MQ=255)
               tGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATc                                                                                                                                        >  1:9256304/1‑150 (MQ=255)
               tGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATc                                                                                                                                        >  2:5456182/1‑150 (MQ=255)
                gCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCt                                                                                                                                       >  2:4453092/1‑150 (MQ=255)
                  aTTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTc                                                                                                                                     >  1:2722252/1‑150 (MQ=255)
                   ttttCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTcc                                                                                                                                    >  1:6459983/1‑150 (MQ=255)
                   ttggcGCCTTAGACGCGGTTTTCTTCCGCACGCTGCCGACCCAACACTCCCGGATCAATAAAAGTTTTACCCCTGCTGCCACCGTGATCGGCCTTTGCCGGCGACGCACCCCCGGGCGCGCGGTAACCCGGGCGCCACGGCGAATCTTcc                                                                                                                                    <  2:7560947/146‑1 (MQ=255)
                         ccGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGACTCTTCCCACTTg                                                                                                                              >  2:201547/1‑150 (MQ=255)
                         ccGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTg                                                                                                                              >  2:10417638/1‑150 (MQ=255)
                          cGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGACCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGc                                                                                                                             >  1:5553198/1‑150 (MQ=255)
                           gTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCg                                                                                                                            >  2:4758392/1‑150 (MQ=255)
                            tAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCgg                                                                                                                           >  1:2386633/1‑150 (MQ=255)
                               aCGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTcc                                                                                                                        >  1:5449994/1‑149 (MQ=255)
                                cgcgGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCAc                                                                                                                       >  1:3952459/1‑150 (MQ=255)
                                cgcgGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCAc                                                                                                                       >  2:701928/1‑150 (MQ=255)
                                cgcgGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCAc                                                                                                                       >  2:9071305/1‑150 (MQ=255)
                                 gcgGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACg                                                                                                                      >  1:1628690/1‑150 (MQ=255)
                                    ggTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACgggg                                                                                                                   >  2:7837022/1‑150 (MQ=255)
                                     gTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGc                                                                                                                  >  2:10891127/1‑150 (MQ=255)
                                              gaaCGCTCACCACCAAACACTTCGGCAAAAAAAAAATTTTAACCACTTAGGCCACCTTTATCGGCATTTTCCGGCTTAACACCACCGGGCGCGGGTCAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGt                                                                                                         <  2:4889456/148‑1 (MQ=255)
                                                aCGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGttt                                                                                                       <  2:10028496/150‑1 (MQ=255)
                                                        gacgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATcgc                                                                                               >  2:4173039/1‑150 (MQ=255)
                                                        gacgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATcgc                                                                                               >  2:7939080/1‑150 (MQ=255)
                                                        gacgCAACACTTCCGGATCAATAAAACTTTTACCACAGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTACCCACTAGCGGTCACGGGGCTGACTCTGTTTGTCATcgc                                                                                               >  2:6223092/1‑150 (MQ=255)
                                                         ccgCAACACTCCCGTAAAAATAAAACTTTCACCCCTTATGCCACGGTGTTCGGCCTTTTCGGGCGTCACACCCCCGGGCGCGGGGTAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATcgcg                                                                                              <  1:1837587/149‑1 (MQ=255)
                                                          cgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCgcgg                                                                                             >  1:7098566/1‑150 (MQ=255)
                                                            cAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCgcgggc                                                                                           >  1:7669859/1‑150 (MQ=255)
                                                               cccTTCCGGATAAATAAAACTTTTACCACGGATGCCACCGTGATGGGCCTTTTCCGGCGTCACATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGcgggcggg                                                                                        <  2:5859780/148‑1 (MQ=255)
                                                               cacTTCCGCATAAATAAAACTTTTACCACGGATGCCACGGTGTTCGGCATTTTCCGGCGTAACAACCCCGGGCGCGCGGCAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGcgggcggg                                                                                        <  2:5981519/150‑1 (MQ=255)
                                                                  ttCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCgggcggggcg                                                                                     >  1:7322607/1‑150 (MQ=255)
                                                                   tCCGGATCAATAAAACTTTTACCACTGATGCCACGGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGa                                                                                    <  2:1623880/150‑1 (MQ=255)
                                                                    ccGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGAc                                                                                   >  1:1572079/1‑150 (MQ=255)
                                                                      ggATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACgg                                                                                 >  2:1164111/1‑150 (MQ=255)
                                                                          cAATAAAACTTTTACCACTGAGGCCACCGTGATCGGCCTTTGCCGGCGTCTCACCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCAc                                                                             <  2:9812773/150‑1 (MQ=255)
                                                                           aaTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCAcc                                                                            >  1:10677608/1‑150 (MQ=255)
                                                                            aTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGTCGGCCACCg                                                                           >  1:6270474/1‑150 (MQ=255)
                                                                             tAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGt                                                                          >  1:3287149/1‑150 (MQ=255)
                                                                             tAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGt                                                                          >  2:3145392/1‑150 (MQ=255)
                                                                              aaaaCTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTc                                                                         >  2:3788240/1‑150 (MQ=255)
                                                                                   ttttCCCACTGAGGCCACCGTGATCGGCATTTTCCGGCGTCTCACCACCCGGCGCGGGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCAc                                                                    <  2:3932078/150‑1 (MQ=255)
                                                                                    tttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACg                                                                   >  1:472599/1‑150 (MQ=255)
                                                                                     ttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACga                                                                  >  1:117895/1‑150 (MQ=255)
                                                                                     ttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACga                                                                  >  2:10580133/1‑150 (MQ=255)
                                                                                     ttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACga                                                                  >  2:9218294/1‑150 (MQ=255)
                                                                                     tcaCCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCACATCACCCGGCGCGGGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACga                                                                  <  1:895581/148‑1 (MQ=255)
                                                                                        ccACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGAt                                                               >  1:8137809/1‑150 (MQ=255)
                                                                                        ccACTGATGCCACCATGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGAt                                                               >  1:7008116/1‑150 (MQ=255)
                                                                                          ccTGAGGCCCCCGTGTGGGGCGTTTCCGGGCGCCTCACCCCCGGGCGCGGGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGAttt                                                             <  1:5022729/149‑1 (MQ=255)
                                                                                          aCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGAttt                                                             >  1:10360075/1‑150 (MQ=255)
                                                                                            tGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTc                                                           <  2:4584812/150‑1 (MQ=255)
                                                                                            tGAGGCCACCGTGATCGGCCTGTCCCGGCGTCTCAGCCCCGGGCGCGCGTGAACCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTc                                                           <  1:8007639/150‑1 (MQ=255)
                                                                                              ctcccACCGTGCTCGGCCTTTTCCGGCGTCACACCCCCCGGCGCGCGGCAACCCGGGCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCt                                                         <  1:979916/147‑1 (MQ=255)
                                                                                              aggCCACGGTGATCGGCCTTTCCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCt                                                         <  1:1523070/148‑1 (MQ=255)
                                                                                              aggCCACGGTGATCGGCCTTTCCCGGCGACTCACCACCCGGCGCGCGGCAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCt                                                         <  1:2083950/148‑1 (MQ=255)
                                                                                              aggCCACCGTGATCGGCCTTTGCCGGCTTCACATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCt                                                         <  2:10702588/148‑1 (MQ=255)
                                                                                              aTGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCt                                                         >  1:10989863/1‑150 (MQ=255)
                                                                                                 ccACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACGTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGa                                                      <  2:7044297/150‑1 (MQ=255)
                                                                                                 ccACCGTGAGCGGCCTTGCCCGGCGTCGCATCCCCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGa                                                      <  1:2566808/150‑1 (MQ=255)
                                                                                                  cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAc                                                     <  2:3680689/150‑1 (MQ=255)
                                                                                                  cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAc                                                     >  2:8694069/1‑150 (MQ=255)
                                                                                                   aCCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAcg                                                    >  2:1748206/1‑150 (MQ=255)
                                                                                                      gTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAcgccc                                                 <  2:9270780/150‑1 (MQ=255)
                                                                                                       tGATCGGCCTTTGCCGGGGTCTCACCCCCGGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAcgcccg                                                <  2:2300866/150‑1 (MQ=255)
                                                                                                        gATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAcgcccgc                                               >  2:8057159/1‑150 (MQ=255)
                                                                                                            ggCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGccccc                                           >  1:7597157/1‑150 (MQ=255)
                                                                                                               cTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCttt                                        >  2:9450656/1‑150 (MQ=255)
                                                                                                                ttttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTc                                       >  2:8300637/1‑150 (MQ=255)
                                                                                                                ttttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTc                                       >  2:9959096/1‑150 (MQ=255)
                                                                                                                ttttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTc                                       >  2:4079148/1‑150 (MQ=255)
                                                                                                                ttttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGCCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTc                                       >  2:335294/1‑150 (MQ=255)
                                                                                                                 tttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTcc                                      <  2:7196359/150‑1 (MQ=255)
                                                                                                                     cGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGt                                  >  2:2849053/1‑150 (MQ=255)
                                                                                                                         gTCTCATCACCGGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTAcgc                              <  2:8770353/150‑1 (MQ=255)
                                                                                                                          tctcATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTAcgcg                             >  2:6828604/1‑150 (MQ=255)
                                                                                                                           ctcATCACCCGGCGCGCGGTAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCgg                            <  2:8998356/150‑1 (MQ=255)
                                                                                                                            tcatcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCggg                           >  1:8825896/1‑150 (MQ=255)
                                                                                                                            tcatcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCggg                           >  2:2958441/1‑150 (MQ=255)
                                                                                                                              atcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGtt                         >  1:8550524/1‑150 (MQ=255)
                                                                                                                                caCCCGGCGCGCGGCAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGTGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTtgt                       <  2:2476268/150‑1 (MQ=255)
                                                                                                                                  cccGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                     >  1:10822100/1‑150 (MQ=255)
                                                                                                                                  cccGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                     <  1:7873383/150‑1 (MQ=255)
                                                                                                                                  cccGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                     >  2:1106323/1‑150 (MQ=255)
                                                                                                                                   ccGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTggg                    <  1:6765463/150‑1 (MQ=255)
                                                                                                                                         cgcgGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGtt              <  1:5934287/150‑1 (MQ=255)
                                                                                                                                         cgcgGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGtt              <  2:5172513/150‑1 (MQ=255)
                                                                                                                                          gcgGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGttt             >  2:5296913/1‑150 (MQ=255)
                                                                                                                                            gggAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCGCGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTAt           <  1:3052855/150‑1 (MQ=255)
                                                                                                                                              gAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCa         <  2:5600726/150‑1 (MQ=255)
                                                                                                                                              gAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCa         >  2:1302962/1‑150 (MQ=255)
                                                                                                                                                aaCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa       <  1:4229170/150‑1 (MQ=255)
                                                                                                                                                   cGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatc    >  1:93192/1‑150 (MQ=255)
                                                                                                                                                     tGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcat  >  1:3746886/1‑150 (MQ=255)
                                                                                                                                                     tGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcat  <  2:6475559/150‑1 (MQ=255)
                                                                                                                                                     tGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcat  <  1:10014153/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCAT  >  NC_000913/2718169‑2718467

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: