Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,176,823 G→T 92.9% 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,8230GT92.9% 1175.5 / 72.3 365D516Y (GAC→TAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (8/18);  new base T (145/194);  total (153/212)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 365 total aligned reads displayed.

GCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTAG  >  CP009273/4176676‑4176970
                                                                                                                                                   |                                                                                                                                                   
gCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTAc                                                                                                                                                   >  2:848191/1‑150 (MQ=255)
   gTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACtgt                                                                                                                                                <  1:169062/150‑1 (MQ=255)
          tACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACg                                                                                                                                         >  1:591659/1‑150 (MQ=255)
             cGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGAACTACATGACCCGTTACTGTGTTAACGCCa                                                                                                                                      <  2:1995526/150‑1 (MQ=255)
                tGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAg                                                                                                                                   >  2:4682830/1‑150 (MQ=255)
                tGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAg                                                                                                                                   >  1:4062903/1‑150 (MQ=255)
                tGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAg                                                                                                                                   <  1:60068/150‑1 (MQ=255)
                 gACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAgg                                                                                                                                  >  2:3118436/1‑150 (MQ=255)
                    gCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGa                                                                                                                               <  2:3336160/150‑1 (MQ=255)
                      tGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAg                                                                                                                             <  2:3233830/150‑1 (MQ=255)
                        gAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGc                                                                                                                           >  2:1600122/1‑150 (MQ=255)
                          aGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCAt                                                                                                                         >  2:18818/1‑150 (MQ=255)
                            cccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATgg                                                                                                                       <  2:38377/150‑1 (MQ=255)
                            cccAGCTGGAAGCGCGTGCACTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGGAGGCATgg                                                                                                                       <  1:2505810/150‑1 (MQ=255)
                             ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGt                                                                                                                      <  2:1233596/150‑1 (MQ=255)
                             ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGt                                                                                                                      <  1:1415778/150‑1 (MQ=255)
                             ccAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGGACTACATGACCCGTGACTGTGTTAACGCCACAGGCGAAGGCATGGt                                                                                                                      >  2:702842/1‑150 (MQ=255)
                              cAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTg                                                                                                                     >  2:4401683/1‑150 (MQ=255)
                                gCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCt                                                                                                                   <  2:3261130/150‑1 (MQ=255)
                                gCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCt                                                                                                                   <  2:2899438/150‑1 (MQ=255)
                                   ggAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGAc                                                                                                                >  1:1033810/1‑150 (MQ=255)
                                     aaGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTg                                                                                                              >  1:2584363/1‑150 (MQ=255)
                                       gcgcGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGc                                                                                                            <  2:4348960/150‑1 (MQ=255)
                                        cgcgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGcc                                                                                                           >  1:2719961/1‑150 (MQ=255)
                                         gcgTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGccc                                                                                                          <  2:210810/150‑1 (MQ=255)
                                             gcgcTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGaaa                                                                                                      >  1:1956/1‑150 (MQ=255)
                                                cTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaa                                                                                                   >  1:3686933/1‑150 (MQ=255)
                                                 tgatgaTGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaag                                                                                                  <  1:2291709/150‑1 (MQ=255)
                                                 tgatgaTGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAaagaag                                                                                                  >  1:748643/1‑150 (MQ=255)
                                                            aCCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTg                                                                                       >  2:4487813/1‑150 (MQ=255)
                                                            aCCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTg                                                                                       >  1:853762/1‑150 (MQ=255)
                                                            aCCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTg                                                                                       <  2:3692286/150‑1 (MQ=255)
                                                             ccAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGt                                                                                      >  2:4731896/1‑150 (MQ=255)
                                                              caacaaCATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTa                                                                                     >  1:1341896/1‑150 (MQ=255)
                                                                acaacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATc                                                                                   <  2:4656236/150‑1 (MQ=255)
                                                                acaacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATc                                                                                   <  2:279517/150‑1 (MQ=255)
                                                                acaacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATc                                                                                   <  1:1440549/150‑1 (MQ=255)
                                                                acaacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATc                                                                                   >  1:1761484/1‑150 (MQ=255)
                                                                  aacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGc                                                                                 >  1:3292360/1‑150 (MQ=255)
                                                                  aacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGc                                                                                 >  1:917081/1‑150 (MQ=255)
                                                                  aacaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGc                                                                                 >  1:3136021/1‑150 (MQ=255)
                                                                   acaTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGct                                                                                >  1:1624179/1‑150 (MQ=255)
                                                                    caTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGctc                                                                               <  1:3947583/150‑1 (MQ=255)
                                                                    caTCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGctc                                                                               <  2:757988/150‑1 (MQ=255)
                                                                      tCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctg                                                                             <  2:2234892/150‑1 (MQ=255)
                                                                      tCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctg                                                                             <  1:957/150‑1 (MQ=255)
                                                                        cTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggt                                                                           >  1:2084638/1‑150 (MQ=255)
                                                                        cTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggt                                                                           >  1:268018/1‑150 (MQ=255)
                                                                         tGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtc                                                                          >  1:2913100/1‑150 (MQ=255)
                                                                          gTCCCCGGCGAACTGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtct                                                                         <  1:4168299/150‑1 (MQ=255)
                                                                            ccccGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCtctggtctgg                                                                       <  2:3799202/150‑1 (MQ=255)
                                                                               cggcgaacggcgaaCCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCtt                                                                    <  2:1500953/150‑1 (MQ=255)
                                                                               cggcgaacggcgaaCCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCtt                                                                    <  2:4446114/150‑1 (MQ=255)
                                                                                ggcgaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTtc                                                                   <  1:3471705/150‑1 (MQ=255)
                                                                                   gaacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTctct                                                                <  2:803394/150‑1 (MQ=255)
                                                                                    aacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAACGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTg                                                               <  1:506480/150‑1 (MQ=255)
                                                                                    aacggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTg                                                               >  1:298918/1‑150 (MQ=255)
                                                                                     acggcgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGc                                                              <  2:1786855/150‑1 (MQ=255)
                                                                                         cgaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgc                                                          >  2:3898239/1‑150 (MQ=255)
                                                                                          gaacCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcg                                                         >  2:4183410/1‑150 (MQ=255)
                                                                                            acCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgcg                                                       <  2:876403/150‑1 (MQ=255)
                                                                                             ccAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATgcgcgc                                                      >  2:1488761/1‑150 (MQ=255)
                                                                                              cAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGcgcgcg                                                     >  2:2724211/1‑150 (MQ=255)
                                                                                              cAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGcgcgcg                                                     >  2:2624525/1‑150 (MQ=255)
                                                                                               aaTCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGt                                                    >  2:4213120/1‑150 (MQ=255)
                                                                                                atcatcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGtt                                                   >  1:4422854/1‑150 (MQ=255)
                                                                                                atcatcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGtt                                                   >  2:4054067/1‑150 (MQ=255)
                                                                                                 tcatcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTa                                                  <  2:1100223/150‑1 (MQ=255)
                                                                                                    tcGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAg                                               <  1:2149667/150‑1 (MQ=255)
                                                                                                     cGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGt                                              <  1:3828292/150‑1 (MQ=255)
                                                                                                      gTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTg                                             >  2:3649518/1‑150 (MQ=255)
                                                                                                      gTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTg                                             <  2:1419265/150‑1 (MQ=255)
                                                                                                         ccGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGt                                          <  1:3976914/150‑1 (MQ=255)
                                                                                                              tcAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAc                                     >  2:4295673/1‑150 (MQ=255)
                                                                                                               cAGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCAcc                                    <  1:3606541/150‑1 (MQ=255)
                                                                                                                aGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCg                                   >  1:1291594/1‑150 (MQ=255)
                                                                                                                aGGACGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCg                                   >  1:3684223/1‑150 (MQ=255)
                                                                                                                    cGTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTa                               <  2:3571540/150‑1 (MQ=255)
                                                                                                                     gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                              >  2:3082840/1‑150 (MQ=255)
                                                                                                                     gTTGTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTAt                              <  2:1304554/150‑1 (MQ=255)
                                                                                                                        gTACTGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGaa                           <  1:4614641/150‑1 (MQ=255)
                                                                                                                            tGGGTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAg                       <  2:2972681/150‑1 (MQ=255)
                                                                                                                             gggTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGa                      >  2:3164021/1‑150 (MQ=255)
                                                                                                                               gTCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATg                    <  2:2595821/150‑1 (MQ=255)
                                                                                                                                tCTGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGc                   <  1:4731896/150‑1 (MQ=255)
                                                                                                                                  tGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTa                 >  1:4505693/1‑150 (MQ=255)
                                                                                                                                  tGTACTACATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTa                 <  2:3292360/150‑1 (MQ=255)
                                                                                                                                     actacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACg              <  1:3015868/150‑1 (MQ=255)
                                                                                                                                     actacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACg              <  2:1017115/150‑1 (MQ=255)
                                                                                                                                       tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt            >  2:2463829/1‑150 (MQ=255)
                                                                                                                                       tacATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGt            >  1:3511264/1‑150 (MQ=255)
                                                                                                                                        acATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTg           <  2:2410409/150‑1 (MQ=255)
                                                                                                                                         cATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGa          >  1:1820431/1‑150 (MQ=255)
                                                                                                                                         cATGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGa          <  1:1137787/150‑1 (MQ=255)
                                                                                                                                           tGACCCGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAAt        <  2:2524202/150‑1 (MQ=255)
                                                                                                                                              cccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtag     >  2:2702098/1‑150 (MQ=255)
                                                                                                                                              cccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtag     <  2:591885/150‑1 (MQ=255)
                                                                                                                                               ccGTTACTGTGTTAACGCCAATGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt    <  1:3497081/150‑1 (MQ=255)
                                                                                                                                               ccGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagt    <  1:1107894/150‑1 (MQ=255)
                                                                                                                                                cGTTACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATtagta   >  2:371928/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCGTTCCGTCTCAGGACGTTGTACTGGGTCTGTACTACATGACCCGTGACTGTGTTAACGCCAAAGGCGAAGGCATGGTGCTGACTGGCCCGAAAGAAGCAGAACGTCTGTATCGCTCTGGTCTGGCTTCTCTGCATGCGCGCGTTAAAGTGCGTATCACCGAGTATGAAAAAGATGCTAACGGTGAATTAGTAG  >  CP009273/4176676‑4176970

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: