Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,176,823 G→T 85.7% D516Y (GAC→TAC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,176,8230GT85.7% 681.9 / 111.3 258D516Y (GAC→TAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base T (92/129);  total (113/145)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 261 total aligned reads displayed.

TGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTAGCG  >  CP009273/4176674‑4176972
                                                                                                                                                     |                                                                                                                                                     
tGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTg                                                                                                                                                       <  2:336407/150‑1 (MQ=255)
  gCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGAc                                                                                                                                                     >  2:2231955/1‑150 (MQ=255)
    tGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACtg                                                                                                                                                   <  2:1503297/150‑1 (MQ=255)
         aCGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTa                                                                                                                                              <  1:1706298/150‑1 (MQ=255)
         aCGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTa                                                                                                                                              <  1:1706362/150‑1 (MQ=255)
            tACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACg                                                                                                                                           <  1:4956875/150‑1 (MQ=255)
             aCCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGc                                                                                                                                          >  2:4858073/1‑150 (MQ=255)
                gCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCaa                                                                                                                                       >  2:1945591/1‑150 (MQ=255)
                 cTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCaaa                                                                                                                                      <  1:1151417/150‑1 (MQ=255)
                  tGACGCTGGAAGCCCAGCTGGAGGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAg                                                                                                                                     >  1:3868160/1‑150 (MQ=255)
                  tGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAg                                                                                                                                     <  2:220585/150‑1 (MQ=255)
                   gACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAgg                                                                                                                                    <  1:2153312/150‑1 (MQ=255)
                            aGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCAt                                                                                                                           <  2:218576/150‑1 (MQ=255)
                             gCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATg                                                                                                                          >  2:475386/1‑150 (MQ=255)
                             gCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATg                                                                                                                          <  2:1323174/150‑1 (MQ=255)
                              cccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATgg                                                                                                                         <  2:4562020/150‑1 (MQ=255)
                               ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGt                                                                                                                        <  2:1029908/150‑1 (MQ=255)
                                cAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTg                                                                                                                       >  1:4488769/1‑150 (MQ=255)
                                  gCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCt                                                                                                                     <  2:2525826/150‑1 (MQ=255)
                                  gCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCt                                                                                                                     <  2:3111483/150‑1 (MQ=255)
                                  gCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCt                                                                                                                     >  2:1801977/1‑150 (MQ=255)
                                    tGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGa                                                                                                                   >  2:4001847/1‑150 (MQ=255)
                                     ggAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGAc                                                                                                                  <  2:3058902/150‑1 (MQ=255)
                                     ggAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGAc                                                                                                                  <  1:958426/150‑1 (MQ=255)
                                      gAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACt                                                                                                                 <  1:43302/150‑1 (MQ=255)
                                      gAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACt                                                                                                                 >  2:1066649/1‑150 (MQ=255)
                                              tGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGaa                                                                                                         <  1:2878607/150‑1 (MQ=255)
                                               gcgcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGaaa                                                                                                        <  2:3172364/150‑1 (MQ=255)
                                                 gcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaaga                                                                                                      <  1:4434834/150‑1 (MQ=255)
                                                 gcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaaga                                                                                                      <  1:3374837/150‑1 (MQ=255)
                                                 gcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaaga                                                                                                      >  2:3664125/1‑150 (MQ=255)
                                                 gcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGGTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaaga                                                                                                      <  1:3597186/150‑1 (MQ=255)
                                                   tgatgaTGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaag                                                                                                    >  1:1485337/1‑150 (MQ=255)
                                                    gatgatGTCTACCAACAACGTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGc                                                                                                   >  2:1525447/1‑150 (MQ=255)
                                                    gatgatGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGc                                                                                                   <  1:2162058/150‑1 (MQ=255)
                                                     atgatgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCa                                                                                                  <  1:1366936/150‑1 (MQ=255)
                                                             tACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCt                                                                                          >  2:1371143/1‑150 (MQ=255)
                                                               ccAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGt                                                                                        >  2:444846/1‑150 (MQ=255)
                                                                  acaacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATc                                                                                     <  1:2996049/150‑1 (MQ=255)
                                                                    aacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTACCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACTCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGc                                                                                   >  1:3982227/1‑150 (MQ=255)
                                                                     acaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGct                                                                                  <  2:548366/150‑1 (MQ=255)
                                                                       aTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGctct                                                                                >  2:1638394/1‑150 (MQ=255)
                                                                         ccTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctgg                                                                              <  1:1235613/150‑1 (MQ=255)
                                                                         ccTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctgg                                                                              <  1:845958/150‑1 (MQ=255)
                                                                          cTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGACGCAGAACGTCTGTATCGCtctggt                                                                             >  1:3229864/1‑150 (MQ=255)
                                                                          cTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggt                                                                             >  1:623351/1‑150 (MQ=255)
                                                                            gTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtct                                                                           <  1:1271151/150‑1 (MQ=255)
                                                                                ccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGCCCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCt                                                                       <  2:1840357/150‑1 (MQ=255)
                                                                                ccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCt                                                                       <  1:3044621/150‑1 (MQ=255)
                                                                                 cggcgaacggcgaaCCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCtt                                                                      <  2:4398728/150‑1 (MQ=255)
                                                                                  ggcgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtc                                                                     <  1:1989695/150‑1 (MQ=255)
                                                                                  ggcgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtc                                                                     <  1:1989802/150‑1 (MQ=255)
                                                                                       acggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGc                                                                <  1:4975474/150‑1 (MQ=255)
                                                                                       acggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGc                                                                <  1:2107178/150‑1 (MQ=255)
                                                                                        cggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCa                                                               >  1:4455646/1‑150 (MQ=255)
                                                                                          gcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATg                                                             >  1:2488983/1‑150 (MQ=255)
                                                                                          gcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGAACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATg                                                             <  2:2710332/150‑1 (MQ=255)
                                                                                             aacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgc                                                          <  2:2896356/150‑1 (MQ=255)
                                                                                             aacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgc                                                          <  2:2582581/150‑1 (MQ=255)
                                                                                             aacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgc                                                          <  2:5160375/150‑1 (MQ=255)
                                                                                              acCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgcg                                                         <  1:585220/150‑1 (MQ=255)
                                                                                               ccAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgcgc                                                        >  1:3737649/1‑150 (MQ=255)
                                                                                                cAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGcgcgcg                                                       >  1:2095306/1‑150 (MQ=255)
                                                                                                 aaTCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGt                                                      <  1:2951553/150‑1 (MQ=255)
                                                                                                  atcatcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGtt                                                     >  2:1431266/1‑150 (MQ=255)
                                                                                                    catcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaa                                                   >  1:3448793/1‑150 (MQ=255)
                                                                                                    catcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaa                                                   >  2:3331442/1‑150 (MQ=255)
                                                                                                     atcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaaa                                                  <  1:774672/150‑1 (MQ=255)
                                                                                                     atcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaaa                                                  >  1:5057949/1‑150 (MQ=255)
                                                                                                      tcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAg                                                 <  2:1912478/150‑1 (MQ=255)
                                                                                                      tcGGTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAg                                                 <  1:2890840/150‑1 (MQ=255)
                                                                                                          tCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCg                                             >  2:2439458/1‑150 (MQ=255)
                                                                                                            cGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTa                                           >  2:672028/1‑150 (MQ=255)
                                                                                                             gTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTAt                                          <  2:2276509/150‑1 (MQ=255)
                                                                                                             gTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTAt                                          <  2:2814232/150‑1 (MQ=255)
                                                                                                              tctcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATc                                         <  1:2756169/150‑1 (MQ=255)
                                                                                                              tctcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATc                                         <  2:4024180/150‑1 (MQ=255)
                                                                                                               ctcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCa                                        <  2:1665079/150‑1 (MQ=255)
                                                                                                               ctcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCa                                        >  2:5097078/1‑150 (MQ=255)
                                                                                                                 cAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAcc                                      <  2:3382528/150‑1 (MQ=255)
                                                                                                                 cAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAcc                                      >  2:4726751/1‑150 (MQ=255)
                                                                                                                      cGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTa                                 >  1:1878866/1‑150 (MQ=255)
                                                                                                                      cGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTa                                 >  1:195167/1‑150 (MQ=255)
                                                                                                                      cGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTa                                 >  1:2027536/1‑150 (MQ=255)
                                                                                                                       gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                                <  1:100499/150‑1 (MQ=255)
                                                                                                                        ttGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATg                               <  1:1047573/150‑1 (MQ=255)
                                                                                                                         tGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGa                              <  1:3456570/150‑1 (MQ=255)
                                                                                                                             cTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGaaaaa                          <  2:3251676/150‑1 (MQ=255)
                                                                                                                                   cTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGACAAAGATGCt                    >  2:4145117/1‑150 (MQ=255)
                                                                                                                                    tGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTCTGAACAAGATGCTa                   >  1:4988553/1‑150 (MQ=255)
                                                                                                                                       actacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACg                <  1:4818235/150‑1 (MQ=255)
                                                                                                                                        ctacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACgg               >  2:381202/1‑150 (MQ=255)
                                                                                                                                           cATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGa            <  1:2227411/150‑1 (MQ=255)
                                                                                                                                              gACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAAtt         >  1:2367852/1‑150 (MQ=255)
                                                                                                                                              gACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAAtt         >  2:4239638/1‑150 (MQ=255)
                                                                                                                                                 ccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt      <  1:2657508/150‑1 (MQ=255)
                                                                                                                                                 ccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt      <  2:4822620/150‑1 (MQ=255)
                                                                                                                                                  cGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagta     >  1:4659833/1‑150 (MQ=255)
                                                                                                                                                   gTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagtag    >  2:4247315/1‑150 (MQ=255)
                                                                                                                                                     gACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTAGCg  <  2:2961236/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTAGCG  >  CP009273/4176674‑4176972

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: