Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,555,008 T→C 8.8% D60G (GAC→GGC)  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,0080TC 8.8% 2325.8 / 29.6 895D60G (GAC→GGC) ddpBD‑ala‑D‑ala transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (437/327);  new base C (5/95);  total (472/423)
Fisher's exact test for biased strand distribution p-value = 3.45e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 906 total aligned reads displayed.

GCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAG  >  CP009273/1554859‑1555157
                                                                                                                                                     |                                                                                                                                                     
gcagcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGt                                                                                                                                                       <  2:981892/150‑1 (MQ=255)
gcagcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGt                                                                                                                                                       <  1:561206/150‑1 (MQ=255)
gcagcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGt                                                                                                                                                       <  1:43529/150‑1 (MQ=255)
 cagcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTc                                                                                                                                                      <  1:2591/150‑1 (MQ=255)
  agcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTcc                                                                                                                                                     <  2:361509/150‑1 (MQ=255)
  agcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTcc                                                                                                                                                     <  2:1063196/150‑1 (MQ=255)
  agcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTcc                                                                                                                                                     <  1:233312/150‑1 (MQ=255)
   gcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTc                                                                                                                                                      >  2:933945/1‑148 (MQ=255)
    cagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAg                                                                                                                                                   <  1:1020134/150‑1 (MQ=255)
    cagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAg                                                                                                                                                   <  2:1047055/150‑1 (MQ=255)
       ggCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTcc                                                                                                                                                >  2:887226/1‑150 (MQ=255)
         cACCAAATGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCa                                                                                                                                              <  1:13318/150‑1 (MQ=255)
          aCCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAg                                                                                                                                             <  1:905789/150‑1 (MQ=255)
          aCCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAg                                                                                                                                             >  1:920711/1‑150 (MQ=255)
           ccAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGt                                                                                                                                            >  2:223043/1‑150 (MQ=255)
            cAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGtt                                                                                                                                           <  1:659354/150‑1 (MQ=255)
             aaaaGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTg                                                                                                                                          <  1:366973/150‑1 (MQ=255)
               aaGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCt                                                                                                                                        >  1:1044661/1‑150 (MQ=255)
                aGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTg                                                                                                                                       <  2:575219/150‑1 (MQ=255)
                  ccAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGc                                                                                                                                     <  2:834982/150‑1 (MQ=255)
                   cAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCg                                                                                                                                    <  2:399413/150‑1 (MQ=255)
                   cAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGGTGCTGGc                                                                                                                                     >  2:110513/1‑149 (MQ=255)
                    aGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCCGGTGCCGGCg                                                                                                                                    >  1:449517/1‑149 (MQ=255)
                     gCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGaa                                                                                                                                  <  1:62401/150‑1 (MQ=255)
                       tCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGTTGCTGGCGAata                                                                                                                                >  2:1023343/1‑150 (MQ=255)
                        ccAGCGTCTCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCCGTTGCTGGCGGatat                                                                                                                               >  2:788607/1‑150 (MQ=255)
                             gTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTcc                                                                                                                          <  1:998336/150‑1 (MQ=255)
                             gTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTcc                                                                                                                          >  2:278310/1‑150 (MQ=255)
                             gTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTcc                                                                                                                          >  2:865993/1‑150 (MQ=255)
                             gTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTc                                                                                                                           >  1:973474/1‑149 (MQ=255)
                              tCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCa                                                                                                                         <  1:1002950/150‑1 (MQ=255)
                                  cGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGAt                                                                                                                     <  2:197671/150‑1 (MQ=255)
                                   gggAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATa                                                                                                                    <  2:31237/150‑1 (MQ=255)
                                     gAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGc                                                                                                                  <  2:432124/150‑1 (MQ=255)
                                      aaaaaTATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCa                                                                                                                 <  2:17754/150‑1 (MQ=255)
                                        aaaTATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATc                                                                                                               <  2:504027/150‑1 (MQ=255)
                                        aaaTATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATc                                                                                                               <  2:13788/150‑1 (MQ=255)
                                         aaTATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCg                                                                                                              <  1:1082010/150‑1 (MQ=255)
                                           tataCGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCt                                                                                                            >  2:1176386/1‑150 (MQ=255)
                                             taCGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTgg                                                                                                          <  2:615814/150‑1 (MQ=255)
                                               cGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgc                                                                                                         <  1:1098746/149‑1 (MQ=255)
                                                 cAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGcgcg                                                                                                      >  1:935046/1‑150 (MQ=255)
                                                    ctctTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTc                                                                                                   <  1:538539/150‑1 (MQ=255)
                                                       ttCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCAcc                                                                                                >  1:725123/1‑150 (MQ=255)
                                                         ccAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCAccc                                                                                               <  2:701320/149‑1 (MQ=255)
                                                           aGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCCGGCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCCCGGGCCCCCGcc                                                                                            >  2:95438/1‑150 (MQ=255)
                                                             cACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                          >  1:1067068/1‑150 (MQ=255)
                                                             cACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAg                                                                                          <  2:472655/150‑1 (MQ=255)
                                                              aCCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTTCTCGATAACAACGCCGGAGCGGCCACCCGcccg                                                                                          >  2:754368/1‑147 (MQ=255)
                                                               ccGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCCGGCGCGGGCACCCGCCAGt                                                                                         >  1:315884/1‑149 (MQ=255)
                                                                        cGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCg                                                                               >  1:70888/1‑150 (MQ=255)
                                                                        cGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGCCCCAGATGCTGGCGAATATTTTTCACGATAGCATTGCTGGCGCGGTCACCCCCCAGCAACCGCGCCg                                                                               >  2:1118051/1‑150 (MQ=255)
                                                                          gTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTTCCGACCCAGTTGCCGGCGAATATTTTCCACGATAGCAGCGCTAGCGCGGGCACCCCCCAGTAAGCACGCCggg                                                                             >  1:859577/1‑149 (MQ=255)
                                                                             cGAATGGATGTTCCCAGGCCCCCATGAAACAGGTCGCTGACGTAGCGTTAAACCTGTACATACGGTGGCTGGTCCGGTCCCAGTTGCTGGCAAATATTTCCCACGATAGCATCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATc                                                                           <  2:384012/149‑1 (MQ=255)
                                                                              gAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGATCAc                                                                         >  2:21124/1‑150 (MQ=255)
                                                                              gAATGGATGTCCCCAGGTCACCAGAAAACAGGTCGCTGACGTAGGGGAAAAACTGTACATACAGGGGCTGGTCCAGCCCCAGTTGCGGGCGAATATTTCCCACGATAGCATCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCAc                                                                         <  2:713490/150‑1 (MQ=255)
                                                                                 tGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAg                                                                      >  2:73398/1‑150 (MQ=255)
                                                                                 tGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGAACACCAg                                                                      >  2:113711/1‑150 (MQ=255)
                                                                                  ggATGTCCCCAGGTCCCCATGAAACAGGCCGCTGACGTAGCGGTAAAACTGTACATACGGTGGCTGGCCCAGTCCCAGTTGCGGGAAAATATTTCCCACGATAGCACCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAgg                                                                     <  2:1085648/150‑1 (MQ=255)
                                                                                   gATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAgg                                                                     >  1:810152/1‑149 (MQ=255)
                                                                                   gATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGGAACCGCGCCGGATCACCAggg                                                                    >  2:120573/1‑150 (MQ=255)
                                                                                    aTGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGa                                                                   >  1:603073/1‑150 (MQ=255)
                                                                                          ccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTTCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGAt                                                             >  2:652078/1‑150 (MQ=255)
                                                                                           cAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATg                                                            >  1:289630/1‑150 (MQ=255)
                                                                                             ggTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGCTCAGTCCCCGTTGCTGGCGAAAATTTTTCACGATAGCATCGCTGGCGCGGTCCCCCGCCAGTCACCGCGCCGGATCACCAGGGATCAGATGc                                                           >  2:633828/1‑149 (MQ=255)
                                                                                               tCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGc                                                           >  1:475535/1‑147 (MQ=255)
                                                                                                    aTGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGGGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                    >  1:213952/1‑149 (MQ=255)
                                                                                                    aTGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataa                                                    >  1:507571/1‑149 (MQ=255)
                                                                                                     tGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACGCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGataat                                                   >  1:753852/1‑149 (MQ=255)
                                                                                                      gAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGCGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCCGGCGCGGTCACCCGCCACTAACCGCGCCGGATCACCAGGGATTAGATGCGAGataata                                                  >  1:787460/1‑149 (MQ=255)
                                                                                                          cAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGCCCCAGTTGCTGGCAAATTTTTCCCACGATAGCATCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGt                                             <  1:159845/150‑1 (MQ=255)
                                                                                                            ggTCGCTGCCGTAGCGGAAAACCTGTACATACGGTGGCTGGCCCAGCCCCGGTTGCTGGCAAATTTTTCCCACGATAGCACCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTaa                                           <  2:597808/150‑1 (MQ=255)
                                                                                                                gCTGACGTGGCGGTAAACCTGTACATCAAGTGGCTGGTCCGGCCCCAGTTGCGGGCAAATATTTCCCACGATAGCACCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCAc                                       <  1:887274/150‑1 (MQ=255)
                                                                                                                   gACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCCGTTGCTGGCGAATATTTTCCACGATAGCATCGCTCGCGCGGGCACCCGCCAGTAACCGCGCCCGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACg                                      >  2:899199/1‑148 (MQ=255)
                                                                                                                   gACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACg                                      >  1:403773/1‑148 (MQ=255)
                                                                                                                   gACGTAGCGGAAAAACTGTACAGACAGTGGCTGTTCCAGTCCCAGTTGCTGGCGAATATTTCCCACGATAGCACCGCTGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                    <  1:643532/150‑1 (MQ=255)
                                                                                                                    aCGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACg                                      >  1:333329/1‑147 (MQ=255)
                                                                                                                      gTAGGGGAAAAGCTATACATGCGGTGGCGGGTCAAGCCCAAGTTGCTGGGGAATATTTTCCACAATAGCACCGCGGGCGGGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAc                                 <  1:687579/150‑1 (MQ=255)
                                                                                                                      gTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAc                                 >  2:412338/1‑150 (MQ=255)
                                                                                                                      gTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAc                                 >  1:579848/1‑150 (MQ=255)
                                                                                                                      gTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGAAAGCATCCCTGGCGCGGTCACCCCCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAc                                 >  1:832573/1‑150 (MQ=255)
                                                                                                                       tAGCGATAAAACTTTACATACAGTGGCTCGTCCGGCCCCAGTTGCGGGCAAATTTTTCCCACGATACCATCGCTGGCGGGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcg                                <  2:252639/150‑1 (MQ=255)
                                                                                                                           ggTAAAACTGTACATACAGTGGCTGGTCCAGTTCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccg                            >  1:1128685/1‑150 (MQ=255)
                                                                                                                           ggTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccg                            >  1:1146781/1‑150 (MQ=255)
                                                                                                                            gTAAAACTGTACATACAGTGGCTGGTCCAGTTCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTTGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGGTCCGATGCGAGAAAATAAACGTAATCACCCCGAcccccgc                           >  2:685215/1‑150 (MQ=255)
                                                                                                                             tAAAACTGTACATACAGGGGCTGGCCCAGTCCCAGTTGCGGGCAAATTTTTCCCACAATAGCACCGCTGGCGGGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccgcc                          <  1:948482/150‑1 (MQ=255)
                                                                                                                              aaaaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgcccgc                           >  1:489450/1‑148 (MQ=255)
                                                                                                                               aaaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAAAGTAATAACGCAGACGCCCGccac                        >  2:552253/1‑150 (MQ=255)
                                                                                                                                aaCTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCCGGTCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGccacc                       >  2:1585/1‑150 (MQ=255)
                                                                                                                                  cTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCcaccac                     >  2:816554/1‑150 (MQ=255)
                                                                                                                                  cTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCcaccac                     >  1:852826/1‑150 (MQ=255)
                                                                                                                                     tacatacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAg                  >  1:675166/1‑150 (MQ=255)
                                                                                                                                     taaatacgGGGGCTGGCCCAGCCCCGGTTGCGGGCAAATTTTCCCCACGATAGCCCCGCGGGCGGGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGa                   <  2:520006/146‑1 (MQ=255)
                                                                                                                                        atacaGTGGCTGGGCCAGTCCCAGTTGCTGGCGAATATTTTTCACGATAGCAACGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCAc               >  2:805499/1‑150 (MQ=255)
                                                                                                                                        atacaGTGGCTGGCCCAGCCCCGGTTGCTGGCAAATATTTCCCACGATAGCCTCGCTGGCGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCAc               <  2:448661/150‑1 (MQ=255)
                                                                                                                                           cggTGGCTGGCCCAGCCCCGGTGGCTGGCAAATATTTTCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCa             <  1:517523/147‑1 (MQ=255)
                                                                                                                                              tGGCTGCCCCAGCCCCAGTTGCGGGCGAATATTTCCCACGATACCATCGCGGCCGCGTTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccc         <  1:190365/150‑1 (MQ=255)
                                                                                                                                               ggtgggTCAAGCCCCATTTGCGGGCAAATTTTTCCCACGATAGCTCCGCGGGCGGGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGcccc        <  1:11998/146‑1 (MQ=255)
                                                                                                                                                   ggTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCcag     >  1:385641/1‑149 (MQ=255)
                                                                                                                                                    ctCCAGTCCCGGTCGCGGGCAATTTTTTCCCACAAAGCCACCGCGGGCCCGCTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCcagc    <  2:335051/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAG  >  CP009273/1554859‑1555157

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