Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,176,823 G→T 94.3% 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,8230GT94.3% 1238.3 / 57.2 386D516Y (GAC→TAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (10/12);  new base T (152/212);  total (162/224)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 386 total aligned reads displayed.

GCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTA  >  CP009273/4176676‑4176969
                                                                                                                                                   |                                                                                                                                                  
gCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTAc                                                                                                                                                  >  2:4414130/1‑150 (MQ=255)
  tGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACtg                                                                                                                                                >  1:4357263/1‑150 (MQ=255)
     tCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTgtgt                                                                                                                                             <  2:2352532/150‑1 (MQ=255)
      cACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGtt                                                                                                                                            >  2:561967/1‑150 (MQ=255)
       aCGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTa                                                                                                                                           <  1:3271803/150‑1 (MQ=255)
       aCGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTa                                                                                                                                           <  1:820422/150‑1 (MQ=255)
          tACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACg                                                                                                                                        <  1:3625931/150‑1 (MQ=255)
            ccGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGcc                                                                                                                                      <  2:3359888/150‑1 (MQ=255)
             cGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCa                                                                                                                                     <  1:1033042/150‑1 (MQ=255)
               cTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCaaa                                                                                                                                   <  2:912031/150‑1 (MQ=255)
                 gACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAgg                                                                                                                                 <  2:664745/150‑1 (MQ=255)
                   cGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCg                                                                                                                               >  2:2391666/1‑150 (MQ=255)
                      tGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAg                                                                                                                            >  1:2698192/1‑150 (MQ=255)
                         aaGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCa                                                                                                                         >  1:4751174/1‑150 (MQ=255)
                         aaGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACAAGACCCGTTACTGAGTTAACGCCAAAGGCGAAGGCa                                                                                                                         >  1:717569/1‑150 (MQ=255)
                          aGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCAt                                                                                                                        >  2:1644547/1‑150 (MQ=255)
                           gCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATg                                                                                                                       <  1:2407225/150‑1 (MQ=255)
                                  tGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGa                                                                                                                >  1:497531/1‑150 (MQ=255)
                                    gAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACt                                                                                                              <  2:3122916/150‑1 (MQ=255)
                                     aaGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTg                                                                                                             <  2:4206255/150‑1 (MQ=255)
                                       gcgcGTGCGCTGATGATGTCTACCAACAACCTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGc                                                                                                           <  1:2891208/150‑1 (MQ=255)
                                       gcgcGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGc                                                                                                           <  2:4552073/150‑1 (MQ=255)
                                        cgcgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGcc                                                                                                          <  1:4758923/150‑1 (MQ=255)
                                        cgcgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGcc                                                                                                          >  1:1258423/1‑150 (MQ=255)
                                          cgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATTGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCg                                                                                                        <  2:3421755/150‑1 (MQ=255)
                                          cgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCg                                                                                                        >  1:1084843/1‑150 (MQ=255)
                                            tGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGaa                                                                                                      <  2:13473/150‑1 (MQ=255)
                                             gcgcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGaaa                                                                                                     <  1:1119316/150‑1 (MQ=255)
                                              cgcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaag                                                                                                    <  1:1380378/150‑1 (MQ=255)
                                                cTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaa                                                                                                  <  1:4069985/150‑1 (MQ=255)
                                                cTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaa                                                                                                  >  1:260839/1‑150 (MQ=255)
                                                 tgatgaTGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaag                                                                                                 <  2:147128/150‑1 (MQ=255)
                                                   atgatgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCa                                                                                               >  2:4071132/1‑150 (MQ=255)
                                                   atgatgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCa                                                                                               <  2:2492401/150‑1 (MQ=255)
                                                    tgatgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAg                                                                                              >  1:455741/1‑150 (MQ=255)
                                                       tgTCTACCAACAACATCCTGTCCCGGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAAc                                                                                           >  1:1737704/1‑150 (MQ=255)
                                                       tgTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAAc                                                                                           <  2:3651975/150‑1 (MQ=255)
                                                          cTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTc                                                                                        >  1:2104529/1‑150 (MQ=255)
                                                           tACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCt                                                                                       >  1:1482331/1‑150 (MQ=255)
                                                            aCCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTg                                                                                      <  1:2015708/150‑1 (MQ=255)
                                                             ccAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGt                                                                                     <  2:3601275/150‑1 (MQ=255)
                                                              caacaaCATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTa                                                                                    >  2:3172736/1‑150 (MQ=255)
                                                               aacaacATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGACGCAGAACGTCTGTAt                                                                                   >  1:3989854/1‑150 (MQ=255)
                                                                     aTCCTGTCCCCGGCGACCGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAACGGCGAAGGCATGGTGCTGACGGGCCCGAAAGAAGCAGACCGTCTGTATCGctct                                                                             >  1:653784/1‑150 (MQ=255)
                                                                     aTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGctct                                                                             <  1:4681396/150‑1 (MQ=255)
                                                                       ccTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctgg                                                                           <  1:4099141/150‑1 (MQ=255)
                                                                       ccTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctgg                                                                           >  2:85373/1‑150 (MQ=255)
                                                                        cTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggt                                                                          >  1:88172/1‑150 (MQ=255)
                                                                        cTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggt                                                                          <  2:4181043/150‑1 (MQ=255)
                                                                         tGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTCTCGCtctggtc                                                                         >  1:2526778/1‑150 (MQ=255)
                                                                          gTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtct                                                                        <  2:2688929/150‑1 (MQ=255)
                                                                          gTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtct                                                                        >  2:550019/1‑150 (MQ=255)
                                                                           tCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtctg                                                                       <  2:388198/150‑1 (MQ=255)
                                                                            ccccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtctgg                                                                      <  2:1530777/150‑1 (MQ=255)
                                                                                 gcgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTTTGTATCGCTCTGGTCTGGCTnct                                                                 <  2:2831636/150‑1 (MQ=255)
                                                                                 gcgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtct                                                                 <  1:2152810/150‑1 (MQ=255)
                                                                                  cgaacggcgaacCAAACATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCAGTTACTGTGTTAACTCCAAAGGCGAAGGCATGGTGCTGACTGGCCCTAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtctc                                                                <  1:355361/150‑1 (MQ=255)
                                                                                      cggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCa                                                            >  2:3697191/1‑150 (MQ=255)
                                                                                         cgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgc                                                         >  1:182153/1‑150 (MQ=255)
                                                                                          gaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcg                                                        >  2:3000815/1‑150 (MQ=255)
                                                                                            acCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgcg                                                      >  2:2111375/1‑150 (MQ=255)
                                                                                            acCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgcg                                                      <  2:4581133/150‑1 (MQ=255)
                                                                                                  catcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaa                                                <  1:3365489/150‑1 (MQ=255)
                                                                                                   atcGGTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTaaa                                               <  2:1904759/150‑1 (MQ=255)
                                                                                                    tcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAg                                              <  1:1963866/150‑1 (MQ=255)
                                                                                                     cGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGt                                             <  2:4120843/150‑1 (MQ=255)
                                                                                                     cGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGt                                             >  2:3280046/1‑150 (MQ=255)
                                                                                                         ccGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGCAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGt                                         <  1:1850648/150‑1 (MQ=255)
                                                                                                         ccGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGt                                         <  2:2710244/150‑1 (MQ=255)
                                                                                                           gTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTAt                                       <  1:3683934/150‑1 (MQ=255)
                                                                                                                 ggACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGa                                 <  2:1249440/150‑1 (MQ=255)
                                                                                                                  gACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAg                                <  2:4334010/150‑1 (MQ=255)
                                                                                                                   aCGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGt                               <  2:4755707/150‑1 (MQ=255)
                                                                                                                     gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                             >  1:3184573/1‑150 (MQ=255)
                                                                                                                      ttGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATg                            <  1:4256541/150‑1 (MQ=255)
                                                                                                                       tGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGa                           <  2:981533/150‑1 (MQ=255)
                                                                                                                        gTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGaa                          <  1:3486662/150‑1 (MQ=255)
                                                                                                                        gTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGaa                          >  2:2638570/1‑150 (MQ=255)
                                                                                                                             gggTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGa                     >  1:4171813/1‑150 (MQ=255)
                                                                                                                             gggTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGa                     <  2:4706965/150‑1 (MQ=255)
                                                                                                                              ggTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGAt                    <  2:2204618/150‑1 (MQ=255)
                                                                                                                                   gTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTaa               >  2:756324/1‑150 (MQ=255)
                                                                                                                                   gTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTaa               <  1:2321084/150‑1 (MQ=255)
                                                                                                                                   gTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTaa               >  2:3553008/1‑150 (MQ=255)
                                                                                                                                    tactacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAAc              >  2:443936/1‑150 (MQ=255)
                                                                                                                                      ctacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACgg            <  2:4212692/150‑1 (MQ=255)
                                                                                                                                      ctacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACgg            >  1:473618/1‑150 (MQ=255)
                                                                                                                                       tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt           <  1:2762638/150‑1 (MQ=255)
                                                                                                                                       tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt           <  2:634858/150‑1 (MQ=255)
                                                                                                                                       tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt           <  2:1807581/150‑1 (MQ=255)
                                                                                                                                       tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt           <  2:3156966/150‑1 (MQ=255)
                                                                                                                                        acATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTg          >  1:3125691/1‑150 (MQ=255)
                                                                                                                                        acATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAGATGCTAACGGTGa         <  1:971388/150‑1 (MQ=255)
                                                                                                                                          aTGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCTTGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGaa        <  1:2989786/150‑1 (MQ=255)
                                                                                                                                              cccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtag    >  1:4446271/1‑150 (MQ=255)
                                                                                                                                              cccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtag    >  1:3786273/1‑150 (MQ=255)
                                                                                                                                              cccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtag    >  2:3518622/1‑150 (MQ=255)
                                                                                                                                               ccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt   <  2:2656772/150‑1 (MQ=255)
                                                                                                                                               ccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt   >  2:3401976/1‑150 (MQ=255)
                                                                                                                                                cGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagta  <  1:340447/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTA  >  CP009273/4176676‑4176969

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: