Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,102 C→T 5.0% H113H (CAC→CAT rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947542,1020CT 5.0% 636.3 / 13.6 199H113H (CAC→CATrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (94/95);  new base T (6/4);  total (100/99)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

CCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGAACGACG  >  NC_002947/541953‑542250
                                                                                                                                                     |                                                                                                                                                    
ccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCcac                                                                                                                                                      >  2:1489356/1‑150 (MQ=255)
  aGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAt                                                                                                                                                    <  2:1453188/150‑1 (MQ=255)
    tCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCt                                                                                                                                                  >  1:1819052/1‑150 (MQ=255)
      aaGGACTACGAGTGCCTGTGCGGCACGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                                                >  1:4271732/1‑150 (MQ=255)
      aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACCTCTgg                                                                                                                                                >  1:2723120/1‑150 (MQ=255)
      aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                                                <  2:3244329/150‑1 (MQ=255)
      aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                                                >  1:1409322/1‑150 (MQ=255)
      aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                                                >  1:1409342/1‑150 (MQ=255)
       aGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGt                                                                                                                                               >  2:1377240/1‑150 (MQ=255)
         gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                                             >  2:3829967/1‑150 (MQ=255)
         gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                                             <  2:4785147/150‑1 (MQ=255)
         gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                                             <  2:5258731/150‑1 (MQ=255)
          aCTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTcc                                                                                                                                            <  2:1006815/150‑1 (MQ=255)
            tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                                                                                                                          >  1:953933/1‑150 (MQ=255)
              tgaatGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGaa                                                                                                                                        >  1:1778378/5‑150 (MQ=255)
              cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                                                                                                        <  1:4829244/150‑1 (MQ=255)
              cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                                                                                                        >  2:1267651/1‑150 (MQ=255)
               gaatGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAg                                                                                                                                       >  1:918944/4‑150 (MQ=255)
               gAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAg                                                                                                                                       >  2:4211373/1‑150 (MQ=255)
               gAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAg                                                                                                                                       >  2:3646930/1‑150 (MQ=255)
                   gCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGc                                                                                                                                   <  1:3247122/150‑1 (MQ=255)
                     cTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTg                                                                                                                                 <  1:4406910/150‑1 (MQ=255)
                     cTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTg                                                                                                                                 <  1:1246145/150‑1 (MQ=255)
                           ggCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTcc                                                                                                                           <  1:4921872/150‑1 (MQ=255)
                            gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCATATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                                                                                                          >  1:1268839/1‑150 (MQ=255)
                            gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                                                                                                          >  1:4915232/1‑150 (MQ=255)
                              aaGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGt                                                                                                                        >  2:4714640/1‑150 (MQ=255)
                                         ccTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCt                                                                                                             >  2:1822044/1‑150 (MQ=255)
                                         ccTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCt                                                                                                             >  1:4189230/1‑150 (MQ=255)
                                          cTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTg                                                                                                            >  2:3133327/1‑150 (MQ=255)
                                           tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa                                                                                                           >  2:1182796/1‑150 (MQ=255)
                                           tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa                                                                                                           >  2:3923388/1‑150 (MQ=255)
                                           tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAGCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa                                                                                                           >  2:3303981/1‑150 (MQ=255)
                                             aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                                                                                                         >  2:3390688/1‑150 (MQ=255)
                                             aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                                                                                                         >  2:975218/1‑150 (MQ=255)
                                                      ggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGAcc                                                                                                <  1:654012/150‑1 (MQ=255)
                                                        cgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCt                                                                                              >  1:4694145/1‑150 (MQ=255)
                                                            aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATcgc                                                                                                                          <  1:2128623/118‑3 (MQ=255)
                                                            aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCtcaag                                                                                          <  1:3035398/150‑5 (MQ=255)
                                                               tGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGat                                                                                       >  2:1252540/1‑150 (MQ=255)
                                                                gCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGata                                                                                      <  1:583762/150‑1 (MQ=255)
                                                                gCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGata                                                                                      <  2:4435336/150‑1 (MQ=255)
                                                                       gTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGc                                                                                                              >  2:2106026/1‑119 (MQ=255)
                                                                        tGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAgcgc                                                                              >  1:855222/1‑150 (MQ=255)
                                                                              gTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGctc                                                                        >  1:1491481/1‑150 (MQ=255)
                                                                               ttGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                                             <  2:971628/80‑1 (MQ=255)
                                                                               ttGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                                             >  1:971641/1‑80 (MQ=255)
                                                                               ttGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGctct                                                                       >  2:3170149/1‑150 (MQ=255)
                                                                                 gAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCCCATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTAc                                                                     <  2:2945168/150‑1 (MQ=255)
                                                                                      tGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGctctcgaagta                                                                <  2:453381/150‑8 (MQ=255)
                                                                                       gCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgag                                                               <  2:4146163/150‑1 (MQ=255)
                                                                                        cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGATCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                                              <  2:2219427/150‑1 (MQ=255)
                                                                                        cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                                              <  1:2272304/150‑1 (MQ=255)
                                                                                        cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                                              <  2:2467225/150‑1 (MQ=255)
                                                                                        cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                                              <  2:440127/150‑1 (MQ=255)
                                                                                           tGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCATGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCt                                                           <  1:1154073/150‑1 (MQ=255)
                                                                                               gAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGt                                                       >  1:3818962/1‑150 (MQ=255)
                                                                                                   gTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCTCAAGGATATCGAGCGCGTTCTGTATTTTGAAAACTacattgtt                                                   >  1:3000334/1‑143 (MQ=255)
                                                                                                   gTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGtt                                                   <  1:2003284/150‑1 (MQ=255)
                                                                                                   gTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGtt                                                   <  1:1320561/150‑1 (MQ=255)
                                                                                                    ttCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTa                                                  <  1:1577171/150‑1 (MQ=255)
                                                                                                    ttCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTa                                                  <  1:4510155/150‑1 (MQ=255)
                                                                                                     tcgtcgTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTAt                                                 <  1:2053098/150‑1 (MQ=255)
                                                                                                     tcgtcgTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTAt                                                 <  1:2316846/150‑1 (MQ=255)
                                                                                                     tcgtcgTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTAt                                                 >  2:2619811/1‑150 (MQ=255)
                                                                                                      cgtcgtGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATc                                                <  2:4950282/150‑1 (MQ=255)
                                                                                                        tcgtGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGa                                              >  1:2326517/1‑150 (MQ=255)
                                                                                                           tGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                                              <  2:1910428/131‑1 (MQ=255)
                                                                                                           tGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                                              >  1:1910428/1‑131 (MQ=255)
                                                                                                             aGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACCTGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCgg                                         >  2:3652312/1‑150 (MQ=255)
                                                                                                             aGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCgg                                         <  2:4716682/150‑1 (MQ=255)
                                                                                                              gcgcATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATc                                                <  1:1229363/142‑1 (MQ=255)
                                                                                                                 cctGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCAt                                     >  2:5376814/3‑150 (MQ=255)
                                                                                                                 cATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCAt                                     <  2:771880/150‑1 (MQ=255)
                                                                                                                 cATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCAt                                     <  2:1883286/150‑1 (MQ=255)
                                                                                                                 cATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCAt                                     >  2:3692991/1‑150 (MQ=255)
                                                                                                                  aTGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATg                                    >  2:5035546/1‑150 (MQ=255)
                                                                                                                   tGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGa                                   <  1:5276080/150‑1 (MQ=255)
                                                                                                                   tGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGa                                   <  1:2549163/150‑1 (MQ=255)
                                                                                                                        cacaTCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTAcc                              >  2:4140666/1‑150 (MQ=255)
                                                                                                                            tCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTgg                          <  1:1370938/150‑1 (MQ=255)
                                                                                                                             cGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGa                         <  1:3646930/150‑1 (MQ=255)
                                                                                                                                  tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg                    <  2:982902/150‑1 (MQ=255)
                                                                                                                                  tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg                    <  1:1198517/150‑1 (MQ=255)
                                                                                                                                  tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg                    <  1:5046678/150‑1 (MQ=255)
                                                                                                                                  tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg                    <  1:3956552/150‑1 (MQ=255)
                                                                                                                                  tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg                    <  1:3469049/150‑1 (MQ=255)
                                                                                                                                    gCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGc                                                                                                              >  2:2368649/1‑58 (MQ=255)
                                                                                                                                    gCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGc                                                                                                              <  1:2368649/58‑1 (MQ=255)
                                                                                                                                    gCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGc                  >  1:684933/1‑150 (MQ=255)
                                                                                                                                     ccTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGcc                 >  2:1669588/1‑150 (MQ=255)
                                                                                                                                     ccTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGcc                 <  1:2445237/150‑1 (MQ=255)
                                                                                                                                      cTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGACCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCa                >  2:70652/1‑150 (MQ=255)
                                                                                                                                       tCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAg               >  2:1763653/1‑150 (MQ=255)
                                                                                                                                               tGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGaa       <  2:158133/150‑1 (MQ=255)
                                                                                                                                               tGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGaa       <  1:1489356/150‑1 (MQ=255)
                                                                                                                                                 cccACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGAacg     <  1:975298/150‑1 (MQ=255)
                                                                                                                                                  ccACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGAacga    >  2:4275014/1‑150 (MQ=255)
                                                                                                                                                   cacaTCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGAacgac   <  1:4211373/150‑1 (MQ=255)
                                                                                                                                                    acaTCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGAacgacg  >  2:3444684/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGGGCCAGCTGCTGAACGACG  >  NC_002947/541953‑542250

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: