Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,151,750 G→T 91.2% 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,7500GT91.2% 1065.1 / 92.4 367R178L (CGC→CTC) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (14/18);  new base T (157/177);  total (172/195)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

Only 100 of 369 total aligned reads displayed.

TGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCGG  >  CP009273/1151603‑1151896
                                                                                                                                                   |                                                                                                                                                  
tGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCActcg                                                                                                                                                  <  1:1111851/150‑1 (MQ=255)
   aaaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAg                                                                                                                                               >  2:1392821/1‑150 (MQ=255)
    aaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGt                                                                                                                                              >  1:3830623/1‑150 (MQ=255)
    aaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGt                                                                                                                                              <  2:3849778/150‑1 (MQ=255)
     aCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTg                                                                                                                                             <  2:4390507/150‑1 (MQ=255)
      cGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGa                                                                                                                                            <  2:2297003/150‑1 (MQ=255)
      cGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGa                                                                                                                                            <  1:936391/150‑1 (MQ=255)
         ttGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACaa                                                                                                                                         <  2:2888302/150‑1 (MQ=255)
          tGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAt                                                                                                                                        <  2:3998049/150‑1 (MQ=255)
          tGAAGAGCCAACAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAt                                                                                                                                        <  2:4274447/150‑1 (MQ=255)
               agCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTcac                                                                                                                                   <  2:1882383/150‑1 (MQ=255)
               agCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTcac                                                                                                                                   >  1:3043288/1‑150 (MQ=255)
               agCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTcac                                                                                                                                   >  2:1270782/1‑150 (MQ=255)
                       gcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGGACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCa                                                                                                                           >  2:11567/1‑150 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                          >  2:3078866/1‑150 (MQ=255)
                        cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                                                          >  2:1311712/1‑150 (MQ=255)
                         agAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAtt                                                                                                                         >  2:2594392/1‑150 (MQ=255)
                         agAAACTTGGCTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCAtt                                                                                                                         <  2:2409520/150‑1 (MQ=255)
                          gAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACAGTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTa                                                                                                                        <  2:2028786/150‑1 (MQ=255)
                           aaaCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTAc                                                                                                                       >  1:2590974/1‑150 (MQ=255)
                             aCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACtt                                                                                                                     <  2:1569581/150‑1 (MQ=255)
                              cTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTg                                                                                                                    >  1:3795242/1‑150 (MQ=255)
                                     tttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCAtt                                                                                                             <  1:3097250/150‑1 (MQ=255)
                                      ttCCTGGCGACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTg                                                                                                            <  1:3519051/150‑1 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  2:4734788/1‑150 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         <  1:4714661/150‑1 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         <  1:3318572/150‑1 (MQ=255)
                                         cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                                                         >  2:301107/1‑150 (MQ=255)
                                          tGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCt                                                                                                        >  2:1960392/1‑150 (MQ=255)
                                           ggCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCtt                                                                                                       <  2:4618439/150‑1 (MQ=255)
                                           ggCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCtt                                                                                                       >  1:4456628/1‑150 (MQ=255)
                                            gCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTa                                                                                                      <  1:373105/150‑1 (MQ=255)
                                              tACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAAGTCACAGGAAGCAATACTTTCCGCATTGCGCTTAAg                                                                                                    >  1:2484448/1‑150 (MQ=255)
                                                cccGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAgca                                                                                                  <  2:1120587/150‑2 (MQ=255)
                                                  cgcgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCCCAGGATGCATTACTTGCCGCATTGCGCTTAAGCGcc                                                                                                >  2:376165/1‑150 (MQ=255)
                                                  cgcgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGcc                                                                                                >  2:2746958/1‑150 (MQ=255)
                                                      aGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGtt                                                                                            <  2:3552344/150‑1 (MQ=255)
                                                       gCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTc                                                                                           <  2:3831558/150‑1 (MQ=255)
                                                          tGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGcc                                                                                        <  1:2883640/150‑1 (MQ=255)
                                                           gAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCt                                                                                       <  2:1285715/150‑1 (MQ=255)
                                                               gTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGgcg                                                                                   <  2:4449714/150‑1 (MQ=255)
                                                                    cTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAATCGCCGGTTCGCCTGGCGCGGcc                                                                              >  2:1282811/1‑149 (MQ=255)
                                                                     tGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCAc                                                                             >  2:760303/1‑150 (MQ=255)
                                                                     tGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCAc                                                                             >  2:1911641/1‑150 (MQ=255)
                                                                       gCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAAAGTCACAGGATGCATTACTTGCCGCATTGCGCTTAATCGCCGGTTCGCCTGGCGCGGCACTg                                                                           >  1:2103640/1‑150 (MQ=255)
                                                                        cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTgg                                                                          <  1:360991/150‑1 (MQ=255)
                                                                        cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTgg                                                                          >  2:4193111/1‑150 (MQ=255)
                                                                         aaCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGc                                                                         <  1:825685/150‑1 (MQ=255)
                                                                                 gTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTc                                                                 <  2:2183526/150‑1 (MQ=255)
                                                                                   aGTCGTAGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAg                                                               >  1:2300370/1‑150 (MQ=255)
                                                                                    gTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAgg                                                              >  1:670031/1‑150 (MQ=255)
                                                                                    gTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCCGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAgg                                                              <  2:3448560/150‑1 (MQ=255)
                                                                                       gTTGTCGGTTACCTTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgag                                                           >  1:2087175/1‑150 (MQ=255)
                                                                                        ttGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgaga                                                          <  1:2680500/150‑1 (MQ=255)
                                                                                         tGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGAt                                                         <  1:3927798/150‑1 (MQ=255)
                                                                                            cGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAAc                                                      >  2:1265329/1‑150 (MQ=255)
                                                                                              gTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTg                                                    <  2:166118/150‑1 (MQ=255)
                                                                                              gTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTg                                                    >  2:1843046/1‑150 (MQ=255)
                                                                                                tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc                                                  >  2:4492878/1‑150 (MQ=255)
                                                                                                tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc                                                  <  1:3134380/150‑1 (MQ=255)
                                                                                                 aCATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCa                                                 <  1:1540087/150‑1 (MQ=255)
                                                                                                  cATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAg                                                <  2:2676728/150‑1 (MQ=255)
                                                                                                     tACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCt                                             >  2:4366834/1‑150 (MQ=255)
                                                                                                      aCCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTc                                            <  1:2179764/150‑1 (MQ=255)
                                                                                                       ccTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCg                                           >  1:1005342/1‑150 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          <  2:4161527/150‑1 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          <  2:3331306/150‑1 (MQ=255)
                                                                                                        cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt                                          >  2:4356463/1‑150 (MQ=255)
                                                                                                         ttGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTg                                         <  2:927497/150‑1 (MQ=255)
                                                                                                            cgccgccgcCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGaaa                                      <  1:1703570/150‑1 (MQ=255)
                                                                                                              ccgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACa                                    >  2:2101154/1‑150 (MQ=255)
                                                                                                              ccgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACa                                    <  2:1431849/150‑1 (MQ=255)
                                                                                                              ccgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACa                                    <  2:4082334/150‑1 (MQ=255)
                                                                                                               cgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACAt                                   >  2:537717/1‑150 (MQ=255)
                                                                                                               cgccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACAt                                   <  1:4166303/150‑1 (MQ=255)
                                                                                                                gccgccgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACAtt                                  <  1:4569748/150‑1 (MQ=255)
                                                                                                                    ccgGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCTCAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTgtgt                              <  1:4345835/150‑1 (MQ=255)
                                                                                                                     cgGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTc                             >  2:336007/1‑150 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            >  1:2965808/1‑150 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            <  1:1859370/150‑1 (MQ=255)
                                                                                                                      ggAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCa                            <  1:4473780/150‑1 (MQ=255)
                                                                                                                       gAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAg                           >  2:588289/1‑150 (MQ=255)
                                                                                                                       gAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAg                           <  1:1818298/150‑1 (MQ=255)
                                                                                                                        aaCAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAgg                          >  1:3421563/1‑150 (MQ=255)
                                                                                                                          cAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCg                        >  2:873867/1‑150 (MQ=255)
                                                                                                                           aGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGt                       <  2:2377269/150‑1 (MQ=255)
                                                                                                                             tACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTg                     <  1:4640460/150‑1 (MQ=255)
                                                                                                                              aCGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTgg                    <  2:4556279/150‑1 (MQ=255)
                                                                                                                                gCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCa                  <  1:532654/150‑1 (MQ=255)
                                                                                                                                 ccGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCat                 <  1:1266530/150‑1 (MQ=255)
                                                                                                                                  cGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCata                >  1:1786281/1‑150 (MQ=255)
                                                                                                                                  cGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCata                >  2:3243368/1‑150 (MQ=255)
                                                                                                                                     gACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAg             <  1:1329326/150‑1 (MQ=255)
                                                                                                                                       ccTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCg           <  1:2799045/150‑1 (MQ=255)
                                                                                                                                           gCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGcc       <  1:1179064/150‑1 (MQ=255)
                                                                                                                                            cTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCa      <  1:696293/150‑1 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt     >  1:1676025/1‑150 (MQ=255)
                                                                                                                                             tttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCAt     <  2:3162129/150‑1 (MQ=255)
                                                                                                                                              ttCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATc    <  1:3237143/150‑1 (MQ=255)
                                                                                                                                                cACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAACATTGTGTCAGGCGTTGGCATATAGCGTGCCATCgg  <  2:3176535/150‑1 (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: