Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,125 T→G 19.3% T605P (ACT→CCT)  rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,1250TG19.3% 58.7 / 8.5 31T605P (ACT→CCT) rpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base T (14/11);  new base G (3/3);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGGTGTTCATGTCGATACCGAAACGCATGCGCAGCACCTTGG  >  NC_002947/470980‑471262
                                                                                                                                                 |                                                                                                                                         
gtCTGAAGTCGAGTGTAGACGGGCATTGCCTCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCggg                                                                                                                                       >  2:350808/1‑150 (MQ=255)
             tgtAGACGGGCATTGCCTCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgc                                                                                                                          <  1:350808/150‑1 (MQ=255)
                    gggCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGTAAAACCCCGGCCTGGGCCGGGGGTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGc                                                                                                                   >  2:104114/1‑150 (MQ=255)
                              gCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGcc                                                                                                                  >  2:253820/1‑141 (MQ=255)
                              gCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGcc                                                                                                                  <  1:253820/141‑1 (MQ=255)
                                          cGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGGTGGCGCAGCTTGCGCAAc                                                                                                                     <  2:36315/126‑1 (MQ=255)
                                          cGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGGTGGCGCAGCTTGCGCAAc                                                                                                                     >  1:36315/1‑126 (MQ=255)
                                                  gggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGa                                                                                                   <  1:130569/136‑1 (MQ=255)
                                                  gggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGa                                                                                                   >  2:130569/1‑136 (MQ=255)
                                                            tCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTATTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCGGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAAGCGTTACGTCa                                                                           >  1:583722/1‑150 (MQ=255)
                                                                   ccACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttc                                                                                                               <  2:212179/107‑1 (MQ=255)
                                                                   ccACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttc                                                                                                               >  1:212179/1‑107 (MQ=255)
                                                                     aCTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTg                                                                  <  1:330699/150‑1 (MQ=255)
                                                                        atgacaaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGccc                                                               >  1:532346/6‑150 (MQ=255)
                                                                         aaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCa                                                              >  1:96588/1‑150 (MQ=255)
                                                                            aaaaaCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACct                                                           >  2:371632/1‑150 (MQ=255)
                                                                             aaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACctc                                                          >  2:546428/1‑150 (MQ=255)
                                                                             aaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACctc                                                          >  2:527896/1‑150 (MQ=255)
                                                                               aaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCtt                                                        >  2:48536/1‑150 (MQ=255)
                                                                                aCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc                                                       <  1:104114/150‑1 (MQ=255)
                                                                                aCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTc                                                       <  1:305798/150‑1 (MQ=255)
                                                                                  cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa                                                     >  1:160970/1‑150 (MQ=255)
                                                                                  cccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa                                                     >  1:131833/1‑150 (MQ=255)
                                                                                  cccGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGa                                                     <  2:532346/150‑1 (MQ=255)
                                                                                    cGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg                                                   <  2:160970/150‑1 (MQ=255)
                                                                                              cGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCgg                                         >  2:66044/1‑150 (MQ=255)
                                                                                                gggTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGGCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGgtg                                       <  1:527896/150‑1 (MQ=255)
                                                                                                                           gaGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGGTGTTCATGTCGATACCGAAACGCATgcgc            <  1:48536/150‑1 (MQ=255)
                                                                                                                                   gTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGGTGTTCATGTCGATACCGAAACGCATGCGCAGCACCtt    >  2:453355/1‑150 (MQ=255)
                                                                                                                                    tGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGGTGTTCATGTCGATACCGAAACGCATGCGCAGCACCTTg   <  1:546428/150‑1 (MQ=255)
                                                                                                                                     gCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGGTGTTCATGTCGATACCGAAACGCATGCGCAGCACCTTgg  <  1:371632/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
GTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGGGTGTGGTCGGTGTTCATGTCGATACCGAAACGCATGCGCAGCACCTTGG  >  NC_002947/470980‑471262

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: