Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,172,507 C→A 100% I445I (ATC→ATA rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,172,5070CA97.4% 851.2 / 9.0 269I445I (ATC→ATArpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (144/118);  total (147/122)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 271 total aligned reads displayed.

TGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAA  >  CP009273/4172362‑4172653
                                                                                                                                                 |                                                                                                                                                  
 gTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGacca                                                                                                                                               >  2:3388517/1‑150 (MQ=255)
  tCTGCGGTTGGTCGGATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTAGGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCTATGATATAGaccac                                                                                                                                              >  2:2983497/1‑150 (MQ=255)
   cTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGaccacc                                                                                                                                             >  2:91920/1‑150 (MQ=255)
   cTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGaccacc                                                                                                                                             >  2:4284815/1‑150 (MQ=255)
      cGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCg                                                                                                                                          >  1:386445/1‑150 (MQ=255)
           ggTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAAc                                                                                                                                     >  1:977078/1‑150 (MQ=255)
            gTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAAcc                                                                                                                                    <  2:5067270/150‑1 (MQ=255)
               gTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAACAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACcgtc                                                                                                                                 <  2:4121294/150‑1 (MQ=255)
                 aTGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcgt                                                                                                                               >  2:1730950/1‑150 (MQ=255)
                  tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTa                                                                                                                              <  1:508442/150‑1 (MQ=255)
                   gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAt                                                                                                                             <  2:3509507/150‑1 (MQ=255)
                    aaGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtc                                                                                                                            >  1:798487/1‑150 (MQ=255)
                    aaGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtc                                                                                                                            <  1:4719973/150‑1 (MQ=255)
                     aGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCACAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtcc                                                                                                                           >  2:4789910/1‑150 (MQ=255)
                      gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccg                                                                                                                          >  2:4620756/1‑150 (MQ=255)
                      gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccg                                                                                                                          >  2:2740993/1‑150 (MQ=255)
                      gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccg                                                                                                                          >  1:2315217/1‑150 (MQ=255)
                      gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccg                                                                                                                          <  2:5012847/150‑1 (MQ=255)
                         cAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                       >  2:773850/1‑150 (MQ=255)
                            ccGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttccgt                                                                                                                    <  2:664163/150‑1 (MQ=255)
                               ttCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGAGATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTgg                                                                                                                 >  2:4717256/1‑150 (MQ=255)
                                tctcTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGc                                                                                                                >  1:1101346/1‑150 (MQ=255)
                                 ctctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                                                                                                               <  1:623221/150‑1 (MQ=255)
                                 ctctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                                                                                                               >  2:215441/1‑150 (MQ=255)
                                  tctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGa                                                                                                              <  2:2430765/150‑1 (MQ=255)
                                   ctgctgCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                             <  2:2905368/150‑1 (MQ=255)
                                   ctgctgCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                             >  1:3488546/1‑150 (MQ=255)
                                    tgctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaaa                                                                                                            <  1:1264320/150‑1 (MQ=255)
                                    tgctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaaa                                                                                                            >  2:779285/1‑150 (MQ=255)
                                     gctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                           <  1:3120286/150‑1 (MQ=255)
                                      ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                                                                                                          >  1:306567/1‑150 (MQ=255)
                                       tgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAAAGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATgg                                                                                                         <  2:1445392/150‑1 (MQ=255)
                                         cgcgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCg                                                                                                       >  2:2124587/1‑150 (MQ=255)
                                           cgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                                     >  1:1823039/1‑150 (MQ=255)
                                           cgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                                     <  1:3306920/150‑1 (MQ=255)
                                            gaagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaa                                                                                                    >  2:3583511/1‑150 (MQ=255)
                                             aagatATCGACGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa                                                                                                   <  2:4332852/150‑1 (MQ=255)
                                             aagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa                                                                                                   >  1:5156367/1‑150 (MQ=255)
                                              agaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaaa                                                                                                  >  1:1869768/1‑150 (MQ=255)
                                              agaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACCGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGACGTATCCGTTCCGTTGGCGAAATGGCGGaaaa                                                                                                  >  2:2665577/1‑150 (MQ=255)
                                                   tCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                             >  2:3576179/1‑150 (MQ=255)
                                                    cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                            <  2:1042668/150‑1 (MQ=255)
                                                     gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTc                                                                                           <  2:2998983/150‑1 (MQ=255)
                                                       aGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTcca                                                                                         <  2:1575459/150‑2 (MQ=255)
                                                         gTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCcgcg                                                                                       <  1:552521/150‑1 (MQ=255)
                                                         gTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCcgcg                                                                                       <  2:1228726/150‑1 (MQ=255)
                                                          ttCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGt                                                                                      >  1:3397014/1‑150 (MQ=255)
                                                                 aTCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTg                                                                               >  2:2437598/1‑150 (MQ=255)
                                                                 aTCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTg                                                                               >  2:1928305/1‑150 (MQ=255)
                                                                   ccTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGTTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGt                                                                             >  1:2659189/1‑150 (MQ=255)
                                                                   ccTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGt                                                                             >  1:4522768/1‑150 (MQ=255)
                                                                      gAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACg                                                                          >  2:3619938/1‑150 (MQ=255)
                                                                          aaaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTa                                                                      >  1:2621231/1‑150 (MQ=255)
                                                                             gacgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTagag                                                                   >  2:3798803/1‑150 (MQ=255)
                                                                                 acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTg                                                               <  2:3315165/150‑1 (MQ=255)
                                                                                 acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTg                                                               <  1:1706733/150‑1 (MQ=255)
                                                                                  catcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGc                                                              <  2:4444000/150‑1 (MQ=255)
                                                                                       ttGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGa                                                         >  2:717514/1‑150 (MQ=255)
                                                                                        tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGaa                                                        <  1:2543448/150‑1 (MQ=255)
                                                                                        tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGaa                                                        >  1:1478272/1‑150 (MQ=255)
                                                                                          aTGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAag                                                      <  2:4515072/150‑1 (MQ=255)
                                                                                           tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTCCGTGTAGAGCGTGCGGTGAAaga                                                     >  2:2562523/1‑150 (MQ=255)
                                                                                            gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAagag                                                    >  2:1617824/1‑150 (MQ=255)
                                                                                             ttATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGc                                                   <  1:382161/150‑1 (MQ=255)
                                                                                                tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                                >  2:639602/1‑150 (MQ=255)
                                                                                                tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                                >  1:2800179/1‑150 (MQ=255)
                                                                                                 gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtct                                               >  2:2199030/1‑150 (MQ=255)
                                                                                                   aaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgt                                             <  1:4893556/150‑1 (MQ=255)
                                                                                                   aaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgt                                             <  1:1053198/150‑1 (MQ=255)
                                                                                                    aaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgtc                                            <  1:2380345/150‑1 (MQ=255)
                                                                                                     aaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgtct                                           >  1:3814162/1‑150 (MQ=255)
                                                                                                     aaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgtct                                           <  1:4110394/150‑1 (MQ=255)
                                                                                                      aGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                          <  1:3258472/150‑1 (MQ=255)
                                                                                                       gCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTctct                                         <  1:3576184/150‑1 (MQ=255)
                                                                                                        cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTg                                        <  1:4489756/150‑1 (MQ=255)
                                                                                                             cGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGa                                   <  2:5120914/150‑1 (MQ=255)
                                                                                                              gATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGAt                                  >  2:3282941/1‑150 (MQ=255)
                                                                                                                tatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCt                                <  2:2315217/150‑1 (MQ=255)
                                                                                                                tatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCt                                >  2:374754/1‑150 (MQ=255)
                                                                                                                tatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCt                                <  2:977078/150‑1 (MQ=255)
                                                                                                                  tCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTgg                              <  1:1450685/150‑1 (MQ=255)
                                                                                                                    cGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                            >  2:2024170/1‑150 (MQ=255)
                                                                                                                    cGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACAAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                            >  2:165023/1‑150 (MQ=255)
                                                                                                                      tAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAc                          <  2:3732103/150‑1 (MQ=255)
                                                                                                                      tAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAc                          >  1:4284535/1‑150 (MQ=255)
                                                                                                                      tAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAc                          >  1:363261/1‑150 (MQ=255)
                                                                                                                        aCGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAccc                        <  2:1254796/150‑1 (MQ=255)
                                                                                                                         cGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCt                       >  1:3558020/1‑150 (MQ=255)
                                                                                                                           gTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                     <  1:3281197/150‑1 (MQ=255)
                                                                                                                            tAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGAt                    <  1:2207653/150‑1 (MQ=255)
                                                                                                                            tAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGAt                    >  1:1654226/1‑150 (MQ=255)
                                                                                                                                  cGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcaca              <  1:2437598/150‑1 (MQ=255)
                                                                                                                                   gAAGTCGATGATATATACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCTGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAg             >  2:1640224/1‑150 (MQ=255)
                                                                                                                                   gAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAg             >  2:1979663/1‑150 (MQ=255)
                                                                                                                                   gAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAg             >  1:891680/1‑150 (MQ=255)
                                                                                                                                   gAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAg             <  1:1752254/150‑1 (MQ=255)
                                                                                                                                    aaGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAgg            <  2:3151882/150‑1 (MQ=255)
                                                                                                                                     aGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGa           >  1:2144836/1‑150 (MQ=255)
                                                                                                                                        cGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGatat        <  2:546606/150‑1 (MQ=255)
                                                                                                                                              tataGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCaa  >  1:1397605/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAA  >  CP009273/4172362‑4172653

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: