Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,555,047 T→C 6.8% D47G (GAT→GGT)  ddpB ← D‑ala‑D‑ala transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,555,0470TC 6.8% 2263.9 / 9.6 836D47G (GAT→GGT) ddpBD‑ala‑D‑ala transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (470/245);  new base C (2/92);  total (478/358)
Fisher's exact test for biased strand distribution p-value = 4.65e-36
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 845 total aligned reads displayed.

AAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACT  >  CP009273/1554898‑1555196
                                                                                                                                                     |                                                                                                                                                     
  aaTATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCg                                                                                                                                                     <  2:407541/150‑1 (MQ=255)
   atatACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGc                                                                                                                                                    >  2:32838/1‑150 (MQ=255)
     ataCGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGGTGCTGGCGAATATTTTCCACGATAGCATCGCTg                                                                                                                                                  >  2:1188234/1‑150 (MQ=255)
        cGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgcg                                                                                                                                               <  1:15765/150‑1 (MQ=255)
        cGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgc                                                                                                                                                <  2:1028406/149‑1 (MQ=255)
        cGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgc                                                                                                                                                <  1:1140420/149‑1 (MQ=255)
         gCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATc                                                                                                                                                      <  1:563006/142‑1 (MQ=255)
         gCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgcgc                                                                                                                                              >  1:734633/1‑150 (MQ=255)
         gCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgcgc                                                                                                                                              >  2:90252/1‑150 (MQ=255)
             ctctTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTc                                                                                                                                          <  1:662871/150‑1 (MQ=255)
                  ccAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCAccc                                                                                                                                      <  2:701320/149‑1 (MQ=255)
                   cAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCg                                                                                                                                     <  2:1063076/149‑1 (MQ=255)
                       cccGGACGCCCGGTACGAATGGATGTCCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGAAAAACTGTCCATACAGTGGCTGGCCCAGCCCCAGTTGCTGGCAAATATTTCCCACGATAGCACCGCTGGCGCGGTCACCCGCCAGt                                                                                                                                <  2:128704/149‑1 (MQ=255)
                       aCCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCCGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                                                                 >  2:184842/1‑149 (MQ=255)
                       aCCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                                                                 >  2:267074/1‑149 (MQ=255)
                        ccGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCCGGCGCGGGCACCCGCCAGt                                                                                                                                >  1:315884/1‑149 (MQ=255)
                          ggACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGATGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCAACCGCCCGTAAc                                                                                                                             >  2:457795/1‑150 (MQ=255)
                             cGCCGGGTACGAAGGGATGTCCCCAGGTCACCATGAAACAGGTCGCTGACGTAGGGGTAAACCTGTACATACGGGGGCGGGCCCAGCCCCAGTTGCTGGCGAATATTTTCCAGGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACcgc                                                                                                                          <  2:813310/150‑1 (MQ=255)
                              gCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACcgcg                                                                                                                         >  1:521653/1‑150 (MQ=255)
                                ccGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGcc                                                                                                                       >  1:930783/1‑150 (MQ=255)
                                ccGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGcc                                                                                                                       >  2:599437/1‑150 (MQ=255)
                                  ggTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCg                                                                                                                      >  1:501605/1‑149 (MQ=255)
                                    tACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGAt                                                                                                                   >  1:414902/1‑150 (MQ=255)
                                      caaaTGGTTTTCCCCGGGCCACCATAAAAAAGGTCCCGGAGGTAGCGAAAAAACTGTACATACGGTGGCGGGCCCGGCCCCGGTGGCGGGCGAATTTTTCCCACGATAGCACCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCa                                                                                                                 <  2:1102790/148‑1 (MQ=255)
                                       gAATGGTTGTTCCAAGGTCACCAAGAAACAGGTCGCGGCCGTAGGGGTAAAACTGTACAGACAGTGGCTGGCCCTGTCCCAGTTGCGGGCGAATATTTCCCAGGAAAGCATCGCTGGGGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                <  1:69884/150‑1 (MQ=255)
                                       gAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGCCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                <  1:476126/150‑1 (MQ=255)
                                       gAATGGATGTCCCCAGGTCCCCATGAAACAGGTCGCTGACGTGGCGGTAAACCTGTCAATACAGTGGCTGGTCCAGCCCCGGTTGCTGGCGAATATTTCCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                <  1:998417/150‑1 (MQ=255)
                                          tGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGAACACCAg                                                                                                             >  2:113711/1‑150 (MQ=255)
                                          tGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGAt                                                                                                                   >  1:941886/1‑144 (MQ=255)
                                          tGGATGTTCCAAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAg                                                                                                             <  1:1177533/150‑1 (MQ=255)
                                           ggATGTCCCCAGGTCCCCATGAAACAGGCCGCTGACGTAGCGGTAAAACTGTACATACGGTGGCTGGCCCAGTCCCAGTTGCGGGAAAATATTTCCCACGATAGCACCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAgg                                                                                                            <  2:1085648/150‑1 (MQ=255)
                                               ttccccAGGTCACCATGAAACGGGTCGCTGACGTAGCGGTAAACCTGTACATACGGTGGCTGGTCCGGCCCCGGTTGCGGGCAAATATTTCCCCCGATAGCATCGCGGGCGGGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGAt                                                                                                         <  2:210300/146‑1 (MQ=255)
                                                ttCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGAt                                                                                                         >  1:537523/1‑148 (MQ=255)
                                                ttCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGAt                                                                                                         >  1:803079/1‑148 (MQ=255)
                                                   ccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGAt                                                                                                    >  1:1094809/1‑150 (MQ=255)
                                                   ccAGGTCACCAAGAAACAGGTCGCTGCCGTGGCGGAAAAACTGTACATACAGTGGCTGGCCCAGCCCCGGTTGCTGGCGAATATTTCCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGAt                                                                                                    <  2:168861/150‑1 (MQ=255)
                                                    cAGCTCCCCATAAAAGAGGTCGCTGCGGTGGCGTAAAAACTGTACATACGGTGCCGGATCCAGTCCCGGTTGCGGGCGAATTTTTTCCAGGATAGGATCGGTGGCGGGCTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATg                                                                                                   <  1:1121291/150‑1 (MQ=255)
                                                     aGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGc                                                                                                  >  2:716261/1‑150 (MQ=255)
                                                      ggTCACCATGAAACAGGTCGCTGACGTAGGGGTAAAACTGTACATACAGTGGCTGGTCCAGCCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCg                                                                                                 <  1:1072577/150‑1 (MQ=255)
                                                      ggTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGc                                                                                                  >  1:376832/1‑149 (MQ=255)
                                                       gTCCCCATGAAACGGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGCCCAGCCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCga                                                                                                <  1:60601/150‑1 (MQ=255)
                                                         cACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGc                                                                                                  >  1:761466/1‑146 (MQ=255)
                                                         cACCATAAAACAGGTCGCTGACGTGGCGTTAAACCTGTACATACAGTGGCTGGTCCAGTCCCATTTGCTGGCGAATTTTTCCCACGATAGCATCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCgaga                                                                                              <  1:120639/150‑1 (MQ=255)
                                                          acaaTGAAACGGGTCGCGGACGTGGGGGTAAACCTGTACATACAGTGGCTGGCCCGGTCCCAGTGGCGGGCGAATATTTCCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGat                                                                                             <  1:178390/147‑1 (MQ=255)
                                                           ccATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGat                                                                                             >  2:978097/1‑149 (MQ=255)
                                                             aTGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                                                           >  1:360321/1‑149 (MQ=255)
                                                               gAAACGGGTCGTGGCGGTAGCGGAAAAATTGTACATACGGGGGCGGGTCCAGCCCCGGTTGCTGGCAAATTTTTCCCAGGATAGCATGGCTGGCGGGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGAtaataa                                                                                        <  2:560648/150‑1 (MQ=255)
                                                                aaaCGGGTCGCTGACGGAGCGGTAAAACTGTACATACGGTGGCTGGTCCAGCCCCGGTTGCTGGCAAATATTTCCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGAtaataa                                                                                        <  1:90252/149‑1 (MQ=255)
                                                                 aaCAGGTCGCTGACGTAGGGGTAAAACTGTAAATCAAGGGGCTGGCCCAGCCCCGGTTGCGGGGAAATTTTTCCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAAc                                                                                      <  2:758524/150‑1 (MQ=255)
                                                                    gggTCGTTGACGTAGGGGAAAACCTGTCCATACGGGGGCTGGTCAAGGCCCGGTTGCTGGCAAATTTTTCCCACGATACCCCGGGGGGGGGGGCCATCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTa                                                                                   <  2:832771/149‑1 (MQ=255)
                                                                      gTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGAAAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAAAGTAAt                                                                                 >  2:245764/1‑150 (MQ=255)
                                                                         gCTGACGTGGCGGTAAACCTGTACATCAAGTGGCTGGTCCGGCCCCAGTTGCGGGCAAATATTTCCCACGATAGCACCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCAc                                                                              <  1:887274/150‑1 (MQ=255)
                                                                           tGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACg                                                                             >  1:84969/1‑149 (MQ=255)
                                                                            gccGTAGGGATAAACCTATACATACAGTGGCTGGCCCAGTCCCATTTGGTGGCAAATATTTACCGGGATAACATCGGTGGCGGGGACAACCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                                                           <  2:1045950/148‑1 (MQ=255)
                                                                            gACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCGGGCGAATATTTTCCACGATAGCACCGCTGGCGCGGCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                                                           <  2:521471/150‑1 (MQ=255)
                                                                             aCGTAGCGGAAAACCTGTACATAAAGTGGCTGGCCCAGCCCCAGTTGCGGGCAAATATTTTCCACGATAGCATCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                                                           <  1:719645/149‑1 (MQ=255)
                                                                              cGTAGCGGAAAAACTGTCCATACGGGGGCTGCCCCAGCCCCAGTTGCGGGCGAATTTTTCCCACGATAGCATCGCGGGCGCGGCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGa                                                                         <  2:681392/150‑1 (MQ=255)
                                                                               gTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAc                                                                        >  1:579848/1‑150 (MQ=255)
                                                                                tAGCGATAAAACTGTACATACGGTGGCTGGCCCAGTCCCATTTGCTGGCAAATTTTTCCCAAGATGGCATCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcg                                                                       <  2:360227/150‑1 (MQ=255)
                                                                                   cGGAAAACCTGTACATCCGGGGGCTGGCCCAGCCCCATTTGCGGGCGAATATTTCCCACGATAGCACCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgccc                                                                    <  1:32838/150‑1 (MQ=255)
                                                                                       aaaaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGcca                                                                >  1:1006503/1‑150 (MQ=255)
                                                                                         aaCTGTACATACAGTGGCTGGTCCAGTCCCCGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGCGATAATAAACGTAATCACGCCGACGCCCGccacc                                                              >  2:1117153/1‑150 (MQ=255)
                                                                                         aaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGccac                                                               >  2:48097/1‑149 (MQ=255)
                                                                                          aCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGccaccc                                                             >  2:299336/1‑149 (MQ=255)
                                                                                            tGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCcaccac                                                            >  1:679440/1‑149 (MQ=255)
                                                                                               acatacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAg                                                         >  1:854296/1‑149 (MQ=255)
                                                                                                     aGTGGCTGGTCCCGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCCGGGGTCAGAAGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGc                                                  >  1:510917/1‑150 (MQ=255)
                                                                                                     aGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGGAATCACGCAGACGCCCGCCACCACGAGCACCAGc                                                  >  1:397776/1‑150 (MQ=255)
                                                                                                       tGGCTGGTCCGGCCCCGGTTGCGGGCAAATTTTTCCCACGATAGCACCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccc                                                <  2:579848/150‑1 (MQ=255)
                                                                                                       tGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccc                                                >  2:44145/1‑150 (MQ=255)
                                                                                                         gCTGGCCCAGCCCCAGTTGCGGGCAAATTTTTCCCACGATAGCTCCGCGGGCGGGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccccc                                              <  1:820029/150‑1 (MQ=255)
                                                                                                          cTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGcccc                                               >  1:644010/1‑148 (MQ=255)
                                                                                                            ggTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCcag                                            >  2:506565/1‑149 (MQ=255)
                                                                                                               ccAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCagcagc                                        >  2:1150644/1‑150 (MQ=255)
                                                                                                               ccAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCagcagc                                        >  2:735515/1‑150 (MQ=255)
                                                                                                               ccAGTCCCAGTTGCGGGCGAATTTTTTCCAGGATGGCATGGCTGGGGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCagcagc                                        <  1:433657/150‑1 (MQ=255)
                                                                                                                cAGTCCCAGTTGCGGGCAAATATTTTCCACGATAGCATCGCGGGCGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCg                                       <  1:221188/150‑1 (MQ=255)
                                                                                                                     ccAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGc                                  >  2:636789/1‑150 (MQ=255)
                                                                                                                      cAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCg                                 >  2:654859/1‑150 (MQ=255)
                                                                                                                      cAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCg                                 >  1:241747/1‑150 (MQ=255)
                                                                                                                      cAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCg                                 >  1:1064136/1‑150 (MQ=255)
                                                                                                                       aGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCg                                 >  2:143533/1‑149 (MQ=255)
                                                                                                                            cgggCGAATATTTCCCACAAAAGCACCGGGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAAt                           <  1:448586/148‑1 (MQ=255)
                                                                                                                            cgggCAAATTTTTCCCAGGATAGCATCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTaaaa                            <  1:29220/147‑1 (MQ=255)
                                                                                                                              ggaaattatTTCCCGCGATACCATCGCTGGGGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATAc                         <  2:930783/144‑1 (MQ=255)
                                                                                                                              ggCAAATATTTCCCACGATAGCATCGCGGGCGGGCCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATAc                         <  2:417531/150‑1 (MQ=255)
                                                                                                                                cGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTc                       >  2:793101/1‑150 (MQ=255)
                                                                                                                                  aaTTTTTCCCAGGATAGCACCGCGGGGGGGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCa                     <  1:205691/150‑1 (MQ=255)
                                                                                                                                  aaTTTTTCCCACGATAGCATCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCa                     <  1:598297/150‑1 (MQ=255)
                                                                                                                                  aaTATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataat                                                                                          <  1:27864/81‑1 (MQ=255)
                                                                                                                                      ttttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAg                 >  2:385538/1‑150 (MQ=255)
                                                                                                                                       tttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAgg                >  1:954188/1‑150 (MQ=255)
                                                                                                                                       tttCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAgg                >  2:873894/1‑150 (MQ=255)
                                                                                                                                         cccACGATAGCATCGCGGGCGGGCTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGt               <  2:84969/148‑1 (MQ=255)
                                                                                                                                            aCGATAGCATCGCGGGCGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACg           <  1:622856/150‑1 (MQ=255)
                                                                                                                                             cGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGt          >  1:1041487/1‑150 (MQ=255)
                                                                                                                                             cGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGt          >  2:620951/1‑150 (MQ=255)
                                                                                                                                              gATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGt          >  2:952050/1‑149 (MQ=255)
                                                                                                                                                    aTCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTAc   >  2:223428/1‑150 (MQ=255)
                                                                                                                                                    aTCGCTGGCGCGCTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTAc   <  1:978097/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACT  >  CP009273/1554898‑1555196

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