Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,151,711 C→A 6.3% A165E (GCG→GAG)  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,7110CA 6.3% 1429.3 / 44.3 378A165E (GCG→GAG) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (173/180);  new base A (16/8);  total (190/188)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 394 total aligned reads displayed.

CGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAA  >  CP009273/1151562‑1151860
                                                                                                                                                     |                                                                                                                                                     
cGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTgc                                                                                                                                                       >  2:3436065/1‑150 (MQ=255)
 gATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTgcg                                                                                                                                                      <  2:537196/150‑1 (MQ=255)
 gATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTgcg                                                                                                                                                      <  2:2803360/150‑1 (MQ=255)
 gATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTgcg                                                                                                                                                      >  1:2064785/1‑150 (MQ=255)
   tgctgcCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGcgcc                                                                                                                                                    <  2:2342164/150‑1 (MQ=255)
   tgctgcCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGcgcc                                                                                                                                                    >  2:1386851/1‑150 (MQ=255)
     ctgcCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGcgccgc                                                                                                                                                  >  2:1317323/1‑150 (MQ=255)
       gcCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGcgccgccg                                                                                                                                                <  1:4566853/150‑1 (MQ=255)
         cTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCgccgccgcc                                                                                                                                              <  1:4332159/150‑1 (MQ=255)
          ttACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGccgccgccg                                                                                                                                             >  2:933970/1‑150 (MQ=255)
           tACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCgg                                                                                                                                            >  2:2430938/1‑150 (MQ=255)
             cTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGaa                                                                                                                                          >  1:541683/1‑150 (MQ=255)
               aaCCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACa                                                                                                                                        >  2:3423269/1‑150 (MQ=255)
                aCCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAg                                                                                                                                       <  2:3802392/150‑1 (MQ=255)
                  cGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTa                                                                                                                                     <  2:850343/150‑1 (MQ=255)
                  cGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTa                                                                                                                                     >  1:3162129/1‑150 (MQ=255)
                  cGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTa                                                                                                                                     >  1:3831558/1‑150 (MQ=255)
                   gACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTAc                                                                                                                                    >  1:2183526/1‑150 (MQ=255)
                    aCGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACg                                                                                                                                   >  2:796486/1‑150 (MQ=255)
                      gccgcGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGcc                                                                                                                                 <  2:2261196/150‑1 (MQ=255)
                             cTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCt                                                                                                                          >  1:2377269/1‑150 (MQ=255)
                              tAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGAcccg                                                                                                                         >  1:1215266/1‑148 (MQ=255)
                              tAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTg                                                                                                                         >  2:3259981/1‑150 (MQ=255)
                              tAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTg                                                                                                                         >  1:1524709/1‑150 (MQ=255)
                              tAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTg                                                                                                                         >  1:3911568/1‑150 (MQ=255)
                                 cGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCt                                                                                                                      >  2:3366486/1‑150 (MQ=255)
                                 cGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCt                                                                                                                      <  2:4218748/150‑1 (MQ=255)
                                  gCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCtt                                                                                                                     <  2:1822813/150‑1 (MQ=255)
                                  gCATTGCTGAAAACGCGTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCtt                                                                                                                     <  2:4284156/150‑1 (MQ=255)
                                   cATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCttt                                                                                                                    <  2:1657305/150‑1 (MQ=255)
                                    aTTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCATGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTc                                                                                                                   >  1:1387677/1‑150 (MQ=255)
                                      tGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAc                                                                                                                 >  1:741493/1‑150 (MQ=255)
                                      tGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCAc                                                                                                                 <  1:463224/150‑1 (MQ=255)
                                         tGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCActcg                                                                                                              <  1:1111851/150‑1 (MQ=255)
                                            aaaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAg                                                                                                           <  1:1953523/150‑1 (MQ=255)
                                            aaaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAg                                                                                                           <  1:2659421/150‑1 (MQ=255)
                                             aaCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGt                                                                                                          >  2:277502/1‑150 (MQ=255)
                                              aCGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTg                                                                                                         <  2:4390507/150‑1 (MQ=255)
                                               cGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGa                                                                                                        <  1:936391/150‑1 (MQ=255)
                                               cGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGa                                                                                                        <  1:3479719/150‑1 (MQ=255)
                                                 cTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACa                                                                                                      >  2:4166303/1‑150 (MQ=255)
                                                  ttGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACaa                                                                                                     >  1:4697155/1‑150 (MQ=255)
                                                   tGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCCTTCACTCGAAGTGACAAt                                                                                                    <  2:888139/150‑1 (MQ=255)
                                                      agagCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTc                                                                                                 <  2:1043950/150‑1 (MQ=255)
                                                       gagCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTca                                                                                                <  2:4619734/150‑1 (MQ=255)
                                                        agCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTcac                                                                                               >  2:1270782/1‑150 (MQ=255)
                                                          ccaccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAg                                                                                             <  1:1763087/150‑1 (MQ=255)
                                                          ccaccaGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAg                                                                                             >  2:1404044/1‑150 (MQ=255)
                                                                gcagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGGACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCa                                                                                       >  2:11567/1‑150 (MQ=255)
                                                                 cagAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAt                                                                                      >  2:3078866/1‑150 (MQ=255)
                                                                  agAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAtt                                                                                     >  2:2594392/1‑150 (MQ=255)
                                                                     aaCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACt                                                                                  <  2:1883263/150‑1 (MQ=255)
                                                                      aCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACtt                                                                                 >  2:3921704/1‑150 (MQ=255)
                                                                      aCTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACtt                                                                                 >  2:15037/1‑150 (MQ=255)
                                                                       cTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTg                                                                                <  1:1804270/150‑1 (MQ=255)
                                                                       cTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTg                                                                                >  1:3795266/1‑150 (MQ=255)
                                                                           gTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGc                                                                            >  2:2395090/1‑150 (MQ=255)
                                                                              tttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCAtt                                                                         <  1:4446765/150‑1 (MQ=255)
                                                                              tttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCAtt                                                                         >  2:128708/1‑150 (MQ=255)
                                                                              tttCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTAACTTTCGCCGCCGCCCGAACAGTAAGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTTCCGCAtt                                                                         >  1:2561410/1‑150 (MQ=255)
                                                                                  cTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTgcgc                                                                     <  1:3318572/150‑1 (MQ=255)
                                                                                    ggCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCtt                                                                   >  1:4458577/1‑150 (MQ=255)
                                                                                     gCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTa                                                                  <  1:373105/150‑1 (MQ=255)
                                                                                      cTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTaa                                                                 >  2:723170/1‑150 (MQ=255)
                                                                                         cccGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAgca                                                              <  2:3510092/150‑2 (MQ=255)
                                                                                           cgcgAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTTCGCCCCCGCCGGAACAGTACGCCGTGTCCTGGCTTTCACTCGAAGTGAAAAAGTAACACGATTCAATACTTTCCGCCTTGCGCGTAAGCGcc                                                            >  1:418091/1‑150 (MQ=255)
                                                                                           cgcgAGCCTGAACGTTTACTGGCAACATTACGGAGTCGTTGTCGGTTACATGACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTGGCCGCATTGCGCTTAAGCGcc                                                            <  2:3694758/150‑1 (MQ=255)
                                                                                              gAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCCCAGGATGCATTACTTGCCGCAGTGCGCTTAAGCGCCGGt                                                         >  2:4738786/1‑150 (MQ=255)
                                                                                                gCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTc                                                       <  2:3831562/150‑1 (MQ=255)
                                                                                                   tGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGcc                                                    <  1:2883640/150‑1 (MQ=255)
                                                                                                    gAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCt                                                   <  2:1285656/150‑1 (MQ=255)
                                                                                                     aaCGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTg                                                  <  2:2973713/150‑1 (MQ=255)
                                                                                                       cGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGgc                                                >  2:4002501/1‑150 (MQ=255)
                                                                                                        gTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGgcg                                               <  2:4449714/150‑1 (MQ=255)
                                                                                                                gCAACATTACGTAGTCGTTGTCGGTTACAATAACTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCAATACTTTCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCAcag                                       >  1:2103596/1‑148 (MQ=255)
                                                                                                                 cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTgg                                      <  2:2688602/150‑1 (MQ=255)
                                                                                                                 cAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTgg                                      <  1:360991/150‑1 (MQ=255)
                                                                                                                  aaCATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGc                                     <  1:825685/150‑1 (MQ=255)
                                                                                                                      ttACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCgttg                                 <  1:3366486/150‑1 (MQ=255)
                                                                                                                          gTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTc                             <  2:4142871/150‑1 (MQ=255)
                                                                                                                          gTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTc                             <  1:2298370/150‑1 (MQ=255)
                                                                                                                          gTAGTCGTTGTCGGTTAAATTACCATGCGCCGCCGCCGGAACAGTACGCCCTGACCTGGCTTTCACTCGAAGTGACAAAGTCACAGGATGCATTACTTGCCGCAATGCGCTTAAGCGCCGGTTCGCCAGGCGCGGCACTGGCGTTTTTTc                             >  2:1985924/1‑150 (MQ=255)
                                                                                                                             gTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAgg                          >  1:670031/1‑150 (MQ=255)
                                                                                                                                 ttGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGgaga                      >  2:3156013/1‑150 (MQ=255)
                                                                                                                                     cGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAAc                  >  2:1265329/1‑150 (MQ=255)
                                                                                                                                      ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACt                 >  1:1617815/1‑150 (MQ=255)
                                                                                                                                      ggTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACt                 <  1:128708/150‑1 (MQ=255)
                                                                                                                                       gTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTg                >  2:2146082/1‑150 (MQ=255)
                                                                                                                                       gTTACATTACCTTGAGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTg                >  2:1843372/1‑150 (MQ=255)
                                                                                                                                         tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc              >  2:4493744/1‑150 (MQ=255)
                                                                                                                                         tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc              <  1:3134352/150‑1 (MQ=255)
                                                                                                                                         tACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGc              >  1:2465675/1‑150 (MQ=255)
                                                                                                                                            aTTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAgg           >  1:3059689/1‑150 (MQ=255)
                                                                                                                                              tACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCt         >  2:4366834/1‑150 (MQ=255)
                                                                                                                                               aCCATGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTc        >  2:1947100/1‑150 (MQ=255)
                                                                                                                                                ccTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCg       >  1:2843688/1‑150 (MQ=255)
                                                                                                                                                ccTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCg       >  1:1005342/1‑150 (MQ=255)
                                                                                                                                                 cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt      >  2:4356463/1‑150 (MQ=255)
                                                                                                                                                 cTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGt      <  2:3331306/150‑1 (MQ=255)
                                                                                                                                                  ttGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACTCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTg     <  2:2993797/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCCTGAACGTTTACTGGCAACATTACGTAGTCGTTGTCGGTTACATTACCTTGCGCCGCCGCCGGAACAGTACGCCGTGACCTGGCTTTCACGCGAAGTGACAATGTCACAGGATGCATTACTTGCCGCATTGCGCTTAAGCGCCGGTTCGCCTGGCGCGGCACTGGCGTTGTTTCAGGGAGATAACTGGCAGGCTCGTGAAA  >  CP009273/1151562‑1151860

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: