Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,180 A→T 10.7% T473S (ACC→TCC) ‡ rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,1800AT10.7% 645.4 / 51.0 212T473S (ACC→TCC) ‡rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (105/84);  new base T (11/12);  total (116/96)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

TCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGG  >  NC_002947/543031‑543328
                                                                                                                                                     |                                                                                                                                                    
  aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGAcc                                                                                                                                                    <  2:2932088/150‑1 (MQ=255)
          aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                                                                                                            >  2:5144729/1‑150 (MQ=255)
          aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                                                                                                            >  2:4651166/1‑150 (MQ=255)
          aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                                                                                                            >  1:2065677/1‑150 (MQ=255)
          aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                                                                                                            >  2:2580462/1‑150 (MQ=255)
           ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                                                                                                                           >  2:1773843/1‑150 (MQ=255)
           ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                                                                                                                           >  1:1373312/1‑150 (MQ=255)
               cgctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                                                                                                                       >  2:1749073/6‑150 (MQ=255)
                gctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCt                                                                                                                                      >  2:5350005/5‑150 (MQ=255)
                  ttCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCg                                                                                                                                    >  1:3892406/1‑150 (MQ=255)
                   gcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGa                                                                                                                                   >  2:4482745/2‑150 (MQ=255)
                        gTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAgcac                                                                                                                              <  1:5350005/150‑3 (MQ=255)
                         tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                                                                                              <  1:2836670/149‑1 (MQ=255)
                         tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                                                                                              >  2:2836620/1‑149 (MQ=255)
                             ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                                                                                                                         >  2:2898366/1‑150 (MQ=255)
                                 tCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtga                                                                                                                     <  2:4537019/150‑1 (MQ=255)
                                   cAGGCGTTTGAACCCGTACTGATCGAAGCTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                                                                                              <  2:325619/139‑1 (MQ=255)
                                    aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga                                                                                                                  >  1:4445135/1‑150 (MQ=255)
                                     ggCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTgatgat                                                                                                                 <  1:424493/150‑1 (MQ=255)
                                              aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa                                                                                                        <  1:50691/150‑1 (MQ=255)
                                              aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa                                                                                                        >  2:949120/1‑150 (MQ=255)
                                                ccacgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACcaaca                                                                                                      <  2:690928/145‑1 (MQ=255)
                                                ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaaca                                                                                                      >  1:2175035/1‑150 (MQ=255)
                                                  gTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                                                                                                                                  >  2:4621075/1‑120 (MQ=255)
                                                  acgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCaacaac                                                                                                    >  1:1110036/4‑150 (MQ=255)
                                                      tGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaac                                                                                                       <  2:3450750/143‑1 (MQ=255)
                                                        aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGACGGCCGTTCACGTGCCGCTGACCCTGGAAGccc                                                                                                                                          <  1:5407048/106‑1 (MQ=255)
                                                         tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                                                                                                                                  <  1:2141422/113‑1 (MQ=255)
                                                         tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGt                                                                                             >  1:462270/1‑150 (MQ=255)
                                                          cGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTc                                                                                            <  2:1744827/150‑1 (MQ=255)
                                                           gAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCg                                                                                           >  2:973236/1‑150 (MQ=255)
                                                           gAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCg                                                                                           >  1:1842748/1‑150 (MQ=255)
                                                              ggTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgcca                                                                                        <  1:1513780/150‑1 (MQ=255)
                                                                   ggCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCaa                                                                                   >  2:4019365/1‑150 (MQ=255)
                                                                   ggCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                                                                                                                       >  1:4795312/1‑98 (MQ=255)
                                                                    gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                  >  1:5245656/1‑150 (MQ=255)
                                                                      tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg                                                                                >  2:3455468/1‑150 (MQ=255)
                                                                         tCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGa                                                                             <  2:462327/150‑1 (MQ=255)
                                                                         tCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGa                                                                             <  1:4978819/150‑1 (MQ=255)
                                                                           aGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGc                                                                           <  1:2814241/150‑1 (MQ=255)
                                                                           aGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGc                                                                           >  1:272686/1‑150 (MQ=255)
                                                                            gCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGcc                                                                          >  2:4298397/1‑150 (MQ=255)
                                                                              tGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCaa                                                                        <  1:1822174/150‑1 (MQ=255)
                                                                                catccGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtg                                                                                                                      <  2:4283613/99‑1 (MQ=255)
                                                                                catccGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtg                                                                                                                      >  1:4283613/4‑102 (MQ=255)
                                                                                cACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatc                                                                      <  2:1373312/150‑1 (MQ=255)
                                                                                   ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                   >  2:3924164/1‑150 (MQ=255)
                                                                                   ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                   <  1:3509797/150‑1 (MQ=255)
                                                                                       tGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTc                                                               >  2:1543714/1‑150 (MQ=255)
                                                                                         gTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCg                                                             >  2:2720182/1‑150 (MQ=255)
                                                                                         gTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCg                                                             >  2:4109335/1‑150 (MQ=255)
                                                                                               acagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCt                                                                                                                       >  1:4546718/4‑86 (MQ=255)
                                                                                                ctgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAgg                                                      >  2:1074284/1‑150 (MQ=255)
                                                                                                 tgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGa                                                     <  1:740799/150‑1 (MQ=255)
                                                                                                 tgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGa                                                     >  2:1968312/1‑150 (MQ=255)
                                                                                                  gcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGAc                                                    >  2:78927/1‑150 (MQ=255)
                                                                                                   cGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAca                                                   >  2:4759947/1‑149 (MQ=255)
                                                                                                     tACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGtt                                                 <  2:5236690/150‑1 (MQ=255)
                                                                                                      aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                                <  1:4647198/150‑1 (MQ=255)
                                                                                                      aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                                >  1:3308408/1‑150 (MQ=255)
                                                                                                      aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                                >  1:2301451/1‑150 (MQ=255)
                                                                                                       cAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGt                                               >  1:5035630/1‑150 (MQ=255)
                                                                                                        aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                              >  1:822/1‑150 (MQ=255)
                                                                                                        aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                              >  2:1617003/1‑150 (MQ=255)
                                                                                                          cGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACt                                            <  2:146504/150‑1 (MQ=255)
                                                                                                               aCTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAg                                                                            <  2:5309270/113‑1 (MQ=255)
                                                                                                                 ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTg                                     >  1:1309014/1‑150 (MQ=255)
                                                                                                                 ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTg                                     >  2:4916134/1‑150 (MQ=255)
                                                                                                                 ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTg                                     >  1:26313/1‑150 (MQ=255)
                                                                                                                 ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTg                                     >  2:3109995/1‑150 (MQ=255)
                                                                                                                   cGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatg                                                                                                                   <  1:2412369/70‑1 (MQ=255)
                                                                                                                   cGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGta                                   >  2:1663220/1‑150 (MQ=255)
                                                                                                                      cGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtacta                                <  2:3892406/150‑1 (MQ=255)
                                                                                                                       ggTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtactac                               <  2:5424302/150‑1 (MQ=255)
                                                                                                                       ggTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtactac                               <  1:293974/150‑1 (MQ=255)
                                                                                                                       ggTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtactac                               <  1:316989/150‑1 (MQ=255)
                                                                                                                        gTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACa                              <  1:2671548/150‑1 (MQ=255)
                                                                                                                          gACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATg                            >  1:4234414/1‑150 (MQ=255)
                                                                                                                             cAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGAcc                         >  1:2080618/1‑150 (MQ=255)
                                                                                                                             cAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGAcc                         >  2:2850549/1‑150 (MQ=255)
                                                                                                                               gATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaacaa                                                                                                     >  2:1687237/1‑72 (MQ=255)
                                                                                                                               gATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaacaa                                                                                                     <  1:1687237/72‑1 (MQ=255)
                                                                                                                                aTGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCgg                      <  2:2149421/150‑1 (MQ=255)
                                                                                                                                aTGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCgg                      >  1:397993/1‑150 (MQ=255)
                                                                                                                                aTGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCgg                      >  1:565737/1‑150 (MQ=255)
                                                                                                                                   gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                   <  1:5443944/150‑1 (MQ=255)
                                                                                                                                     cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                 >  1:5308123/1‑150 (MQ=255)
                                                                                                                                      gTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGcc                <  2:440392/150‑1 (MQ=255)
                                                                                                                                        tCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCAt              <  1:344879/150‑1 (MQ=255)
                                                                                                                                           cGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCaa           >  1:4087337/1‑150 (MQ=255)
                                                                                                                                            gTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAAt          <  1:1663220/150‑1 (MQ=255)
                                                                                                                                             tGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATg         <  2:397993/150‑1 (MQ=255)
                                                                                                                                              gCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGc        <  1:584507/150‑1 (MQ=255)
                                                                                                                                                 gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGta                                   <  2:1601139/120‑1 (MQ=255)
                                                                                                                                                 gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCTaa     >  2:1926478/1‑150 (MQ=255)
                                                                                                                                                 gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCaa     <  2:4817071/150‑1 (MQ=255)
                                                                                                                                                 gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCaa     <  2:2080618/150‑1 (MQ=255)
                                                                                                                                                    gACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATg         <  1:2744819/143‑1 (MQ=255)
                                                                                                                                                    gACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATg         >  2:2744819/1‑143 (MQ=255)
                                                                                                                                                    gACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAggg  >  2:48789/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGG  >  NC_002947/543031‑543328

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: