Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,172,507 C→A 92.6% 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,5070CA92.6% 523.0 / 34.1 190I445I (ATC→ATArpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base A (96/80);  total (103/87)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 194 total aligned reads displayed.

TCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATC  >  CP009273/4172364‑4172651
                                                                                                                                               |                                                                                                                                                
       ttGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCgg                                                                                                                                       >  1:543291/1‑148 (MQ=255)
            cGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAAcc                                                                                                                                  <  1:1560344/148‑1 (MQ=255)
              tATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgt                                                                                                                                <  1:2828342/148‑1 (MQ=255)
              tATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgt                                                                                                                                >  1:1497885/1‑148 (MQ=255)
               aTGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtc                                                                                                                               >  2:2019638/1‑148 (MQ=255)
                 gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcgt                                                                                                                             >  2:3327446/1‑148 (MQ=255)
                    gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtc                                                                                                                          >  1:2747824/1‑148 (MQ=255)
                    gTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtc                                                                                                                          >  1:2741218/1‑148 (MQ=255)
                     ttCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtcc                                                                                                                         >  1:366644/1‑148 (MQ=255)
                         aCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                     >  1:1828330/1‑148 (MQ=255)
                          ccGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttcc                                                                                                                    >  1:1329626/1‑148 (MQ=255)
                           cGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttccg                                                                                                                   <  2:2170048/148‑1 (MQ=255)
                           cGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttccg                                                                                                                   <  2:2239335/148‑1 (MQ=255)
                            gTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttccgt                                                                                                                  <  1:2695099/148‑1 (MQ=255)
                              tctcTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                                >  1:2917594/1‑148 (MQ=255)
                              tctcTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                                >  1:1600521/1‑148 (MQ=255)
                              tctcTGCTGCGCGAAGAAATCGAACGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTAATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                                >  1:2917647/1‑148 (MQ=255)
                                   gctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                           <  1:2595554/148‑1 (MQ=255)
                                   gctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                           >  1:49463/1‑148 (MQ=255)
                                   gctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                           <  2:2323649/148‑1 (MQ=255)
                                    ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaaa                                                                                                          <  1:1967049/148‑1 (MQ=255)
                                      gcgcGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATg                                                                                                        <  2:429315/148‑1 (MQ=255)
                                           aagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                                   <  1:1918368/148‑1 (MQ=255)
                                           aagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                                   >  1:3155751/1‑148 (MQ=255)
                                              aaaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaaa                                                                                                <  2:184648/148‑1 (MQ=255)
                                                aTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc                                                                                              >  2:3512420/1‑148 (MQ=255)
                                                aTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc                                                                                              >  2:697844/1‑148 (MQ=255)
                                                  cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCTTTCCGTTGGCGAAATGGCGGAAAACCAg                                                                                            <  2:2568393/148‑1 (MQ=255)
                                                  cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAg                                                                                            <  1:3199750/148‑1 (MQ=255)
                                                    aaGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                          >  1:2623068/1‑148 (MQ=255)
                                                    aaGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                          <  2:1142422/148‑1 (MQ=255)
                                                          ccGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGt                                                                                    <  1:1226287/148‑1 (MQ=255)
                                                            ggTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTg                                                                                  <  1:789480/148‑1 (MQ=255)
                                                            ggTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTg                                                                                  <  1:789566/148‑1 (MQ=255)
                                                             gTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTgg                                                                                 >  1:2390723/1‑148 (MQ=255)
                                                              tATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGc                                                                                >  2:3090748/1‑148 (MQ=255)
                                                               aTCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGcc                                                                               <  2:1804433/148‑1 (MQ=255)
                                                                 ccTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTg                                                                             >  2:2505589/1‑148 (MQ=255)
                                                                 ccTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTg                                                                             >  2:1145341/1‑148 (MQ=255)
                                                                  cTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTgg                                                                            <  1:2353406/148‑1 (MQ=255)
                                                                     aGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTAc                                                                         >  1:1552245/1‑148 (MQ=255)
                                                                      gCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACg                                                                        >  2:2178286/1‑148 (MQ=255)
                                                                         aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgtgt                                                                     >  1:2336315/1‑148 (MQ=255)
                                                                         aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgtgt                                                                     <  1:3327446/148‑1 (MQ=255)
                                                                         aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgtgt                                                                     <  1:1906484/148‑1 (MQ=255)
                                                                          aGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTa                                                                    <  2:1233697/148‑1 (MQ=255)
                                                                               acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCg                                                               >  2:1142720/1‑148 (MQ=255)
                                                                               acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCg                                                               >  2:2849214/1‑148 (MQ=255)
                                                                                catcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGt                                                              >  2:178690/1‑148 (MQ=255)
                                                                                   catTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCg                                                           >  1:1322916/1‑148 (MQ=255)
                                                                                   catTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCg                                                           <  2:1801826/148‑1 (MQ=255)
                                                                                      tGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTg                                                        >  1:1641263/1‑148 (MQ=255)
                                                                                         tGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGaaa                                                     >  2:2392494/1‑148 (MQ=255)
                                                                                          gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAag                                                    <  1:429763/148‑1 (MQ=255)
                                                                                          gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAag                                                    >  1:1153095/1‑148 (MQ=255)
                                                                                             aTGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGc                                                 >  1:1666833/1‑148 (MQ=255)
                                                                                               gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgt                                               <  2:840946/148‑1 (MQ=255)
                                                                                                 aaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtct                                             <  2:1406821/148‑1 (MQ=255)
                                                                                                   aaGCTTATCTATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTACGTTGGCGAAATGGCGGAAATCCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgt                                           <  2:509948/148‑1 (MQ=255)
                                                                                                    aGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgtc                                          >  2:2301676/1‑148 (MQ=255)
                                                                                                     gCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgtct                                         <  2:3340636/148‑1 (MQ=255)
                                                                                                     gCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgtct                                         <  2:2970804/148‑1 (MQ=255)
                                                                                                      cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                        >  1:3299165/1‑148 (MQ=255)
                                                                                                      cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                        >  1:2126384/1‑148 (MQ=255)
                                                                                                       tCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTctct                                       <  1:697844/148‑1 (MQ=255)
                                                                                                        cATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTg                                      >  2:1463188/1‑148 (MQ=255)
                                                                                                        cATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTg                                      >  2:590111/1‑148 (MQ=255)
                                                                                                         aTCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTgg                                     >  2:2724356/1‑148 (MQ=255)
                                                                                                         aTCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTgg                                     >  1:2222572/1‑148 (MQ=255)
                                                                                                          tCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTggg                                    <  2:2747824/148‑1 (MQ=255)
                                                                                                           cGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGc                                   >  1:3392044/1‑148 (MQ=255)
                                                                                                               atCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATc                               >  2:1051566/1‑148 (MQ=255)
                                                                                                                    tAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                          >  1:2550184/1‑148 (MQ=255)
                                                                                                                    tAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                          >  1:2489939/1‑148 (MQ=255)
                                                                                                                    tAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                          >  1:2487919/1‑148 (MQ=255)
                                                                                                                      aCGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAc                        >  2:1216712/1‑148 (MQ=255)
                                                                                                                      aCGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAc                        >  2:2657096/1‑148 (MQ=255)
                                                                                                                       cGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAcc                       >  2:3135624/1‑148 (MQ=255)
                                                                                                                       cGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAcc                       >  2:946268/1‑148 (MQ=255)
                                                                                                                          tAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTg                    <  2:3246703/148‑1 (MQ=255)
                                                                                                                           aaaGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                   <  2:2741218/148‑1 (MQ=255)
                                                                                                                           aaaGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                   >  1:1004040/1‑148 (MQ=255)
                                                                                                                           aaaGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                   <  1:1989968/148‑1 (MQ=255)
                                                                                                                           aaaGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                   >  1:2985308/1‑148 (MQ=255)
                                                                                                                           aaaGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                   <  2:2192503/148‑1 (MQ=255)
                                                                                                                              ggCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGc                <  1:1039407/148‑1 (MQ=255)
                                                                                                                               gCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGcc               >  2:1109164/1‑148 (MQ=255)
                                                                                                                                   aGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAg           <  2:1555489/148‑1 (MQ=255)
                                                                                                                                    gTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAgg          <  1:946268/148‑1 (MQ=255)
                                                                                                                                    gTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAgg          <  1:2725465/148‑1 (MQ=255)
                                                                                                                                     tCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAggc         >  2:3249152/1‑147 (MQ=255)
                                                                                                                                      cGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGat        <  2:573918/148‑1 (MQ=255)
                                                                                                                                       gatgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGata       <  1:3020008/148‑1 (MQ=255)
                                                                                                                                         tgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATg     >  1:3308781/1‑148 (MQ=255)
                                                                                                                                         tgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATg     >  1:2162295/1‑148 (MQ=255)
                                                                                                                                          gatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGa    >  1:3515399/1‑148 (MQ=255)
                                                                                                                                            tataGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGCTc  >  2:563434/1‑148 (MQ=255)
                                                                                                                                            tataGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATc  >  2:3016127/1‑148 (MQ=255)
                                                                                                                                            tataGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATc  <  1:3508396/148‑1 (MQ=255)
                                                                                                                                            tataGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATc  >  2:1378381/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATC  >  CP009273/4172364‑4172651

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: