Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 826,945 A→G 5.9% L103P (CTG→CCG)  ybhS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913826,9450AG 5.9% 339.8 / 8.6 101L103P (CTG→CCG) ybhSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (49/46);  new base G (3/3);  total (52/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 106 total aligned reads displayed.

ATTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCA  >  NC_000913/826801‑827092
                                                                                                                                                |                                                                                                                                                   
aTTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACcg                                                                                                                                                <  1:405186/150‑1 (MQ=255)
aTTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACcg                                                                                                                                                <  1:402574/150‑1 (MQ=255)
     ccagatctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAAt                                                                                                                                           >  1:1601287/1‑150 (MQ=255)
        gatctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAAtttt                                                                                                                                        <  1:58028/150‑1 (MQ=255)
        gatctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAAtttt                                                                                                                                        <  1:58069/150‑1 (MQ=255)
        gatctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAAtttt                                                                                                                                        <  1:705967/150‑1 (MQ=255)
           ctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTccc                                                                                                                                     <  2:1451049/150‑1 (MQ=255)
            tgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTcccc                                                                                                                                    >  2:2208856/1‑150 (MQ=255)
            tgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTcccc                                                                                                                                    >  1:2376416/1‑150 (MQ=255)
                agatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCt                                                                                                                                >  1:2050389/1‑150 (MQ=255)
                     ccTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtt                                                                                                                           >  1:1545661/1‑150 (MQ=255)
                     ccTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtt                                                                                                                           >  1:1886114/1‑150 (MQ=255)
                     ccTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtt                                                                                                                           >  1:224974/1‑150 (MQ=255)
                     ccTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtt                                                                                                                           >  1:970237/1‑150 (MQ=255)
                       ttCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtttt                                                                                                                         <  2:560469/150‑1 (MQ=255)
                       ttCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAtttt                                                                                                                         <  2:261764/150‑1 (MQ=255)
                             aTACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTTCCCGCCTGCATTTTGGCGAt                                                                                                                   <  2:2187884/150‑1 (MQ=255)
                               aCCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCa                                                                                                                 <  2:331675/150‑1 (MQ=255)
                                 cccTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCGGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGt                                                                                                               >  1:234019/1‑150 (MQ=255)
                                 cccTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCGGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGt                                                                                                               >  1:234025/1‑150 (MQ=255)
                                 cccTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCGGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGt                                                                                                               >  1:234054/1‑150 (MQ=255)
                                 cccTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCa                                                                                                                 >  2:1662673/1‑148 (MQ=255)
                                   cTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTc                                                                                                             >  2:1328852/1‑150 (MQ=255)
                                   cTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTc                                                                                                             >  2:1213343/1‑150 (MQ=255)
                                   cTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTc                                                                                                             >  2:2224427/1‑150 (MQ=255)
                                   cTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTc                                                                                                             >  2:785446/1‑150 (MQ=255)
                                    tGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTcc                                                                                                            <  1:575914/150‑1 (MQ=255)
                                    tGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTcc                                                                                                            <  1:123938/150‑1 (MQ=255)
                                        cAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGAc                                                                                                        <  2:1265694/150‑1 (MQ=255)
                                          aaGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACgg                                                                                                      >  2:300565/1‑150 (MQ=255)
                                           aGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGt                                                                                                     <  1:2231438/150‑1 (MQ=255)
                                             ttAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTa                                                                                                   <  2:868797/150‑1 (MQ=255)
                                                gCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCg                                                                                                <  1:2208856/150‑1 (MQ=255)
                                                 cGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGc                                                                                               <  1:458430/150‑1 (MQ=255)
                                                 cGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGc                                                                                               <  1:1613646/150‑1 (MQ=255)
                                                 cGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGc                                                                                               <  1:1243817/150‑1 (MQ=255)
                                                  ggTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCt                                                                                              >  2:777933/1‑150 (MQ=255)
                                                   gTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTg                                                                                             >  1:2133916/1‑150 (MQ=255)
                                                    tATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGa                                                                                            <  1:1384344/150‑1 (MQ=255)
                                                     aTTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAt                                                                                           >  1:618662/1‑150 (MQ=255)
                                                         ggCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtg                                                                                       >  1:108102/1‑150 (MQ=255)
                                                          gCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtgg                                                                                      >  2:766965/1‑150 (MQ=255)
                                                             cACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGt                                                                                   <  1:744260/150‑1 (MQ=255)
                                                             cACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGt                                                                                   <  1:2207216/150‑1 (MQ=255)
                                                                tGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGa                                                                                >  2:634418/1‑150 (MQ=255)
                                                                    gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGTCGGTTATCGCTGATGGTGGCGTCGTTGTa                                                                            >  2:1355232/1‑150 (MQ=255)
                                                                    gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTa                                                                            >  2:255843/1‑150 (MQ=255)
                                                                    gTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTa                                                                            >  2:2324010/1‑150 (MQ=255)
                                                                      cGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAgg                                                                          <  1:1088021/150‑1 (MQ=255)
                                                                      cGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAgg                                                                          <  1:2274551/150‑1 (MQ=255)
                                                                       ggTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAggg                                                                         >  2:908417/1‑150 (MQ=255)
                                                                       ggTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAggg                                                                         >  2:908402/1‑150 (MQ=255)
                                                                         tGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCg                                                                       >  2:1542381/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:1540894/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:2149634/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:2212264/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:1623467/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:112975/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:644523/1‑150 (MQ=255)
                                                                            tCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCAGTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAc                                                                    >  1:1828854/1‑150 (MQ=255)
                                                                                cTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGt                                                                >  1:1304355/1‑150 (MQ=255)
                                                                                            ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGTTCTGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCAGGGTGTGGGt                                                    <  2:1986368/150‑1 (MQ=255)
                                                                                            ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGt                                                    >  2:486713/1‑150 (MQ=255)
                                                                                            ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGt                                                    <  1:1126811/150‑1 (MQ=255)
                                                                                            ggtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGt                                                    <  2:915135/150‑1 (MQ=255)
                                                                                             gtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTATGGCGAACCGGTCATGGTGTGGGTg                                                   >  1:1997075/1‑150 (MQ=255)
                                                                                             gtgGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTg                                                   >  1:252830/1‑150 (MQ=255)
                                                                                                gCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGaaa                                                >  2:515480/1‑150 (MQ=255)
                                                                                                       cgcgcTACATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAgcg                                         <  1:239595/150‑1 (MQ=255)
                                                                                                        gcgcTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAgcgc                                        >  1:1748709/1‑150 (MQ=255)
                                                                                                            tCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCt                                    >  2:1241342/1‑150 (MQ=255)
                                                                                                            tCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCt                                    >  1:1232525/1‑150 (MQ=255)
                                                                                                             ccATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCtt                                   >  1:1426179/1‑150 (MQ=255)
                                                                                                              cATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTc                                  <  1:520573/150‑1 (MQ=255)
                                                                                                                tCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGc                                <  1:349938/150‑1 (MQ=255)
                                                                                                                tCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGc                                <  1:2298756/150‑1 (MQ=255)
                                                                                                                tCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGc                                <  1:2374928/150‑1 (MQ=255)
                                                                                                                tCTGTTCCGAAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGc                                <  1:963259/150‑1 (MQ=255)
                                                                                                                    ttCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACg                            <  2:1648260/150‑1 (MQ=255)
                                                                                                                    ttCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACg                            <  2:644259/150‑1 (MQ=255)
                                                                                                                    ttCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACg                            <  2:1121440/150‑1 (MQ=255)
                                                                                                                    ttCCACAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACg                            <  2:2311831/150‑1 (MQ=255)
                                                                                                                     tCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                           <  1:2163216/150‑1 (MQ=255)
                                                                                                                     tCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                           <  2:2376416/150‑1 (MQ=255)
                                                                                                                     tCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                           <  1:2109324/150‑1 (MQ=255)
                                                                                                                     tCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                           <  1:884775/150‑1 (MQ=255)
                                                                                                                     tCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGc                           <  1:884783/150‑1 (MQ=255)
                                                                                                                       cGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTg                         >  1:1794343/1‑150 (MQ=255)
                                                                                                                             aTCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAg                   <  1:2132473/150‑1 (MQ=255)
                                                                                                                              tCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGt                  <  1:490647/150‑1 (MQ=255)
                                                                                                                                        tAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGa        >  1:165735/1‑150 (MQ=255)
                                                                                                                                        tAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGa        >  1:165762/1‑150 (MQ=255)
                                                                                                                                        tAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAGGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGa        >  2:1861635/1‑150 (MQ=255)
                                                                                                                                          accaccGGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAcc      <  2:234054/150‑1 (MQ=255)
                                                                                                                                          accaccGGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAcc      <  2:234025/150‑1 (MQ=255)
                                                                                                                                          accaccGGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAcc      <  2:234019/150‑1 (MQ=255)
                                                                                                                                             accaGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGc   >  1:2155415/1‑150 (MQ=255)
                                                                                                                                              ccaGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa  <  2:2383984/150‑1 (MQ=255)
                                                                                                                                              ccaGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa  <  2:2383982/150‑1 (MQ=255)
                                                                                                                                              ccaGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCa  >  1:263983/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCCGCA  >  NC_000913/826801‑827092

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