Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,206 T→C 10.7% R481R (CGT→CGC 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,2060TC10.7% 527.5 / 38.4 188R481R (CGT→CGCrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (102/65);  new base C (13/7);  total (115/73)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 193 total aligned reads displayed.

CTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCA  >  NC_002947/543057‑543355
                                                                                                                                                     |                                                                                                                                                     
cTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGt                                                                                                                                                       <  1:2060776/150‑1 (MQ=255)
 tGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTg                                                                                                                                                      >  2:2313911/1‑150 (MQ=255)
 tGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTg                                                                                                                                                      >  2:4279528/1‑150 (MQ=255)
 tGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTg                                                                                                                                                      >  2:3077956/1‑150 (MQ=255)
        tcAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtgat                                                                                                                                               <  2:1812285/149‑1 (MQ=255)
        ccAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgat                                                                                                                                               <  2:82345/150‑1 (MQ=255)
          aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga                                                                                                                                             >  2:5052315/1‑150 (MQ=255)
           ggCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTgatgat                                                                                                                                            <  2:723466/150‑1 (MQ=255)
             cGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGt                                                                                                                                          <  1:585077/150‑1 (MQ=255)
              gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc                                                                                                                                         <  1:2836770/150‑1 (MQ=255)
                 tGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGAc                                                                                                                                      <  2:3213298/150‑1 (MQ=255)
                 tGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGAc                                                                                                                                      <  2:3213291/150‑1 (MQ=255)
                  gAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGa                                                                                                                                       <  1:1265795/148‑1 (MQ=255)
                  gAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGa                                                                                                                                       >  2:1265795/1‑148 (MQ=255)
                      ccacgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACcaaca                                                                                                                                 >  1:2216421/6‑150 (MQ=255)
                      ccacgcTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACcaaca                                                                                                                                 >  1:3534313/6‑150 (MQ=255)
                      ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaaca                                                                                                                                 >  1:2679137/1‑150 (MQ=255)
                      ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaaca                                                                                                                                 >  1:421113/1‑150 (MQ=255)
                         tACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAacaaca                                                                                                                              <  2:2097300/150‑1 (MQ=255)
                           cTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATc                                                                                                                            <  2:1057216/150‑1 (MQ=255)
                              aTCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTg                                                                                                                         >  1:4898755/1‑150 (MQ=255)
                               tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGt                                                                                                                        >  2:3536553/1‑150 (MQ=255)
                               tCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGt                                                                                                                        >  1:4098232/1‑150 (MQ=255)
                               tCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTgc                                                                                                                        >  1:209003/1‑149 (MQ=255)
                                 gAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCg                                                                                                                      >  1:3827584/1‑150 (MQ=255)
                                 gAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCg                                                                                                                      >  1:4334955/1‑150 (MQ=255)
                                  aaggcaaggcCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTgcacc                                                                                                                     >  2:2143742/1‑146 (MQ=255)
                                   aGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTgcaccc                                                                                                                    >  2:687382/1‑145 (MQ=255)
                                     gTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccag                                                                                                                  >  1:3174095/1‑150 (MQ=255)
                                     gTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccag                                                                                                                  >  1:3775166/1‑150 (MQ=255)
                                      tAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagc                                                                                                                 >  2:1962522/1‑150 (MQ=255)
                                       aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                >  2:2792034/1‑150 (MQ=255)
                                       aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                >  1:2341741/1‑150 (MQ=255)
                                       aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                >  2:4429324/1‑150 (MQ=255)
                                       aaGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCgccagcc                                                                                                                >  1:1856804/1‑150 (MQ=255)
                                          gCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAAc                                                                                                             >  1:2571770/1‑150 (MQ=255)
                                            tATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACgg                                                                                                           >  1:3903502/1‑150 (MQ=255)
                                             aTTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGt                                                                                                          >  1:3887497/1‑150 (MQ=255)
                                                       aCCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatca                                                                                                <  2:4899210/150‑1 (MQ=255)
                                                         ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                              <  1:110285/150‑1 (MQ=255)
                                                         ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                              >  1:4218743/1‑150 (MQ=255)
                                                         ccGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAatcatc                                                                                              >  1:2508502/1‑150 (MQ=255)
                                                         ccGCTCGTCTGTACAGAGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGatcatc                                                                                              <  2:4624354/150‑1 (MQ=255)
                                                               gTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCg                                                                                        >  2:4950606/1‑150 (MQ=255)
                                                                  tgtacagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCg                                                                                     >  2:2706131/7‑150 (MQ=255)
                                                                   gtacagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGc                                                                                    <  1:2706131/145‑1 (MQ=255)
                                                                      cagcGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAgg                                                                                 >  2:516437/3‑150 (MQ=255)
                                                                       tgcGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGa                                                                                >  1:1058247/1‑150 (MQ=255)
                                                                        gcGTACAACGCCGACTTCGACGGTGACAAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGAc                                                                               >  2:4206142/1‑150 (MQ=255)
                                                                          gTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGt                                                                             >  2:1749610/1‑150 (MQ=255)
                                                                          gTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGt                                                                             >  2:1280884/1‑150 (MQ=255)
                                                                           tACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGtt                                                                            <  1:1440/150‑1 (MQ=255)
                                                                           tACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGtt                                                                            <  1:3609758/150‑1 (MQ=255)
                                                                            tcAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCTGATGATGTCGACCAACAACATCCTGCACCCGGCAAACGGCGCACCGATCATCGTTCCTTCGCAGGAcatgg                                                                           <  1:2299223/149‑5 (MQ=255)
                                                                            aCAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTg                                                                           >  2:4439787/1‑150 (MQ=255)
                                                                             ccaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACACCCTCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGt                                                                          <  1:3401393/148‑1 (MQ=255)
                                                                             cAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgc                                                                                                                                                     >  2:3258469/1‑75 (MQ=255)
                                                                             cAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgc                                                                                                                                                     <  1:3258469/75‑1 (MQ=255)
                                                                              aaCGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTa                                                                         >  1:4144696/1‑150 (MQ=255)
                                                                                cGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACt                                                                       >  2:1550789/1‑150 (MQ=255)
                                                                                 gCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTg                                                                      <  2:4374418/150‑1 (MQ=255)
                                                                                 gCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTg                                                                      <  2:3019393/150‑1 (MQ=255)
                                                                                  ccGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTgg                                                                     >  2:2626019/1‑150 (MQ=255)
                                                                                   cGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTggg                                                                    <  1:1627235/150‑1 (MQ=255)
                                                                                   cGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTggg                                                                    >  1:3804275/1‑150 (MQ=255)
                                                                                    gACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGt                                                                   >  1:1603434/1‑150 (MQ=255)
                                                                                    gACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGt                                                                   <  2:3687502/150‑1 (MQ=255)
                                                                                       ttCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTg                                                                >  1:919176/1‑150 (MQ=255)
                                                                                          gACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGtac                                                             >  1:746914/1‑150 (MQ=255)
                                                                                                     gATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCg                                                  >  1:670344/1‑150 (MQ=255)
                                                                                                      aTGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCgg                                                 >  1:284334/1‑150 (MQ=255)
                                                                                                       tGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCggg                                                >  2:337646/1‑150 (MQ=255)
                                                                                                         gCCTTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTTAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                              <  2:4817673/150‑1 (MQ=255)
                                                                                                         gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                              <  1:3510644/150‑1 (MQ=255)
                                                                                                         gCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGaa                                              <  1:3958509/150‑1 (MQ=255)
                                                                                                           cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                            <  1:514686/150‑1 (MQ=255)
                                                                                                           cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                            >  2:5074572/1‑150 (MQ=255)
                                                                                                           cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                            >  1:3442844/1‑150 (MQ=255)
                                                                                                           cGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGc                                            >  2:480946/1‑150 (MQ=255)
                                                                                                            gTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGcc                                           >  2:3797981/1‑150 (MQ=255)
                                                                                                                aCGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCa                                       >  2:2262037/1‑150 (MQ=255)
                                                                                                                    gCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGc                                   <  2:4714279/150‑1 (MQ=255)
                                                                                                                       gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCaa                                <  2:4098232/150‑1 (MQ=255)
                                                                                                                       gCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCaa                                <  2:3433725/150‑1 (MQ=255)
                                                                                                                          gACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAggg                             <  2:421113/150‑1 (MQ=255)
                                                                                                                           aCCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGc                            >  1:2012334/1‑150 (MQ=255)
                                                                                                                             ccTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGa                          >  2:4798558/1‑150 (MQ=255)
                                                                                                                             ccTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGa                          >  2:1386573/1‑150 (MQ=255)
                                                                                                                              cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                         <  2:1999957/150‑1 (MQ=255)
                                                                                                                              cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                         <  1:2167824/150‑1 (MQ=255)
                                                                                                                              cTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGaa                         <  1:888464/150‑1 (MQ=255)
                                                                                                                               tGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAg                        >  2:2232723/1‑150 (MQ=255)
                                                                                                                                    gCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTcgc                   <  2:577006/150‑1 (MQ=255)
                                                                                                                                      ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                 >  1:1526783/1‑150 (MQ=255)
                                                                                                                                      ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                 >  1:461444/1‑150 (MQ=255)
                                                                                                                                      ccAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCgt                 <  1:4798558/150‑1 (MQ=255)
                                                                                                                                         gCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGtt              >  1:2097693/1‑150 (MQ=255)
                                                                                                                                          cTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTc             <  1:2294257/150‑1 (MQ=255)
                                                                                                                                                     tGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCa  <  1:5052315/150‑1 (MQ=255)
                                                                                                                                                     tGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCa  <  1:741244/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAACATCCTGTCGCCAGCCAACGGTGAGCCAATCATCGTTCCGTCGCAGGACGTTGTACTGGGTCTGTACTACATGACCCGGGAAGCCATCAATGCCAAGGGCGAAGGTCGCGTGTTCGCCGACCTGCA  >  NC_002947/543057‑543355

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: