Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,475,495 G→C 23.6% G406G (GGC→GGG PP_3087 ← excinuclease ABC subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,475,4950GC23.6% 21.2 / 6.2 17G406G (GGC→GGGPP_3087excinuclease ABC subunit A
Reads supporting (aligned to +/- strand):  ref base G (11/2);  new base C (2/2);  total (13/4)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TTCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAAGCA  >  NC_002947/3475351‑3475640
                                                                                                                                                |                                                                                                                                                 
ttCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccgg                                                                                                                                               >  3:155822/1‑149 (MQ=255)
      ggCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGCTTGCACCCGGCCGAca                                                                                                                                         <  6:51573/149‑2 (MQ=255)
       gCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACgg                                                                                                                                        >  1:124224/1‑149 (MQ=255)
         aGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCt                                                                                                                                      >  5:132907/1‑149 (MQ=255)
                 cTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACCGTCGGCCGGCTTGCACCCGGCCGAcagtga                                                                                                                                     >  8:73842/1‑137 (MQ=255)
                 cTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACCGTCGGCCGGCTTGCACCCGGCCGAcagtga                                                                                                                                     <  7:73842/142‑6 (MQ=255)
                         accaGCCAGTCGGCGCGGCGCATGGTGTCCAGATCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCtcgtc                                                                                                                                  >  1:33256/1‑137 (MQ=255)
                             gccagTCGGCGCGGCGCATGCTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGCGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCCACGGCACGTCGAGTACGTACaccccc                                                                                                                  >  6:6325/1‑149 (MQ=255)
                                                                                             tttCAGGCGCTGCAGGGCTACAACCAGCGATTCACTGCCGGCCGGTTGCACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACaa                                                    <  1:166763/147‑1 (MQ=255)
                                                                                                 aGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCgg                                              >  3:170296/1‑149 (MQ=255)
                                                                                                   gcgcTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCgggg                                            >  4:186586/1‑149 (MQ=255)
                                                                                                               ttCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTc                                  >  5:186547/1‑147 (MQ=255)
                                                                                                                             cACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGACGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccg                  >  5:274328/1‑149 (MQ=255)
                                                                                                                                tGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGc                                                                     >  3:36538/1‑95 (MQ=255)
                                                                                                                                 gTCGGCCGGGTGCAAGCCGGCCGACAGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAg                         >  3:15254/1‑138 (MQ=255)
                                                                                                                                 gTCGGCCGGCTTGCACCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggc              >  6:13632/1‑149 (MQ=255)
                                                                                                                                             cAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAag    >  8:91746/1‑147 (MQ=255)
                                                                                                                                              aaGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagca  <  3:186586/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAAGCA  >  NC_002947/3475351‑3475640

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