Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,006,613 T→G 18.3% intergenic (+25/+74) ydcJ → / ← PP_5261 metalloprotein/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,006,6130TG18.3% 33.1 / 6.4 22intergenic (+25/+74)ydcJ/PP_5261metalloprotein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (9/9);  new base G (2/2);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

ACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGC  >  NC_002947/6006467‑6006762
                                                                                                                                                  |                                                                                                                                                     
aCCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGC                                                                                                                                                    >  1:128692‑M1/1‑150 (MQ=255)
         cGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGGGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCGCTTCGCGGGCTCGCCCGCTCCCA                                                                                                                                           <  2:267470‑M1/150‑1 (MQ=255)
            tCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAG                                                                                                                                        <  2:794442‑M1/150‑1 (MQ=255)
                          gggTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             >  1:237845‑M1/1‑128 (MQ=255)
                          gggTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                             <  2:237845‑M1/131‑4 (MQ=255)
                                  ggTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgg                                                                                                                  >  2:804188‑M2/1‑142 (MQ=255)
                                      cAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgg                                                                                                                  <  2:623977‑M1/146‑21 (MQ=255)
                                      cAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgg                                                                                                                  >  1:623977‑M1/1‑126 (MQ=255)
                                        ggATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgg                                                                                                            >  1:684890‑M2/1‑136 (MQ=255)
                                            gAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtgggagcgggcga                                                                                                        <  1:303379‑M1/150‑31 (MQ=255)
                                                    gcTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                            >  1:389574‑M1/1‑102 (MQ=255)
                                                    gcTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                            <  2:389574‑M1/106‑5 (MQ=255)
                                                          tGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACA                                                                                                                                         >  1:390762‑M1/1‑103 (MQ=255)
                                                          tGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACA                                                                                                                                         <  2:390762‑M1/103‑1 (MQ=255)
                                                               aaaGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAgg                                                                                                                                 >  2:198808/1‑106 (MQ=255)
                                                               aaaGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAgg                                                                                                                                 <  1:198808/106‑1 (MQ=255)
                                                                                        cTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGGGCCCGCGAAAAGCCCCGTGCAGGTGATCAAAGACGCGGCAGGAt                                                            >  2:472660/1‑150 (MQ=255)
                                                                                         tgtgCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCC                                                                                                                                                   >  1:383587‑M1/1‑62 (MQ=255)
                                                                                         tgtgCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCC                                                                                                                                                   <  2:383587‑M1/62‑1 (MQ=255)
                                                                                              acaGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtg                                                                                                                   <  2:289490‑M1/89‑20 (MQ=255)
                                                                                              acaGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTcgtatgccctgtacctgtg                                                                                                                   >  1:289490‑M1/1‑70 (MQ=255)
                                                                                                                                      ctctTCGCGGGCTCGCCCGCTCCCACAGGTACAGGCCTCGCTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAAc              <  1:472660/150‑1 (MQ=38)
                                                                                                                                                  tCGCCCGCTCCCACAGGTACAGGGCATACGacctgtgggagcgggcgagcccgcgaagaggccggtacagatttacatcgaacccaggttcaacgcctgtgcacaggcctgcagcgaccggcgctcgctttgcgcatacagcgccagctc  <  1:303379‑M2/150‑121 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACCGCGAAGCGTTCGAAGCGGCGCTGGGGTTGCAGGTACAGGATGAGCTGGCGCTGTATGCGCAAAGCGAGCGCCGGTCGCTGCAGGCCTGTGCACAGGCGTTGAACCTGGGTTCGATGTAAATCTGTACCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGTACAGGGCATACGACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCGGTGCAGGTGATCAAAGACGCGGCAGGATTTCCACATCCGGCAAATCCATTGCCGCCACCTCGCCTTCCAGGAACGCACTCAACCGC  >  NC_002947/6006467‑6006762

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