Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,173,145 A→T 100% Q658L (CAG→CTG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,173,1450AT97.7% 717.5 / 4.3 221Q658L (CAG→CTG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base T (109/107);  total (110/111)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 226 total aligned reads displayed.

TCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACG  >  CP009273/4172998‑4173289
                                                                                                                                                   |                                                                                                                                                
tCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCca                                                                                                                                                  >  1:707153/1‑148 (MQ=255)
   aagaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgca                                                                                                                                               <  1:2651554/148‑1 (MQ=255)
   aagaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgca                                                                                                                                               >  1:3100488/1‑148 (MQ=255)
   aagaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgca                                                                                                                                               >  2:2613645/1‑148 (MQ=255)
         gCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                         <  1:2012260/148‑1 (MQ=255)
         gCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                         <  1:2011822/148‑1 (MQ=255)
           aaCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTAt                                                                                                                                       <  2:3064661/148‑1 (MQ=255)
           aaCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTAt                                                                                                                                       <  2:1290333/148‑1 (MQ=255)
            aCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATc                                                                                                                                      >  1:1657781/1‑148 (MQ=255)
             cTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                     >  2:1065907/1‑148 (MQ=255)
             cTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                     >  1:3266721/1‑148 (MQ=255)
             cTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                     >  2:2846855/1‑148 (MQ=255)
              tACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCg                                                                                                                                    >  2:1396416/1‑148 (MQ=255)
              tACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCg                                                                                                                                    <  1:1909518/148‑1 (MQ=255)
               aCGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGt                                                                                                                                   <  2:2396855/148‑1 (MQ=255)
               aCGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGt                                                                                                                                   <  1:2741303/148‑1 (MQ=255)
                cGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTc                                                                                                                                  <  2:219113/148‑1 (MQ=255)
                    aTCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTg                                                                                                                              <  1:1350960/148‑1 (MQ=255)
                       gCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGt                                                                                                                           <  2:2308674/148‑1 (MQ=255)
                         ccAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTcc                                                                                                                         <  2:1964475/148‑1 (MQ=255)
                          cAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTccc                                                                                                                        <  2:2791274/148‑1 (MQ=255)
                          cAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTccc                                                                                                                        <  2:2832205/148‑1 (MQ=255)
                            ggCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTg                                                                                                                      >  2:871449/1‑148 (MQ=255)
                              cGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGAt                                                                                                                    <  1:3081864/148‑1 (MQ=255)
                                aaCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATcc                                                                                                                  <  2:1330847/148‑1 (MQ=255)
                                 aCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccc                                                                                                                 >  2:3213864/1‑148 (MQ=255)
                                  cTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccct                                                                                                                >  2:3160445/1‑147 (MQ=255)
                                     cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTc                                                                                                             <  1:1517517/148‑1 (MQ=255)
                                     cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTc                                                                                                             <  1:3457620/148‑1 (MQ=255)
                                     cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTc                                                                                                             >  1:2993572/1‑148 (MQ=255)
                                      aaCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTcc                                                                                                            <  1:3542058/148‑1 (MQ=255)
                                       aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCt                                                                                                           <  1:1962421/148‑1 (MQ=255)
                                       aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCt                                                                                                           >  2:2812189/1‑148 (MQ=255)
                                       aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCt                                                                                                           >  1:2174108/1‑148 (MQ=255)
                                        cTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTg                                                                                                          <  1:2474196/148‑1 (MQ=255)
                                          tGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGa                                                                                                        <  2:2877135/148‑1 (MQ=255)
                                           ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGaa                                                                                                       <  2:609525/148‑1 (MQ=255)
                                            gATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAac                                                                                                      >  1:867443/1‑148 (MQ=255)
                                             aTGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAaca                                                                                                     >  1:471829/1‑148 (MQ=255)
                                                 agaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGAt                                                                                                 >  1:2842079/1‑148 (MQ=255)
                                                  gaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATg                                                                                                >  2:1391903/1‑148 (MQ=255)
                                                   aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                               >  2:2693249/1‑148 (MQ=255)
                                                   aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                               <  2:1657781/148‑1 (MQ=255)
                                                   aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                               >  2:1309567/1‑148 (MQ=255)
                                                          cTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAAc                                                                                        >  2:2050533/1‑148 (MQ=255)
                                                          cTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAAc                                                                                        >  2:1890652/1‑148 (MQ=255)
                                                           ttCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAAcc                                                                                       >  1:1018061/1‑148 (MQ=255)
                                                               tAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGc                                                                                   <  1:2959566/148‑1 (MQ=255)
                                                                 gaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                 >  2:736110/1‑148 (MQ=255)
                                                                  aagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGccgt                                                                                >  2:1918557/1‑145 (MQ=255)
                                                                  aagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAtt                                                                                >  2:357208/1‑148 (MQ=255)
                                                                  aagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAtt                                                                                >  2:2165810/1‑148 (MQ=255)
                                                                  aagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAtt                                                                                <  2:1883092/148‑1 (MQ=255)
                                                                   agaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTg                                                                               >  2:1586455/1‑148 (MQ=255)
                                                                       cTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATgg                                                                           <  1:355802/148‑1 (MQ=255)
                                                                            aaCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGAATTGATGGGTGCg                                                                      >  1:1813528/1‑148 (MQ=255)
                                                                             aCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGa                                                                     <  2:2857699/148‑1 (MQ=255)
                                                                              cTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                                    >  1:1847379/1‑148 (MQ=255)
                                                                                    gTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                              <  1:1058435/148‑1 (MQ=255)
                                                                                     tAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                             >  2:2656131/1‑148 (MQ=255)
                                                                                        cAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                          >  2:2894416/1‑148 (MQ=255)
                                                                                            ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                      >  2:1112362/1‑148 (MQ=255)
                                                                                            ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                      <  2:451078/148‑1 (MQ=255)
                                                                                            ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                      >  2:3456089/1‑148 (MQ=255)
                                                                                               gAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCg                                                   <  2:348638/148‑1 (MQ=255)
                                                                                               gAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCg                                                   <  2:1676465/148‑1 (MQ=255)
                                                                                                       cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                           >  2:3121385/1‑148 (MQ=255)
                                                                                                       cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                           <  2:2441908/148‑1 (MQ=255)
                                                                                                                ccGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCt                                  <  1:675403/148‑1 (MQ=255)
                                                                                                                 cgcgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTg                                 <  1:500022/148‑1 (MQ=255)
                                                                                                                 cgcgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTg                                 >  1:743074/1‑148 (MQ=255)
                                                                                                                  gcgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGa                                >  2:434693/1‑148 (MQ=255)
                                                                                                                     aCCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                             >  2:56054/1‑148 (MQ=255)
                                                                                                                      ccAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                            >  1:3469760/1‑148 (MQ=255)
                                                                                                                      ccAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                            >  1:1182510/1‑148 (MQ=255)
                                                                                                                         ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                         >  1:468658/1‑148 (MQ=255)
                                                                                                                         ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                         <  1:1428988/148‑1 (MQ=255)
                                                                                                                          gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                        >  1:984859/1‑148 (MQ=255)
                                                                                                                          gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                        >  2:1484295/1‑148 (MQ=255)
                                                                                                                             gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg                     <  1:1309567/148‑1 (MQ=255)
                                                                                                                             gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg                     <  1:2693249/148‑1 (MQ=255)
                                                                                                                              aCTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                    >  2:2967790/1‑148 (MQ=255)
                                                                                                                                tACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttg                  <  1:1890598/148‑1 (MQ=255)
                                                                                                                                   aTGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa               <  1:1396416/148‑1 (MQ=255)
                                                                                                                                   aTGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa               >  1:15822/1‑148 (MQ=255)
                                                                                                                                    tGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc              <  2:2735248/148‑1 (MQ=255)
                                                                                                                                    tGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc              >  2:426887/1‑148 (MQ=255)
                                                                                                                                     ggACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt             >  1:3225553/1‑148 (MQ=255)
                                                                                                                                     ggACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt             >  1:1454138/1‑148 (MQ=255)
                                                                                                                                     ggACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTACGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt             >  1:976542/1‑148 (MQ=255)
                                                                                                                                         gTATCCACCCTGCAGGTGGTATCCGTCGGTGTGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa         <  1:2371446/148‑1 (MQ=255)
                                                                                                                                         gTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa         <  1:3161814/148‑1 (MQ=255)
                                                                                                                                          tATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTAt        <  1:3328952/148‑1 (MQ=255)
                                                                                                                                           aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATg       <  1:514793/148‑1 (MQ=255)
                                                                                                                                            tCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATgg      >  2:2330243/1‑148 (MQ=255)
                                                                                                                                             ccACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGa     <  2:1847379/148‑1 (MQ=255)
                                                                                                                                             ccACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGa     <  2:1198468/148‑1 (MQ=255)
                                                                                                                                             ccACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGa     >  1:987548/1‑148 (MQ=255)
                                                                                                                                              cACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGaa    <  1:3121385/148‑1 (MQ=255)
                                                                                                                                                cccTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg  <  1:1110530/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACG  >  CP009273/4172998‑4173289

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