Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,832 C→A 5.0% 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,8320CA 5.0% 953.3 / 17.4 257E236* (GAA→TAA) mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base C (91/153);  new base A (6/7);  total (97/160)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 265 total aligned reads displayed.

AATAATGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATAT  >  CP009273/2939683‑2939981
                                                                                                                                                     |                                                                                                                                                     
aataatGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTc                                                                                                                                                       <  2:2617012/150‑1 (MQ=255)
aataatGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTc                                                                                                                                                       >  2:1701232/1‑150 (MQ=255)
aataatGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCAAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTc                                                                                                                                                       >  2:1774302/1‑150 (MQ=255)
 ataatGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCg                                                                                                                                                      <  1:948360/150‑1 (MQ=255)
    atGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGccc                                                                                                                                                   <  2:1509864/150‑1 (MQ=255)
      ggAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCa                                                                                                                                                 <  2:1445222/150‑1 (MQ=255)
      ggAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCa                                                                                                                                                 <  2:2347393/150‑1 (MQ=255)
       gAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAg                                                                                                                                                <  1:2030526/150‑1 (MQ=255)
        aaCGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGt                                                                                                                                               <  2:4325617/150‑1 (MQ=255)
          cGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGa                                                                                                                                             <  1:3236392/150‑1 (MQ=255)
          cGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGa                                                                                                                                             <  1:3115002/150‑1 (MQ=255)
           gATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGAc                                                                                                                                            >  1:3524063/1‑150 (MQ=255)
            aTCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACg                                                                                                                                           >  1:4489176/1‑150 (MQ=255)
            aTCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACg                                                                                                                                           >  2:3697383/1‑150 (MQ=255)
             tCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGa                                                                                                                                          <  1:1798349/150‑1 (MQ=255)
                gagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATc                                                                                                                                       >  1:2287807/1‑150 (MQ=255)
                 agGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCg                                                                                                                                      <  1:1218208/150‑1 (MQ=255)
                  ggCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGc                                                                                                                                     >  2:4593819/1‑150 (MQ=255)
                   gCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGcc                                                                                                                                    >  1:166140/1‑150 (MQ=255)
                     aaCTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTg                                                                                                                                  <  2:2284696/150‑1 (MQ=255)
                      aCTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGc                                                                                                                                 <  2:887437/150‑1 (MQ=255)
                      aCTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGc                                                                                                                                 <  1:3009172/150‑1 (MQ=255)
                       cTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCa                                                                                                                                <  2:4952677/150‑1 (MQ=255)
                         gACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATc                                                                                                                              <  2:2885431/150‑1 (MQ=255)
                          aCGCGCGACCAACCAGCGTCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGTAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTTTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCg                                                                                                                             <  1:2589774/150‑1 (MQ=255)
                          aCGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATca                                                                                                                             <  1:1579779/150‑2 (MQ=255)
                           cgcgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGa                                                                                                                            <  2:316396/150‑1 (MQ=255)
                            gcgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGAc                                                                                                                           <  2:2857211/150‑1 (MQ=255)
                             cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                          <  1:2774324/150‑1 (MQ=255)
                             cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                          >  1:1741493/1‑150 (MQ=255)
                             cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                          <  1:4542456/150‑1 (MQ=255)
                             cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                          <  2:3445994/150‑1 (MQ=255)
                                  ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtct                                                                                                                     >  2:1330975/1‑150 (MQ=255)
                                    aaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATCtcttc                                                                                                                   >  1:559663/1‑150 (MQ=255)
                                     accaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttct                                                                                                                  <  1:2423785/150‑1 (MQ=255)
                                     accaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttct                                                                                                                  >  1:1966488/1‑150 (MQ=255)
                                     accaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtcttct                                                                                                                  >  1:4301695/1‑150 (MQ=255)
                                       caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCttt                                                                                                                <  1:611614/150‑1 (MQ=255)
                                       caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCtttt                                                                                                               >  2:3264905/1‑150 (MQ=255)
                                         gCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCttttt                                                                                                              <  2:2181142/150‑1 (MQ=255)
                                           ggCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                            >  2:124021/1‑150 (MQ=255)
                                               ccGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACtt                                                                                                        <  1:2369565/150‑1 (MQ=255)
                                                cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTc                                                                                                       >  2:4717244/1‑150 (MQ=255)
                                                cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTc                                                                                                       <  2:417315/150‑1 (MQ=255)
                                                 gCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCg                                                                                                      <  1:1902847/150‑1 (MQ=255)
                                                  cACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGc                                                                                                     >  2:4109940/1‑150 (MQ=255)
                                                      tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTGCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACg                                                                                                 <  1:2272776/150‑1 (MQ=255)
                                                        ccccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGt                                                                                               <  2:2321438/150‑1 (MQ=255)
                                                            tttCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGAt                                                                                           >  1:470585/1‑150 (MQ=255)
                                                                     aGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCtt                                                                                  >  2:3727919/1‑150 (MQ=255)
                                                                          aaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCaa                                                                             <  2:2913318/150‑1 (MQ=255)
                                                                           aGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                            <  1:1160507/150‑1 (MQ=255)
                                                                           aGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                            <  1:2865785/150‑1 (MQ=255)
                                                                           aGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                            <  1:4086921/150‑1 (MQ=255)
                                                                            gATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtg                                                                           <  2:1431166/150‑1 (MQ=255)
                                                                             aTTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgc                                                                          >  2:2485846/1‑150 (MQ=255)
                                                                               tGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgctg                                                                        <  1:4798813/150‑1 (MQ=255)
                                                                                gCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgctgc                                                                       <  2:842485/150‑1 (MQ=255)
                                                                                gCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgctgc                                                                       <  2:2047814/150‑1 (MQ=255)
                                                                                  ggTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGa                                                                     >  2:1990422/1‑150 (MQ=255)
                                                                                  ggTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGa                                                                     <  2:4701405/150‑1 (MQ=255)
                                                                                     ttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAg                                                                  <  1:3135874/150‑1 (MQ=255)
                                                                                     ttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAg                                                                  <  1:3509782/150‑1 (MQ=255)
                                                                                       aaaGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                                >  1:1012622/1‑150 (MQ=255)
                                                                                       aaaGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                                <  2:3009996/150‑1 (MQ=255)
                                                                                        aagaagACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCAt                                                              <  1:3879763/150‑1 (MQ=255)
                                                                                          gaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATg                                                             <  2:4877952/150‑1 (MQ=255)
                                                                                          gaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATg                                                             <  2:2140900/150‑1 (MQ=255)
                                                                                           aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                            >  1:1558201/1‑150 (MQ=255)
                                                                                             gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAcc                                                          <  1:759628/150‑1 (MQ=255)
                                                                                               cGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGt                                                        <  1:3756128/150‑1 (MQ=255)
                                                                                               cGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGt                                                        >  2:811798/1‑150 (MQ=255)
                                                                                                gAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtt                                                       >  2:3372871/1‑150 (MQ=255)
                                                                                                  aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtttt                                                     <  1:4646022/150‑1 (MQ=255)
                                                                                                  aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCATACCAATGCTGCGATAGGCATGACCGtttt                                                     >  2:4759693/1‑150 (MQ=255)
                                                                                                  aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtttt                                                     >  1:3339687/1‑150 (MQ=255)
                                                                                                       gggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                                <  2:919257/150‑1 (MQ=255)
                                                                                                       gggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                                >  2:4002091/1‑150 (MQ=255)
                                                                                                        ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTg                                               <  2:580995/150‑1 (MQ=255)
                                                                                                          ttCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCa                                             >  1:2806060/1‑150 (MQ=255)
                                                                                                           tCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCAt                                            <  2:1084611/150‑1 (MQ=255)
                                                                                                             tGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAg                                          >  1:209918/1‑150 (MQ=255)
                                                                                                               ttCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCt                                        <  2:1144177/150‑1 (MQ=255)
                                                                                                                tCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTg                                       <  2:4460052/150‑1 (MQ=255)
                                                                                                                tCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTg                                       >  2:782660/1‑150 (MQ=255)
                                                                                                                   agcagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaaa                                    >  1:127242/1‑150 (MQ=255)
                                                                                                                           gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTaan                            <  1:3389628/150‑2 (MQ=255)
                                                                                                                           gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAg                            <  2:559663/150‑1 (MQ=255)
                                                                                                                            cgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGc                           >  1:2311339/1‑150 (MQ=255)
                                                                                                                                 cTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTg                     >  1:2903783/1‑150 (MQ=255)
                                                                                                                                     gCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCAt                 >  2:2365286/1‑150 (MQ=255)
                                                                                                                                          aCTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCAcc             <  1:706265/150‑1 (MQ=255)
                                                                                                                                          aCTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCAcc             >  1:3649306/1‑150 (MQ=255)
                                                                                                                                             gtgtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCaaa          <  1:2973441/150‑1 (MQ=255)
                                                                                                                                              tgtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAg         >  2:5111960/1‑150 (MQ=255)
                                                                                                                                               gtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGt        <  1:3576725/150‑1 (MQ=255)
                                                                                                                                               gtgtTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTc       >  2:95861/1‑150 (MQ=255)
                                                                                                                                                  tttCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGa     <  1:4590392/150‑1 (MQ=255)
                                                                                                                                                  tttCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGat    <  1:1103290/150‑1 (MQ=255)
                                                                                                                                                     cGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGatat  <  2:2084141/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATAATGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATAT  >  CP009273/2939683‑2939981

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: