Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,832 C→A 10.9% E236* (GAA→TAA)  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,8320CA10.9% 712.5 / 40.6 210E236* (GAA→TAA) mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base C (97/90);  new base A (13/10);  total (110/100)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 211 total aligned reads displayed.

AATGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGAT  >  CP009273/2939686‑2939979
                                                                                                                                                  |                                                                                                                                                   
aatGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCg                                                                                                                                                    >  1:2078909/1‑148 (MQ=255)
 atGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGc                                                                                                                                                   >  1:61664/1‑148 (MQ=255)
  tGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGcc                                                                                                                                                  >  1:1578881/1‑148 (MQ=255)
  tGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGcc                                                                                                                                                  >  2:2730941/1‑148 (MQ=255)
   ggAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGccc                                                                                                                                                 <  2:1201529/148‑1 (MQ=255)
   ggAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGccc                                                                                                                                                 <  1:1395339/148‑1 (MQ=255)
      aCGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAg                                                                                                                                              <  2:805944/148‑1 (MQ=255)
        gATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTg                                                                                                                                            <  2:2293823/148‑1 (MQ=255)
        gATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTACAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTg                                                                                                                                            <  1:288944/148‑1 (MQ=255)
         aTCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGa                                                                                                                                           <  2:271812/148‑1 (MQ=255)
         aTCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGa                                                                                                                                           >  1:1023245/1‑148 (MQ=255)
           cAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACg                                                                                                                                         >  1:1426530/1‑148 (MQ=255)
            agagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGa                                                                                                                                        >  2:3101960/1‑148 (MQ=255)
             gagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGaa                                                                                                                                       <  1:659686/148‑1 (MQ=255)
             gagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGaa                                                                                                                                       <  2:1770538/148‑1 (MQ=255)
                gCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCg                                                                                                                                    >  1:2345025/1‑148 (MQ=255)
                   aCTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCt                                                                                                                                 >  1:2946241/1‑148 (MQ=255)
                    cTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTg                                                                                                                                >  1:2546357/1‑148 (MQ=255)
                     tGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGc                                                                                                                               >  1:323562/1‑148 (MQ=255)
                            cgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                        <  1:1822957/148‑1 (MQ=255)
                              accaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                      >  1:3473866/1‑148 (MQ=255)
                               ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACatat                                                                                                                     >  1:727852/1‑148 (MQ=255)
                               ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACatat                                                                                                                     <  2:2009398/148‑1 (MQ=255)
                                  accaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtctt                                                                                                                  >  1:1740717/1‑148 (MQ=255)
                                  accaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtctt                                                                                                                  >  1:34353/1‑148 (MQ=255)
                                    caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttct                                                                                                                >  2:2224328/1‑148 (MQ=255)
                                    caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttct                                                                                                                <  1:2366950/148‑1 (MQ=255)
                                     aGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATcttctt                                                                                                               >  1:2494590/1‑148 (MQ=255)
                                     aGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCttt                                                                                                              <  2:2993136/148‑1 (MQ=255)
                                       cGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttt                                                                                                             <  1:2218570/148‑1 (MQ=255)
                                       cGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttt                                                                                                             <  1:2565337/148‑1 (MQ=255)
                                        ggCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCttttt                                                                                                            >  2:3378478/1‑148 (MQ=255)
                                         gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttttt                                                                                                           >  2:1516445/1‑148 (MQ=255)
                                         gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttttt                                                                                                           <  2:944433/148‑1 (MQ=255)
                                         gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttttt                                                                                                           <  1:1037745/148‑1 (MQ=255)
                                          cACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                          <  2:3136532/148‑1 (MQ=255)
                                          cACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                          <  2:1366997/148‑1 (MQ=255)
                                          cACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCa                                                                                                         >  2:2879359/1‑148 (MQ=255)
                                             cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACt                                                                                                       >  1:3397774/1‑148 (MQ=255)
                                                aCTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCg                                                                                                    <  2:2678757/148‑1 (MQ=255)
                                                  ttGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATATTCTTTTTTCACTTCGcc                                                                                                  <  1:1313117/148‑1 (MQ=255)
                                                   tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                 >  1:526788/1‑148 (MQ=255)
                                                   tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                 <  1:382243/148‑1 (MQ=255)
                                                   tGCCCCTTTCACCGAAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                 >  1:2080412/1‑148 (MQ=255)
                                                         tttCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCg                                                                                           <  1:2417052/148‑1 (MQ=255)
                                                           tCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGAt                                                                                         <  1:519286/148‑1 (MQ=255)
                                                              ccGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAg                                                                                      >  1:3536052/1‑148 (MQ=255)
                                                              ccGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAg                                                                                      <  2:395477/148‑1 (MQ=255)
                                                                    cAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTa                                                                               >  2:3013364/1‑148 (MQ=255)
                                                                      aaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAc                                                                              >  1:570402/1‑148 (MQ=255)
                                                                      aaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAc                                                                              >  1:978033/1‑148 (MQ=255)
                                                                       aaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAcc                                                                             >  2:2171747/1‑148 (MQ=255)
                                                                         gATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCaa                                                                           >  1:2652339/1‑148 (MQ=255)
                                                                           ttGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtg                                                                         <  2:2282426/148‑1 (MQ=255)
                                                                           ttGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtg                                                                         <  1:1847446/148‑1 (MQ=255)
                                                                             gCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgctg                                                                      <  2:1683799/148‑1 (MQ=255)
                                                                                 tttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGa                                                                   <  2:1563192/148‑1 (MQ=255)
                                                                                  ttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGAt                                                                  <  1:3384798/148‑1 (MQ=255)
                                                                                  ttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGAt                                                                  <  1:201405/148‑1 (MQ=255)
                                                                                   tAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATa                                                                 >  1:3309466/1‑148 (MQ=255)
                                                                                   tAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATa                                                                 >  2:2167170/1‑148 (MQ=255)
                                                                                        aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCAt                                                            >  1:1774226/1‑148 (MQ=255)
                                                                                          gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                          >  1:119080/1‑148 (MQ=255)
                                                                                          gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                          <  2:2371632/148‑1 (MQ=255)
                                                                                          gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAc                                                         >  1:1061154/1‑148 (MQ=255)
                                                                                           aCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAc                                                         >  1:781564/1‑148 (MQ=255)
                                                                                           aCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAc                                                         <  1:2151482/148‑1 (MQ=255)
                                                                                             gAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCg                                                       <  1:196827/148‑1 (MQ=255)
                                                                                                 gACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtttt                                                   >  2:2941008/1‑148 (MQ=255)
                                                                                                  aCGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGttttt                                                  >  2:3286012/1‑148 (MQ=255)
                                                                                                     ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTccc                                               >  1:399926/1‑148 (MQ=255)
                                                                                                      gTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                              >  1:3369477/1‑148 (MQ=255)
                                                                                                       ttCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTg                                             <  2:673991/148‑1 (MQ=255)
                                                                                                           gTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATa                                         >  1:1717116/1‑148 (MQ=255)
                                                                                                           gTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATa                                         <  2:1161123/148‑1 (MQ=255)
                                                                                                           gTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATa                                         >  2:43055/1‑148 (MQ=255)
                                                                                                             tCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGc                                       <  2:1049949/148‑1 (MQ=255)
                                                                                                                  cagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaaaa                                 >  2:2706322/1‑148 (MQ=255)
                                                                                                                     cTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAg                               <  2:3397807/148‑1 (MQ=255)
                                                                                                                      tCGCGACCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGt                              <  2:1027650/148‑1 (MQ=255)
                                                                                                                        gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTa                            <  2:1884772/148‑1 (MQ=255)
                                                                                                                        gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTa                            <  2:1882546/148‑1 (MQ=255)
                                                                                                                        gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTa                            <  1:1868444/148‑1 (MQ=255)
                                                                                                                         cgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTaa                           <  1:1836013/148‑1 (MQ=255)
                                                                                                                            aCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCg                        <  2:1426450/148‑1 (MQ=255)
                                                                                                                            aCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCg                        <  2:3193548/148‑1 (MQ=255)
                                                                                                                            aCCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCg                        >  1:694039/1‑148 (MQ=255)
                                                                                                                              cTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGa                      <  2:3149080/148‑1 (MQ=255)
                                                                                                                                   cTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGcc                 >  1:2302963/1‑148 (MQ=255)
                                                                                                                                    ttCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAATCGGACTGCCa                >  2:2553043/1‑148 (MQ=255)
                                                                                                                                    ttCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCAt               <  2:323562/148‑1 (MQ=255)
                                                                                                                                     tCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCAt               >  1:3050863/1‑148 (MQ=255)
                                                                                                                                          gtgtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCAACGGTCGATCAGAACCTTACCAATACTACCATATGAATGACCGTTTTTAACTGCAAAGCTTAAAAATTTAAGCCGACTGCCATCACCa          >  1:1519314/1‑148 (MQ=255)
                                                                                                                                           tgtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCaa         <  1:3175074/148‑1 (MQ=255)
                                                                                                                                             tgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAg       >  1:3060821/1‑148 (MQ=255)
                                                                                                                                              gtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGt      >  2:2755946/1‑148 (MQ=255)
                                                                                                                                                ttCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGTCTGCCATCACCAAAGTCg    >  1:3471833/1‑148 (MQ=255)
                                                                                                                                                ttCGCCCCAGTGACGAATCGCCTGCATCG‑CATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGa   <  1:1275267/148‑1 (MQ=255)
                                                                                                                                                 tCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGa   <  2:2080412/148‑1 (MQ=255)
                                                                                                                                                  cGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGat  <  2:345683/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGAT  >  CP009273/2939686‑2939979

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: ‑