Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,733 C→T 11.9% A657V (GCT→GTT)  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,7330CT11.9% 100.9 / 6.1 42A657V (GCT→GTT) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (22/15);  new base T (2/3);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAG  >  NC_002947/543588‑543878
                                                                                                                                                 |                                                                                                                                                 
tACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTAcc                                                                                                                                               >  1:1043076/1‑150 (MQ=255)
tACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTAcc                                                                                                                                               >  1:410502/1‑150 (MQ=255)
   cgcgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAt                                                                                                                                            <  2:837122/150‑1 (MQ=255)
     cgTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGAtga                                                                                                                                          <  1:221446/150‑1 (MQ=255)
                tGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGaaa                                                                                                                               >  2:773648/1‑150 (MQ=255)
                tGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACAACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTTACTCTACCGATGAAGTGAAGGaaa                                                                                                                               >  2:564045/1‑150 (MQ=255)
                 gAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAAt                                                                                                                              >  1:631784/1‑150 (MQ=255)
                  aaaGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATc                                                                                                                             <  1:982464/150‑1 (MQ=255)
                    agagACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCga                                                                                                                           <  1:677814/150‑1 (MQ=255)
                                  tCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCt                                                                                                             >  1:166548/1‑150 (MQ=255)
                                           aaCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatg                                                                                                                                        >  1:294454/1‑112 (MQ=255)
                                           aaCTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAtgatg                                                                                                                                        <  2:294454/114‑3 (MQ=255)
                                              tGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTTAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTaa                                                                                                 >  1:377571/1‑150 (MQ=255)
                                              tGATGTACACCGGTTTCGAATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAt                                                                                                                                            <  2:1001110/107‑1 (MQ=255)
                                              tGATGTACACCGGTTTCGAATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAGCGTTACCGAt                                                                                                                                            >  1:1001110/1‑107 (MQ=255)
                                                     cacCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAggg                                                                                          >  2:228555/1‑150 (MQ=255)
                                                          accacGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCgaga                                                                                     >  2:806270/5‑150 (MQ=255)
                                                               gCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTac                                                                                >  2:776849/1‑150 (MQ=255)
                                                                  tACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTacaac                                                                             >  2:547261/1‑150 (MQ=255)
                                                                     tCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                          >  2:183370/1‑150 (MQ=255)
                                                                     tCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                          <  2:410502/150‑1 (MQ=255)
                                                                     tCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                          >  1:833909/1‑150 (MQ=255)
                                                                     tCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAg                                                                          <  1:675264/150‑1 (MQ=255)
                                                                                  gCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTCACGCTACCGATGAAGTGAAGGAAATCGAGATCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtgt                                                             >  2:268146/1‑150 (MQ=255)
                                                                                   cGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTtgtg                                                            <  1:990841/150‑1 (MQ=255)
                                                                                               tgtTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGaa                                                <  1:773648/150‑1 (MQ=255)
                                                                                                 ttAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAacg                                              >  2:386981/1‑150 (MQ=255)
                                                                                                             ttATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGGAGCACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCa                                  >  1:910946/1‑150 (MQ=255)
                                                                                                             ttATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCa                                  >  2:14340/1‑150 (MQ=255)
                                                                                                              tATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCaa                                 <  2:1022080/150‑1 (MQ=255)
                                                                                                                   cGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCga                            <  2:166548/150‑1 (MQ=255)
                                                                                                                      aCGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCgatga                         <  1:35311/150‑1 (MQ=255)
                                                                                                                           aaaGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGATACCGAGATCCAGTACGTCGCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGcc                    >  2:38724/1‑150 (MQ=255)
                                                                                                                            aaGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGAAggacc                   >  1:16479/1‑147 (MQ=255)
                                                                                                                              gCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAAc                 <  1:183370/150‑1 (MQ=255)
                                                                                                                                ccGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACct               >  1:939737/1‑150 (MQ=255)
                                                                                                                                  gCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACctct             <  1:510566/150‑1 (MQ=255)
                                                                                                                                       atcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                         <  1:855206/117‑1 (MQ=255)
                                                                                                                                       atcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGaa                                         >  2:855206/1‑117 (MQ=255)
                                                                                                                                       atcGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGaaa        >  2:1114448/1‑150 (MQ=255)
                                                                                                                                            tAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGCGCGATAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGaa   <  2:314448/150‑1 (MQ=255)
                                                                                                                                             aaCGCTACCGATGCAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAg  <  1:228555/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TACCGCGTGGTTGGTCTGAAAGAGACCGTTATCTTCGCCGACCAACTGATGTACACCGGTTTCGCATACTCGACCATTTCCGGCGTTTCCATCGGTGTTAACGACTTCGTTATCCCGGACGAGAAAGCCCGCATCATCGGTAACGCTACCGATGAAGTGAAGGAAATCGAGAGCCAGTACGCCTCCGGCCTGGTAACCCAGGGCGAGAAGTACAACAAGGTCATCGACTTGTGGTCGAAGGCGAACGACGAAGTGTCCAAGGCGATGATGGCCAACCTCTCGAAAGAGAAG  >  NC_002947/543588‑543878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑