Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,101 T→C 12.0% A446A (GCT→GCC 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,1010TC12.0% 457.9 / 39.9 167A446A (GCT→GCCrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (78/69);  new base C (10/10);  total (88/79)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 169 total aligned reads displayed.

CAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGG  >  NC_002947/542954‑543250
                                                                                                                                                   |                                                                                                                                                     
cAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAt                                                                                                                                                     >  2:199609/1‑150 (MQ=255)
cAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGAACCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAt                                                                                                                                                     <  2:612303/150‑1 (MQ=255)
 aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                                                    >  2:2641480/1‑150 (MQ=255)
 aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                                                    <  1:4292537/150‑1 (MQ=255)
  agaagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTc                                                                                                                                                   >  2:4978819/1‑150 (MQ=255)
   gaagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCa                                                                                                                                                  >  2:666674/1‑150 (MQ=255)
      gaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                                                               <  1:3478017/150‑1 (MQ=255)
      gaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                                                               <  1:2603024/150‑1 (MQ=255)
        tGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGc                                                                                                                                             >  1:4844075/1‑150 (MQ=255)
             gAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCg                                                                                                                                        >  1:468047/1‑150 (MQ=255)
              aGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGc                                                                                                                                       <  2:2114693/150‑1 (MQ=255)
               gcgcGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCt                                                                                                                                      >  2:344879/1‑150 (MQ=255)
                cgcgAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTg                                                                                                                                     <  1:3898103/150‑1 (MQ=255)
                  cgAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGt                                                                                                                                   >  1:1252228/1‑150 (MQ=255)
                     gCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCtg                                                                                                                                <  1:666674/150‑1 (MQ=255)
                      cTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCtgt                                                                                                                               <  2:2393797/150‑1 (MQ=255)
                         ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgct                                                                                                                            <  1:2390612/150‑1 (MQ=255)
                           agagGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                                                          <  2:2956380/150‑1 (MQ=255)
                               gtgtGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTAc                                                                                                                      >  2:5157998/1‑150 (MQ=255)
                                 gtgGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACaa                                                                                                                    <  2:933172/150‑1 (MQ=255)
                                     gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                                                <  1:4284630/150‑1 (MQ=255)
                                     gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                                                <  2:1000074/150‑1 (MQ=255)
                                      aCGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCg                                                                                                               <  1:320509/150‑1 (MQ=255)
                                      aCGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCg                                                                                                               <  1:320535/150‑1 (MQ=255)
                                          tctcGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACtt                                                                                                           <  2:1073217/150‑1 (MQ=255)
                                           ctcGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGGTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGc                                                                                                                                                >  2:4480831/1‑112 (MQ=255)
                                            tcGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCg                                                                                                         >  2:4005353/1‑150 (MQ=255)
                                               cTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACg                                                                                                      <  1:3408576/150‑1 (MQ=255)
                                                 gAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGt                                                                                                    >  1:3763299/1‑150 (MQ=255)
                                                        ttCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACCACGCCGACTTCGACGGTGACCAGa                                                                                             >  1:2487981/1‑150 (MQ=255)
                                                         tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGAt                                                                                            <  1:5157998/150‑1 (MQ=255)
                                                         tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGAt                                                                                            <  2:3547321/150‑1 (MQ=255)
                                                         tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGAt                                                                                            >  2:5443944/1‑150 (MQ=255)
                                                           gcgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATgg                                                                                          <  2:3231337/150‑1 (MQ=255)
                                                           gcgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATgg                                                                                          >  2:584507/1‑150 (MQ=255)
                                                                  ccccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCa                                                                                   <  2:2877739/150‑1 (MQ=255)
                                                                   cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                                                  >  2:2740346/1‑150 (MQ=255)
                                                                   cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                                                  >  1:146504/1‑150 (MQ=255)
                                                                    ccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACg                                                                                 >  2:117462/1‑150 (MQ=255)
                                                                    ccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACg                                                                                 >  2:117429/1‑150 (MQ=255)
                                                                     cGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGt                                                                                >  1:1473806/1‑150 (MQ=255)
                                                                         cttctgaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTcct                                                                            <  2:2783369/144‑2 (MQ=255)
                                                                          tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTCAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                                                           >  2:1140347/1‑150 (MQ=255)
                                                                          ctcTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                                                           <  2:3266299/148‑1 (MQ=255)
                                                                            cTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTg                                                                               <  1:263985/144‑1 (MQ=255)
                                                                             tCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGa                                                                        <  1:2068516/150‑1 (MQ=255)
                                                                              gaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTgtc                                                                       <  1:2236796/149‑3 (MQ=255)
                                                                              gaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTgtc                                                                       >  1:2011992/2‑148 (MQ=255)
                                                                              gaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTgtc                                                                       <  1:181501/149‑3 (MQ=255)
                                                                               aaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGAt                                                                                            <  1:2044473/128‑1 (MQ=255)
                                                                               aaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGAt                                                                                            >  2:2044473/1‑128 (MQ=255)
                                                                               aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGAcc                                                                      <  2:2932088/150‑1 (MQ=255)
                                                                                       aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                              >  1:2065677/1‑150 (MQ=255)
                                                                                       aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                                              >  2:2580462/1‑150 (MQ=255)
                                                                                        ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                                             >  1:1373312/1‑150 (MQ=255)
                                                                                        ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                                             >  2:1773843/1‑150 (MQ=255)
                                                                                            cgctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                                         >  2:1749073/6‑150 (MQ=255)
                                                                                             gctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCt                                                        >  2:5350005/5‑150 (MQ=255)
                                                                                               ttCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCg                                                      >  1:5424302/1‑150 (MQ=255)
                                                                                                     gTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAgcac                                                <  1:5350005/150‑3 (MQ=255)
                                                                                                      tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                <  1:2836620/149‑1 (MQ=255)
                                                                                                      tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                <  1:2836670/149‑1 (MQ=255)
                                                                                                      tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                                                >  2:2836620/1‑149 (MQ=255)
                                                                                                          ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                                           >  2:534798/1‑150 (MQ=255)
                                                                                                          ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                                           >  2:2898145/1‑150 (MQ=255)
                                                                                                           gTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGcacgcgtgc                                          <  1:4120652/150‑6 (MQ=255)
                                                                                                           gTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcgc                                          <  2:1654598/150‑1 (MQ=255)
                                                                                                             attcAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtg                                        >  2:3452350/4‑150 (MQ=255)
                                                                                                              tCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtga                                       <  2:4537019/150‑1 (MQ=255)
                                                                                                               ccAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgat                                      <  1:315475/150‑1 (MQ=255)
                                                                                                                cAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatg                                     >  1:5236690/1‑150 (MQ=255)
                                                                                                                cAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatg                                     >  1:2149421/1‑150 (MQ=255)
                                                                                                                 aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga                                    >  1:4445135/1‑150 (MQ=255)
                                                                                                                  ggCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACg                                                                                                      >  2:2060979/1‑83 (MQ=255)
                                                                                                                      tttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCg                               >  2:854721/1‑150 (MQ=255)
                                                                                                                      tttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCg                               >  1:2351567/1‑150 (MQ=255)
                                                                                                                       ttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGa                              <  2:4198731/150‑1 (MQ=255)
                                                                                                                         gAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGAcc                            >  2:4647198/1‑150 (MQ=255)
                                                                                                                         gAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGAcc                            <  1:2664960/150‑1 (MQ=255)
                                                                                                                           acccacgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACcaa                          <  2:4989513/143‑1 (MQ=255)
                                                                                                                           aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa                          >  2:949204/1‑150 (MQ=255)
                                                                                                                           aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa                          <  1:50691/150‑1 (MQ=255)
                                                                                                                             ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAg                                                                                              >  2:4024463/1‑80 (MQ=255)
                                                                                                                             ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAg                                                                                              <  1:4024463/80‑1 (MQ=255)
                                                                                                                               gTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                                                    >  2:4621075/1‑120 (MQ=255)
                                                                                                                               gTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                                                    <  1:4621075/120‑1 (MQ=255)
                                                                                                                               acgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCaacaac                      >  1:1110036/4‑150 (MQ=255)
                                                                                                                                   tGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaac                         >  1:3450750/1‑143 (MQ=255)
                                                                                                                                   tGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATc                   >  2:423733/1‑149 (MQ=255)
                                                                                                                                   tGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATc                   <  1:423733/149‑1 (MQ=255)
                                                                                                                                     aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGACGGCCGTTCACGTGCCGCTGACCCTGGAAGccc                                                            <  1:5407048/106‑1 (MQ=255)
                                                                                                                                      tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGt               >  1:462327/1‑150 (MQ=255)
                                                                                                                                       cGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTc              <  2:1744854/150‑1 (MQ=255)
                                                                                                                                       cGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTc              <  1:4525359/150‑1 (MQ=255)
                                                                                                                                        gAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCg             >  2:973236/1‑150 (MQ=255)
                                                                                                                                        gAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCg             >  1:1842748/1‑150 (MQ=255)
                                                                                                                                           tgTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgcca          <  2:3985041/149‑1 (MQ=255)
                                                                                                                                                ggCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                                         >  1:4795312/1‑98 (MQ=255)
                                                                                                                                                ggCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                                         <  2:4795312/98‑1 (MQ=255)
                                                                                                                                                   tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg  >  2:3455468/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGG  >  NC_002947/542954‑543250

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: