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,5070CA96.7% 649.5 / 11.5 213I445I (ATC→ATA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (117/89);  total (120/93)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 213 total aligned reads displayed.

TGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATG  >  CP009273/4172362‑4172648
                                                                                                                                                 |                                                                                                                                             
tGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGa                                                                                                                                             <  1:3228551/148‑1 (MQ=255)
  tCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTTAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGacc                                                                                                                                           <  2:1645714/148‑1 (MQ=255)
  tCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGacc                                                                                                                                           <  2:2262103/148‑1 (MQ=255)
     gCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGaccacc                                                                                                                                        <  1:1538528/148‑1 (MQ=255)
      cGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCt                                                                                                                                       <  2:2652715/148‑1 (MQ=255)
             tCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAAc                                                                                                                                <  2:2476363/148‑1 (MQ=255)
                tATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgt                                                                                                                             <  1:922222/148‑1 (MQ=255)
                  tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcg                                                                                                                           >  1:1183736/1‑148 (MQ=255)
                  tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcg                                                                                                                           >  2:1181696/1‑148 (MQ=255)
                  tGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcg                                                                                                                           >  2:1898729/1‑148 (MQ=255)
                   gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcgt                                                                                                                          >  1:1168434/1‑148 (MQ=255)
                   gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcgt                                                                                                                          <  1:949961/148‑1 (MQ=255)
                   gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcgt                                                                                                                          >  2:2641131/1‑148 (MQ=255)
                   gAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACcgtcgt                                                                                                                          >  2:3105187/1‑148 (MQ=255)
                    aaGTTCAATCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTa                                                                                                                         >  1:1065482/1‑148 (MQ=255)
                     aGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAt                                                                                                                        <  1:2996560/148‑1 (MQ=255)
                       ttCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtcc                                                                                                                      <  2:1100325/148‑1 (MQ=255)
                        tCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccg                                                                                                                     >  1:2175245/1‑148 (MQ=255)
                          aaCCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTAtccgtt                                                                                                                   >  2:3342834/1‑148 (MQ=255)
                                tctcTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACAACCTCGGCAACCGTCGTATCCGTTCCGTTg                                                                                                             >  2:1976661/1‑148 (MQ=255)
                                 ctctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTgg                                                                                                            >  2:1588340/1‑148 (MQ=255)
                                 ctctGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTgg                                                                                                            <  2:1269585/148‑1 (MQ=255)
                                   ctgctgCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                                                                                                          >  2:455622/1‑148 (MQ=255)
                                   ctgctgCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCg                                                                                                          >  1:132852/1‑148 (MQ=255)
                                     gctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                        >  1:725354/1‑148 (MQ=255)
                                     gctgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGaa                                                                                                        <  1:2562756/148‑1 (MQ=255)
                                       tgcGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAAt                                                                                                      <  1:2645033/148‑1 (MQ=255)
                                         cgcgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATgg                                                                                                    <  2:2737779/148‑1 (MQ=255)
                                         cgcgAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATgg                                                                                                    >  1:1785900/1‑148 (MQ=255)
                                             aagaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGa                                                                                                >  2:2987567/1‑148 (MQ=255)
                                               gaaATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGaaa                                                                                              >  1:159194/1‑148 (MQ=255)
                                                 aaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                            <  2:1275885/148‑1 (MQ=255)
                                                 aaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                            <  2:1276999/148‑1 (MQ=255)
                                                 aaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                            <  2:870201/148‑1 (MQ=255)
                                                 aaTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAc                                                                                            >  2:772560/1‑148 (MQ=255)
                                                  aTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc                                                                                           <  1:3442430/148‑1 (MQ=255)
                                                  aTCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAAcc                                                                                           >  2:2551207/1‑148 (MQ=255)
                                                    cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAg                                                                                         >  2:1901666/1‑148 (MQ=255)
                                                    cGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAg                                                                                         >  2:3298921/1‑148 (MQ=255)
                                                     gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                        <  1:1769788/148‑1 (MQ=255)
                                                     gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                        <  1:1663357/148‑1 (MQ=255)
                                                     gAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGt                                                                                        <  1:3111534/148‑1 (MQ=255)
                                                      aaGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                       >  1:2000574/1‑148 (MQ=255)
                                                      aaGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGtt                                                                                       >  1:1024068/1‑148 (MQ=255)
                                                        ggTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTcc                                                                                     >  2:774031/1‑148 (MQ=255)
                                                        ggTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGCCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTcc                                                                                     >  2:2656549/1‑148 (MQ=255)
                                                              ggTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTg                                                                               >  2:2994858/1‑148 (MQ=255)
                                                                 aTCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGcc                                                                            >  1:2151451/1‑148 (MQ=255)
                                                                 aTCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGcc                                                                            <  1:1465608/148‑1 (MQ=255)
                                                                  tCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCt                                                                           <  2:2998450/148‑1 (MQ=255)
                                                                  tCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCt                                                                           <  1:3347464/148‑1 (MQ=255)
                                                                      gAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTa                                                                       >  2:1682122/1‑148 (MQ=255)
                                                                       aGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTAc                                                                      >  2:2389636/1‑148 (MQ=255)
                                                                       aGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTAc                                                                      >  2:2603460/1‑148 (MQ=255)
                                                                        gCAAAGACGACATCATTTATGTTATGAAAAAGCTCATCTATATCCGTAACGGTAAATGCGAAGTCGATGATATAGACCACCTCGGCAAACGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACg                                                                     <  1:2491288/148‑1 (MQ=255)
                                                                         cAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgt                                                                    <  1:3190449/148‑1 (MQ=255)
                                                                         cAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgt                                                                    >  2:712925/1‑148 (MQ=255)
                                                                          aaaGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACgtg                                                                   <  2:2177191/148‑1 (MQ=255)
                                                                            aGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTa                                                                 >  2:1643674/1‑148 (MQ=255)
                                                                              acgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTaga                                                               >  1:675690/1‑148 (MQ=255)
                                                                               cgacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTagag                                                              <  1:1834123/148‑1 (MQ=255)
                                                                                gacATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGc                                                             <  1:348132/148‑1 (MQ=255)
                                                                                 acATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCg                                                            >  2:523346/1‑148 (MQ=255)
                                                                                      atTGATGTTATGAACAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCgg                                                       <  1:2684168/148‑1 (MQ=255)
                                                                                      atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAGAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCgg                                                       <  2:1776627/148‑1 (MQ=255)
                                                                                      atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCgg                                                       <  1:3421535/148‑1 (MQ=255)
                                                                                      atTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCgg                                                       <  2:3539364/148‑1 (MQ=255)
                                                                                       ttGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGt                                                      >  2:1123401/1‑148 (MQ=255)
                                                                                       ttGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGt                                                      >  2:1158116/1‑148 (MQ=255)
                                                                                       ttGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGt                                                      >  2:1450571/1‑148 (MQ=255)
                                                                                            gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAag                                                 >  2:3219142/1‑148 (MQ=255)
                                                                                            gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAag                                                 >  2:2925802/1‑148 (MQ=255)
                                                                                            gTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAag                                                 >  2:1021079/1‑148 (MQ=255)
                                                                                             ttATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAaga                                                <  1:1246024/148‑1 (MQ=255)
                                                                                              tATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAagag                                               <  2:583708/148‑1 (MQ=255)
                                                                                                tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCg                                             >  1:1032658/1‑148 (MQ=255)
                                                                                                tGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCg                                             >  1:2704619/1‑148 (MQ=255)
                                                                                                 gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgt                                            >  1:2653773/1‑148 (MQ=255)
                                                                                                 gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgt                                            >  1:2545817/1‑148 (MQ=255)
                                                                                                 gAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgt                                            >  1:649632/1‑148 (MQ=255)
                                                                                                  aaaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtc                                           >  1:1134566/1‑148 (MQ=255)
                                                                                                   aaaaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtct                                          >  2:2162422/1‑148 (MQ=255)
                                                                                                     aaGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCgtctgt                                        >  1:2968739/1‑148 (MQ=255)
                                                                                                        cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                     >  2:2831886/1‑148 (MQ=255)
                                                                                                        cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                     <  1:1898729/148‑1 (MQ=255)
                                                                                                        cTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGtctc                                     <  2:1947585/148‑1 (MQ=255)
                                                                                                            tCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTggg                                 >  2:3009981/1‑148 (MQ=255)
                                                                                                               atatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGa                              >  1:281542/1‑148 (MQ=255)
                                                                                                               atatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGa                              >  2:2797566/1‑148 (MQ=255)
                                                                                                                tatCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGAt                             <  1:2937265/148‑1 (MQ=255)
                                                                                                                 atCCGTAACGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATc                            <  2:156665/148‑1 (MQ=255)
                                                                                                                       aaCGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACTGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATa                      <  2:655/148‑1 (MQ=255)
                                                                                                                         cGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAcc                    >  1:1929056/1‑148 (MQ=255)
                                                                                                                         cGGTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAcc                    <  1:1995773/148‑1 (MQ=255)
                                                                                                                          ggTAAAGGCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATAccc                   <  2:2776438/148‑1 (MQ=255)
                                                                                                                                ggCGAAGTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGc             <  2:76903/148‑1 (MQ=255)
                                                                                                                                      gTCGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAgg       <  2:1349533/148‑1 (MQ=255)
                                                                                                                                        cGATGATATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGat     >  2:2004704/1‑148 (MQ=255)
                                                                                                                                           tgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATg  >  2:1816448/1‑148 (MQ=255)
                                                                                                                                           tgatATAGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATg  >  1:48168/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TGTCTGCGGTTGGTCGTATGAAGTTCAACCGTTCTCTGCTGCGCGAAGAAATCGAAGGTTCCGGTATCCTGAGCAAAGACGACATCATTGATGTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATG  >  CP009273/4172362‑4172648

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: ‑