Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,859 Δ1 bp 10.2% coding (679/1098 nt) mltA ← membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,939,8590A.10.2% 604.7 / 80.8 206coding (679/1098 nt)mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base A (101/84);  new base . (11/10);  total (112/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

CGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATA  >  CP009273/2939712‑2940004
                                                                                                                                                   |                                                                                                                                                 
cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGa                                                                                                                                                   >  2:3157575/1‑148 (MQ=255)
cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGa                                                                                                                                                   <  1:417142/148‑1 (MQ=255)
    accaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                                               >  1:3015481/1‑148 (MQ=255)
    accaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                                               >  2:715607/1‑148 (MQ=255)
     ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACatat                                                                                                                                              >  1:1350855/1‑148 (MQ=255)
     ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACatat                                                                                                                                              >  2:360568/1‑148 (MQ=255)
     ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACatat                                                                                                                                              >  1:3517129/1‑148 (MQ=255)
        accaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtctt                                                                                                                                           >  1:2419034/1‑148 (MQ=255)
         ccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttc                                                                                                                                          <  2:186848/148‑1 (MQ=255)
         ccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttc                                                                                                                                          >  2:892760/1‑148 (MQ=255)
           aGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATcttctt                                                                                                                                        <  2:2421118/148‑1 (MQ=255)
             cGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttt                                                                                                                                      >  2:2075773/1‑148 (MQ=255)
             cGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCttttt                                                                                                                                     >  1:2506520/1‑148 (MQ=255)
               gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttttt                                                                                                                                    <  1:48745/148‑1 (MQ=255)
               gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttttt                                                                                                                                    >  2:2243391/1‑148 (MQ=255)
                ctccGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCa                                                                                                                                  <  1:2059722/146‑1 (MQ=255)
                  ccGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACt                                                                                                                                >  2:1089880/1‑148 (MQ=255)
                   cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACtt                                                                                                                               <  2:1715206/148‑1 (MQ=255)
                    gCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACtt                                                                                                                               <  2:3039016/148‑1 (MQ=255)
                    gCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACtt                                                                                                                               >  2:1603738/1‑148 (MQ=255)
                      aCTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCg                                                                                                                             <  2:3246097/148‑1 (MQ=255)
                      aCTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCg                                                                                                                             <  2:1471339/148‑1 (MQ=255)
                         tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGATAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCTGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                                          <  2:273903/148‑1 (MQ=255)
                         tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                                          <  2:2364682/148‑1 (MQ=255)
                             ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGt                                                                                                                      <  1:1161597/148‑1 (MQ=255)
                                ttCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGAt                                                                                                                  <  2:2336007/148‑1 (MQ=255)
                                 tCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGAt                                                                                                                  <  2:2938013/148‑1 (MQ=255)
                                  cACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATc                                                                                                                 >  1:1742162/1‑148 (MQ=255)
                                         gCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCt                                                                                                          <  1:3116825/148‑1 (MQ=255)
                                          cAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCtt                                                                                                         >  1:660948/1‑148 (MQ=255)
                                           aaaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTa                                                                                                        <  1:1528809/148‑1 (MQ=255)
                                            aaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAc                                                                                                       <  2:2291320/148‑1 (MQ=255)
                                             aaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAcc                                                                                                      >  2:3120126/1‑148 (MQ=255)
                                             aaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAcc                                                                                                      >  1:2925811/1‑148 (MQ=255)
                                             aaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAcc                                                                                                      >  1:1148334/1‑148 (MQ=255)
                                              aGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCa                                                                                                     >  1:1907168/1‑148 (MQ=255)
                                                aTTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                                                   <  2:1295615/148‑1 (MQ=255)
                                                aTTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtg                                                                                                  <  2:2506520/148‑1 (MQ=255)
                                                 ttGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtg                                                                                                  <  1:2081672/148‑1 (MQ=255)
                                                 ttGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtg                                                                                                  <  1:1999339/148‑1 (MQ=255)
                                                  tGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgc                                                                                                 <  2:512428/148‑1 (MQ=255)
                                                       tttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGa                                                                                            >  2:1533727/1‑148 (MQ=255)
                                                        ttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGAt                                                                                           >  1:2580703/1‑148 (MQ=255)
                                                            agaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                                                       <  1:2850308/148‑1 (MQ=255)
                                                            agaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                                                       <  2:1599125/148‑1 (MQ=255)
                                                             gaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCAt                                                                                     >  1:343674/1‑148 (MQ=255)
                                                              aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCAt                                                                                     >  1:3205352/1‑148 (MQ=255)
                                                                gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                                                   >  2:821281/1‑148 (MQ=255)
                                                                gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                                                   >  2:822512/1‑148 (MQ=255)
                                                                 aCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAc                                                                                  >  1:1497325/1‑148 (MQ=255)
                                                                     aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtt                                                                              >  1:1816837/1‑148 (MQ=255)
                                                                     aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtt                                                                              >  2:3322104/1‑148 (MQ=255)
                                                                     aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtt                                                                              >  1:12685/1‑148 (MQ=255)
                                                                         cGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTc                                                                          <  1:821045/148‑1 (MQ=255)
                                                                           ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTccc                                                                        <  2:635035/148‑1 (MQ=255)
                                                                             ttCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTg                                                                      >  2:1014643/1‑148 (MQ=255)
                                                                                tGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCAt                                                                   >  2:1674075/1‑148 (MQ=255)
                                                                                 gTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATa                                                                  >  1:606206/1‑148 (MQ=255)
                                                                                       gcagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaa                                                            >  2:155669/1‑148 (MQ=255)
                                                                                       gcagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaa                                                            >  1:1804100/1‑148 (MQ=255)
                                                                                          gcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaaaaa                                                         >  2:757725/1‑148 (MQ=255)
                                                                                             cgcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGtt                                                      <  2:1557786/148‑1 (MQ=255)
                                                                                               cgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTaa                                                    <  2:3485301/148‑1 (MQ=255)
                                                                                                   ccTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCgg                                                <  2:1583307/148‑1 (MQ=255)
                                                                                                       ggCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTg                                            >  1:2403919/1‑148 (MQ=255)
                                                                                                       ggCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGc                                           >  2:1273079/1‑148 (MQ=255)
                                                                                                          ttCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGGTAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCa                                         >  2:2985458/1‑148 (MQ=255)
                                                                                                           tCACTGTGTGTTTCGCCACAGTGACGAATCGCCTGCATCGACATATATTATTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCAt                                        <  1:1574599/148‑1 (MQ=255)
                                                                                                            cACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCAcc                                                                                                           <  2:668313/80‑1 (MQ=255)
                                                                                                                gtgtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCa                                   >  2:3560710/1‑148 (MQ=255)
                                                                                                                 tgtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCaa                                  >  1:381887/1‑148 (MQ=255)
                                                                                                                  gtgtTTCGCCACAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCATTGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCaaa                                 <  1:2211432/148‑1 (MQ=255)
                                                                                                                    gtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGt                               <  1:974991/148‑1 (MQ=255)
                                                                                                                     tttCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTc                              <  1:715607/148‑1 (MQ=255)
                                                                                                                     tttCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTc                              >  2:1250350/1‑148 (MQ=255)
                                                                                                                          ccccAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGatat                         >  2:3494648/1‑148 (MQ=255)
                                                                                                                          ccccAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGAtata                        <  1:1089880/148‑1 (MQ=255)
                                                                                                                            ccAGTGACTAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAc                       <  2:1399547/148‑1 (MQ=255)
                                                                                                                            ccAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAcc                      >  1:683093/1‑148 (MQ=255)
                                                                                                                             cAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAcc                      <  1:2344283/148‑1 (MQ=255)
                                                                                                                              aGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAccc                     >  2:2007089/1‑148 (MQ=255)
                                                                                                                               gTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCa                    <  2:2808462/148‑1 (MQ=255)
                                                                                                                               gTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCa                    <  1:553290/148‑1 (MQ=255)
                                                                                                                                tGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCAc                   >  1:1627745/1‑148 (MQ=255)
                                                                                                                                 gACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGTTAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACt                  <  1:2985458/148‑1 (MQ=255)
                                                                                                                                  aCGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTa                 >  2:406542/1‑148 (MQ=255)
                                                                                                                                   cGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTAc                <  1:2243391/148‑1 (MQ=255)
                                                                                                                                    gAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTAcc               <  1:1964743/148‑1 (MQ=255)
                                                                                                                                     aaTCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTAcc               <  1:44505/147‑1 (MQ=255)
                                                                                                                                     aaTCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTAcc               >  2:3436054/1‑147 (MQ=255)
                                                                                                                                     aaTCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTAcc               >  2:44505/1‑147 (MQ=255)
                                                                                                                                      aTCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAACCACTACCCt             >  1:1145322/1‑148 (MQ=255)
                                                                                                                                      aTCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTg            >  2:1548504/1‑148 (MQ=255)
                                                                                                                                       tCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTg            >  1:1374647/1‑148 (MQ=255)
                                                                                                                                       tCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGc           <  1:2047708/148‑1 (MQ=255)
                                                                                                                                          ccTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGcac         <  2:3184499/148‑1 (MQ=255)
                                                                                                                                            tGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACAt       >  1:3401180/1‑148 (MQ=255)
                                                                                                                                             gCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATc      >  2:218397/1‑148 (MQ=255)
                                                                                                                                             gCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATcc     >  2:973618/1‑148 (MQ=255)
                                                                                                                                                 cGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATa  >  1:793032/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATA  >  CP009273/2939712‑2940004

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: