Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,994,670 G→T 19.1% P213H (CCC→CAC)  PP_2618 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,994,6700GT19.1% 37.6 / 2.1 20P213H (CCC→CAC) PP_2618hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/9);  new base T (2/2);  total (9/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 = 9.75e-01

GATGTGCGTGTCGTTGACCAACACCGACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTGCGCAGCGCCTTGGCCGGGTCCTGCTCCAGTAGCTGGTACACGGTGTTGCGCAGCAGTGAGTCGTCGCGCAGTTGC  >  NC_002947/2994531‑2994813
                                                                                                                                           |                                                                                                                                               
gATGTGCGTGTCGTTGACCAACACCGACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACgg                                                                                                                                       <  1:232869/150‑1 (MQ=255)
        tgtCGTTGACCAACACCGACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACt                                                                                                                               >  1:121800/1‑150 (MQ=255)
         gtCGGTGCCCAACACGAACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCACTACTTCCCGCCAACTGGCGTTTTGATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACtt                                                                                                                              <  2:127997/150‑1 (MQ=38)
                 ccAACACCGACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGc                                                                                                                      <  1:515706/150‑1 (MQ=255)
                       ccGACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGATTCGGCGTGGATGGGTTGGAACGGGATGGACTTGAATTGGCCAGGgc                                                                                                                >  2:307078/1‑150 (MQ=255)
                          aCAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCgca                                                                                                             >  1:286684/1‑150 (MQ=255)
                              cATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATgg                                                                                                                                                 >  2:321615/1‑110 (MQ=255)
                              cATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATgg                                                                                                                                                 <  1:321615/110‑1 (MQ=255)
                               aTGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGCATGGGTTCGAACGGGATGTAATTTAATTGGCCAGGGAGCAGCAgg                                                                                                        >  2:16372/1‑150 (MQ=255)
                                 gATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTa                                                                                                      <  1:63931/150‑1 (MQ=255)
                                                 ggcggAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGata                                                                                      <  2:286684/150‑1 (MQ=255)
                                                     gAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAgg                                                                                  <  1:244261/150‑1 (MQ=255)
                                                       aTGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGTGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAg                                                                                   <  2:537857/147‑1 (MQ=255)
                                                       aTGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGTGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAg                                                                                   >  1:537857/1‑147 (MQ=255)
                                                                  ttGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAgtc                                                                                                   <  1:13762/120‑1 (MQ=255)
                                                                  ttGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAgtc                                                                                                   >  2:13762/1‑120 (MQ=255)
                                                                      cgcaAACTCTCGTTTTCATGGCCTCGATAGTCGGCGATCACCCCTACGAAGCGTATTTCTTCTTAGATCTGTTCTAACGTGATGTACTTTAATTCTCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTTCGCAGCGCCTTg                                                                 <  2:219287/150‑1 (MQ=255)
                                                                            cTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGGGTGGATGGGTTCGAACGAAATGTACTTGAATTGGCCCGGGCGGAGCAGTTAGTCGTCACTAAGGATATCGGTAGTGCGCAGCGCGTTTGctcgg                                                           >  2:91568/1‑146 (MQ=25)
                                                                             tGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGATCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTGCGCAGCGCCTTGGCCGGGt                                                          >  2:222779/1‑150 (MQ=255)
                                                                                            cTCGGTCGTCAGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTGCGCAGCGCCTTGGCCGGGTCCTGCTCCAGTAGCt                                           <  1:307078/150‑1 (MQ=255)
                                                                                                                 cgacgaAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTTAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTGCGCAGCGCCTTGGCCGGGTCCTGCTCCAGTAGCTGGTACACGGTGTTGCgcagca                      <  1:16372/150‑1 (MQ=255)
                                                                                                                                     tGGATGTGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTGCGCAGCGCCTTGGCCGGGTCCTGCTCCAGTAGCTGGTACACGGTGTTGCGCAGCAGTGAGTCGTCGCGCAGTTGc  >  1:488450/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GATGTGCGTGTCGTTGACCAACACCGACAGCATGATCTGGCGGCCACGCGGCGGAATGCGCAGTACTTGCCGCCAACTGGCGTTTTCATGGCCTCGGTAGTCGGCGATCACCGCGACGAAGCGGGTTTCGGCGTGGATGGGTTCGAACGGGATGTACTTGAATTGGCCAGGGCGCAGCAGGTAGTCGTCACTACGGATATAGGTACTGCGCAGCGCCTTGGCCGGGTCCTGCTCCAGTAGCTGGTACACGGTGTTGCGCAGCAGTGAGTCGTCGCGCAGTTGC  >  NC_002947/2994531‑2994813

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