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,5070CA98.1% 808.5 / 5.6 257I445I (ATC→ATArpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (119/133);  total (121/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 258 total aligned reads displayed.

TGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCA  >  CP009273/4172362‑4172652
                                                                                                                                                 |                                                                                                                                                 
tGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGac                                                                                                                                                <  1:2072871/149‑1 (MQ=255)
tGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGac                                                                                                                                                <  1:1442769/149‑1 (MQ=255)
 gTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGacc                                                                                                                                               <  2:2884167/149‑1 (MQ=255)
  tCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGacca                                                                                                                                              >  2:2296717/1‑149 (MQ=255)
      cGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTc                                                                                                                                          <  1:250973/149‑1 (MQ=255)
         ttGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGc                                                                                                                                       <  1:2518759/149‑1 (MQ=255)
           ggTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTACAGGCGAAGTCGATGATATAGACCACCTCGGCaa                                                                                                                                     >  2:1081727/1‑149 (MQ=255)
           ggTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCaa                                                                                                                                     >  2:2507213/1‑149 (MQ=255)
           ggTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCaa                                                                                                                                     <  1:2366160/149‑1 (MQ=255)
           ggTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTAGGCaa                                                                                                                                     >  2:2600146/1‑149 (MQ=255)
            gTCGTATGAAGTTAAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAAc                                                                                                                                    <  1:1294662/149‑1 (MQ=255)
                tATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtc                                                                                                                                >  2:1532818/1‑149 (MQ=255)
                tATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtc                                                                                                                                <  2:4034434/149‑1 (MQ=255)
                   gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTa                                                                                                                             >  2:1876159/1‑149 (MQ=255)
                    aaGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAt                                                                                                                            >  1:853885/1‑149 (MQ=255)
                         cAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgtt                                                                                                                       <  1:1728523/149‑1 (MQ=255)
                          aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttc                                                                                                                      <  1:351573/149‑1 (MQ=255)
                           aCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttcc                                                                                                                     >  1:3014832/1‑149 (MQ=255)
                             cGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgttccgt                                                                                                                   <  2:1093926/149‑1 (MQ=255)
                               ttCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                                 <  1:3441888/149‑1 (MQ=255)
                               ttCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                                 <  1:2772753/149‑1 (MQ=255)
                                tctcTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTgg                                                                                                                <  1:1247650/149‑1 (MQ=255)
                                tctcTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTgg                                                                                                                >  1:676718/1‑149 (MQ=255)
                                    tgctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                            <  1:434274/149‑1 (MQ=255)
                                      ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                          <  2:224584/149‑1 (MQ=255)
                                      ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                          <  1:3071352/149‑1 (MQ=255)
                                      ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                          <  1:1714683/149‑1 (MQ=255)
                                      ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                          <  1:2269599/149‑1 (MQ=255)
                                      ctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                          <  2:1985008/149‑1 (MQ=255)
                                        gcgcGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATgg                                                                                                        >  2:1849022/1‑149 (MQ=255)
                                         cgcgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGc                                                                                                       <  2:1296362/149‑1 (MQ=255)
                                         cgcgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGc                                                                                                       <  1:1881553/149‑1 (MQ=255)
                                              agaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa                                                                                                  >  1:3533041/1‑149 (MQ=255)
                                              agaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa                                                                                                  >  2:360928/1‑149 (MQ=255)
                                               gaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaaa                                                                                                 <  1:3159782/149‑1 (MQ=255)
                                                aaaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                                <  2:2786342/149‑1 (MQ=255)
                                                aaaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                                >  2:1703856/1‑149 (MQ=255)
                                                 aaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc                                                                                               <  1:1317608/149‑1 (MQ=255)
                                                  aTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCa                                                                                              <  2:4094715/149‑1 (MQ=255)
                                                    cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                            <  2:2774264/149‑1 (MQ=255)
                                                    cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                            >  1:3416086/1‑149 (MQ=255)
                                                    cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                            >  1:3645668/1‑149 (MQ=255)
                                                     gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                           <  1:671384/149‑1 (MQ=255)
                                                     gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                           >  1:1551505/1‑149 (MQ=255)
                                                           tCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGt                                                                                     <  1:483518/149‑1 (MQ=255)
                                                                tATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGcc                                                                                <  2:3014191/149‑1 (MQ=255)
                                                                 aTCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCt                                                                               >  2:2112018/1‑149 (MQ=255)
                                                                   ccTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTgg                                                                             <  1:1029127/149‑1 (MQ=255)
                                                                    cTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGt                                                                            <  2:638706/149‑1 (MQ=255)
                                                                        gCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgt                                                                        >  2:810309/1‑149 (MQ=255)
                                                                          aaaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgtgt                                                                      <  2:4058602/149‑1 (MQ=255)
                                                                           aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTa                                                                     >  1:3721504/1‑149 (MQ=255)
                                                                           aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTa                                                                     >  2:2934411/1‑149 (MQ=255)
                                                                           aaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTa                                                                     >  2:3264369/1‑149 (MQ=255)
                                                                            aGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGTGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTag                                                                    >  2:1688409/1‑149 (MQ=255)
                                                                             gacgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTaga                                                                   >  2:1373408/1‑149 (MQ=255)
                                                                             gacgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTaga                                                                   <  2:1376959/149‑1 (MQ=255)
                                                                             gacgacATCAGTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTaga                                                                   <  2:3821355/149‑1 (MQ=255)
                                                                              acgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGTAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTagag                                                                  >  2:2359775/1‑149 (MQ=255)
                                                                              acgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTagag                                                                  <  2:2917598/149‑1 (MQ=255)
                                                                                    tcatTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCg                                                            >  2:2710425/1‑149 (MQ=255)
                                                                                     catTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCgg                                                           <  2:1281141/149‑1 (MQ=255)
                                                                                      atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGt                                                          <  2:3193292/149‑1 (MQ=255)
                                                                                       ttGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTg                                                         >  1:1915140/1‑149 (MQ=255)
                                                                                                 gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                               <  1:1867516/149‑1 (MQ=255)
                                                                                                  aaaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtct                                              <  2:3573205/149‑1 (MQ=255)
                                                                                                   aaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctg                                             >  2:1229959/1‑149 (MQ=255)
                                                                                                       gCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                         <  1:2637160/149‑1 (MQ=255)
                                                                                                       gCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                         >  2:1660184/1‑149 (MQ=255)
                                                                                                        cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTctct                                        >  1:1484471/1‑149 (MQ=255)
                                                                                                         tCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTg                                       >  2:3433089/1‑149 (MQ=255)
                                                                                                           aTCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTggg                                     >  2:412482/1‑149 (MQ=255)
                                                                                                           aTCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTggg                                     >  1:4050307/1‑149 (MQ=255)
                                                                                                                tatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATc                                >  2:1466886/1‑149 (MQ=255)
                                                                                                                  tCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTg                              >  2:1243167/1‑149 (MQ=255)
                                                                                                                  tCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCAGGGCGATCTg                              >  2:3001240/1‑149 (MQ=255)
                                                                                                                   ccGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTgg                             >  2:3438243/1‑149 (MQ=255)
                                                                                                                     gTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGAt                           <  1:1101826/149‑1 (MQ=255)
                                                                                                                       aaCGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAc                         >  1:2757243/1‑149 (MQ=255)
                                                                                                                          ggTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCt                      >  2:390025/1‑149 (MQ=255)
                                                                                                                           gTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTg                     <  2:824613/149‑1 (MQ=255)
                                                                                                                            tAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                    >  2:312296/1‑149 (MQ=255)
                                                                                                                            tAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGa                    <  1:196602/149‑1 (MQ=255)
                                                                                                                                ggCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGcc                <  2:1884883/149‑1 (MQ=255)
                                                                                                                                ggCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGcc                <  1:2507213/149‑1 (MQ=255)
                                                                                                                                  cGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcac              <  2:1304520/149‑1 (MQ=255)
                                                                                                                                   gAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCcaca             <  1:53886/149‑1 (MQ=255)
                                                                                                                                     aGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAgg           >  2:530675/1‑149 (MQ=255)
                                                                                                                                      gTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGa          <  2:879436/149‑1 (MQ=255)
                                                                                                                                      gTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGa          >  2:2135863/1‑149 (MQ=255)
                                                                                                                                       tCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGat         <  1:706307/149‑1 (MQ=255)
                                                                                                                                       tCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGat         <  1:2445172/149‑1 (MQ=255)
                                                                                                                                       tCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGat         >  2:1828240/1‑149 (MQ=255)
                                                                                                                                         gatgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGatat       <  2:3583611/149‑1 (MQ=255)
                                                                                                                                         gatgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGatat       >  1:1192391/1‑149 (MQ=255)
                                                                                                                                          atgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATg      <  1:3821283/149‑1 (MQ=255)
                                                                                                                                           tgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGa     <  2:2756747/149‑1 (MQ=255)
                                                                                                                                            gatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGAt    <  2:599297/149‑1 (MQ=255)
                                                                                                                                              tatCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCa  >  1:3939107/1‑149 (MQ=255)
                                                                                                                                              tatCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCa  >  1:938514/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCA  >  CP009273/4172362‑4172652

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: