Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,151,750 G→T 90.9% R178L (CGC→CTC)  holB → DNA polymerase III, delta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,151,7500GT90.9% 825.0 / 76.4 298R178L (CGC→CTC) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (8/19);  new base T (142/127);  total (152/146)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

Only 100 of 300 total aligned reads displayed.

TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGG  >  CP009273/1151603‑1151896
                                                                                                                                                   |                                                                                                                                                  
  aaaaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGaa                                                                                                                                                <  2:702653/150‑1 (MQ=255)
    aaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGt                                                                                                                                              <  1:3271097/150‑1 (MQ=255)
         ttGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACaa                                                                                                                                         >  2:1172873/1‑150 (MQ=255)
           gAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATg                                                                                                                                       <  2:1432746/150‑1 (MQ=255)
           gAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATg                                                                                                                                       >  1:1922057/1‑150 (MQ=255)
              gagCCACCAGCAGAAACTGGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTCACTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTca                                                                                                                                    >  2:395985/1‑150 (MQ=255)
               agCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTcac                                                                                                                                   <  1:3765147/150‑1 (MQ=255)
                  caccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAgg                                                                                                                                >  2:1827265/1‑150 (MQ=255)
                  caccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAgg                                                                                                                                >  2:3937463/1‑150 (MQ=255)
                   accaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGa                                                                                                                               <  2:67819/150‑1 (MQ=255)
                     cagcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATg                                                                                                                             >  2:3239358/1‑150 (MQ=255)
                      agcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGc                                                                                                                            <  2:4124588/150‑1 (MQ=255)
                       gcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCa                                                                                                                           <  1:1291706/150‑1 (MQ=255)
                       gcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCa                                                                                                                           <  2:3090416/150‑1 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                          <  2:3025453/150‑1 (MQ=255)
                         agAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAtt                                                                                                                         >  1:1556975/1‑150 (MQ=255)
                         agAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAtt                                                                                                                         >  2:3458899/1‑150 (MQ=255)
                          gAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTa                                                                                                                        >  1:646011/1‑150 (MQ=255)
                            aaCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTAAGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACt                                                                                                                      >  1:4279608/1‑150 (MQ=255)
                                tGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGcc                                                                                                                  <  1:305451/150‑1 (MQ=255)
                                tGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGcc                                                                                                                  >  1:4101756/1‑150 (MQ=255)
                                 ggTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCg                                                                                                                 <  1:634741/150‑1 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  2:2809617/1‑150 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  1:1997698/1‑150 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  2:546472/1‑150 (MQ=255)
                                          tGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCt                                                                                                        >  2:189799/1‑150 (MQ=255)
                                           ggCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCtt                                                                                                       <  2:4257636/150‑1 (MQ=255)
                                            gCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTa                                                                                                      <  1:1099167/150‑1 (MQ=255)
                                            gCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTa                                                                                                      <  2:910017/150‑1 (MQ=255)
                                             cTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTaa                                                                                                     >  2:3540045/1‑150 (MQ=255)
                                             cTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTaa                                                                                                     >  2:112150/1‑150 (MQ=255)
                                              tACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAg                                                                                                    >  1:2074188/1‑150 (MQ=255)
                                                cccGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAgcg                                                                                                  <  1:2964889/150‑1 (MQ=255)
                                                  cgcgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGcc                                                                                                <  2:3705015/150‑1 (MQ=255)
                                                    cgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGGACAGTACGCCGTGACCGGGCTTTCACACGAAGTGACAATGTCACAGGCTGCATTACTTGCCGCCTTTCGCTTAAGCGcccg                                                                                              >  1:2314147/1‑148 (MQ=255)
                                                       gCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTc                                                                                           <  1:762688/150‑1 (MQ=255)
                                                          tGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGcc                                                                                        <  2:4112301/150‑1 (MQ=255)
                                                           gAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCt                                                                                       <  1:2361127/150‑1 (MQ=255)
                                                            aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                                                      <  1:92103/150‑1 (MQ=255)
                                                            aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                                                      >  2:1303177/1‑150 (MQ=255)
                                                             aCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTgg                                                                                     <  1:47767/150‑1 (MQ=255)
                                                              cGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGgc                                                                                    >  2:1240141/1‑150 (MQ=255)
                                                              cGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGgc                                                                                    <  2:443496/150‑1 (MQ=255)
                                                              cGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCACTTAAGCGCCGGTTCGCCTGgc                                                                                    <  2:212736/150‑1 (MQ=255)
                                                                 ttACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGcgcg                                                                                 <  1:2531434/150‑1 (MQ=255)
                                                                  tACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTAAAGCGCCGGTTCGCCTGGCGCgg                                                                                >  2:3167176/1‑150 (MQ=255)
                                                                    cTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCa                                                                              >  2:3828963/1‑150 (MQ=255)
                                                                     tGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCAc                                                                             >  2:2977910/1‑150 (MQ=255)
                                                                        cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAATCGCCGGTTCGCCTGGCGCGGCACTgg                                                                          >  1:3255570/1‑150 (MQ=255)
                                                                          aCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCg                                                                        <  2:84566/150‑1 (MQ=255)
                                                                            aTTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATTCATTACTTGCCGCAATGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACAGGCgtt                                                                      >  1:2889165/1‑150 (MQ=255)
                                                                            aTTACGTAGTCGTTGTCGGTTACAATACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTAACTCGAAGTTAAAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAATCGCCGGTTAGCCTGGCGCGGCACTGGCgtt                                                                      >  1:3474713/1‑150 (MQ=255)
                                                                                 gTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTc                                                                 <  2:1982462/150‑1 (MQ=255)
                                                                                   aGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAg                                                               >  2:1064513/1‑150 (MQ=255)
                                                                                       gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTAAGGgag                                                           >  2:2385239/1‑150 (MQ=255)
                                                                                       gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgag                                                           <  1:1765413/150‑1 (MQ=255)
                                                                                       gTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgag                                                           >  2:1551060/1‑150 (MQ=255)
                                                                                         tGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGAt                                                         <  2:3910800/150‑1 (MQ=255)
                                                                                          gTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATa                                                        <  1:697715/150‑1 (MQ=255)
                                                                                             ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACt                                                     <  1:3239358/150‑1 (MQ=255)
                                                                                             ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACt                                                     >  2:3404013/1‑150 (MQ=255)
                                                                                             ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACt                                                     <  2:1604227/150‑1 (MQ=255)
                                                                                                tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc                                                  <  1:3559752/150‑1 (MQ=255)
                                                                                                tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc                                                  >  1:1948030/1‑150 (MQ=255)
                                                                                                   aTTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAgg                                               >  2:1347549/1‑150 (MQ=255)
                                                                                                    ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGc                                              >  2:4239200/1‑150 (MQ=255)
                                                                                                    ttACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGc                                              >  1:4064594/1‑150 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          >  2:1233692/1‑150 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          >  2:2597201/1‑150 (MQ=255)
                                                                                                              ccgccgccgGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACa                                    <  2:2566967/150‑1 (MQ=255)
                                                                                                                 ccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATtg                                 <  2:3332177/150‑1 (MQ=255)
                                                                                                                 ccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATtg                                 <  2:4023556/150‑1 (MQ=255)
                                                                                                                     cgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTc                             <  1:1291093/150‑1 (MQ=255)
                                                                                                                     cgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTc                             <  2:1521109/150‑1 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            <  1:2500932/150‑1 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTAGTGAAACATTGTGTCa                            >  1:3747767/1‑150 (MQ=255)
                                                                                                                        aaCAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAgg                          >  1:3315297/1‑150 (MQ=255)
                                                                                                                        aaCAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAgg                          >  1:3145893/1‑150 (MQ=255)
                                                                                                                          cAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCg                        >  2:3388030/1‑150 (MQ=255)
                                                                                                                          cAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCg                        <  1:3611680/150‑1 (MQ=255)
                                                                                                                          cAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCg                        <  1:708685/150‑1 (MQ=255)
                                                                                                                           aGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGt                       >  2:2022140/1‑150 (MQ=255)
                                                                                                                            gTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGtt                      <  1:1254231/150‑1 (MQ=255)
                                                                                                                              aCGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTgg                    >  2:815100/1‑150 (MQ=255)
                                                                                                                                gCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCa                  >  2:2347122/1‑150 (MQ=255)
                                                                                                                                   gTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCatat               <  1:1240103/150‑1 (MQ=255)
                                                                                                                                     gACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAg             <  2:2494374/150‑1 (MQ=255)
                                                                                                                                     gACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAg             >  2:1546516/1‑150 (MQ=255)
                                                                                                                                     gACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAg             <  2:3101374/150‑1 (MQ=255)
                                                                                                                                     gACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAg             <  1:830968/150‑1 (MQ=255)
                                                                                                                                        ctcgCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGt          >  1:2902580/4‑150 (MQ=255)
                                                                                                                                         tGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTg         <  1:3210113/150‑1 (MQ=255)
                                                                                                                                         tGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTg         <  1:2809478/150‑1 (MQ=255)
                                                                                                                                           gCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGcc       >  2:1685393/1‑150 (MQ=255)
                                                                                                                                            cTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCTCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCa      >  2:2753887/1‑150 (MQ=255)
                                                                                                                                            cTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCa      >  1:2240291/1‑150 (MQ=255)
                                                                                                                                            cTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTAGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCa      <  1:2177892/150‑1 (MQ=255)
                                                                                                                                            cTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCa      <  1:3937463/150‑1 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt     >  1:2981823/1‑150 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt     >  1:1030685/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGG  >  CP009273/1151603‑1151896

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: