Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,549,816 A→G T806A (ACG→GCG)  fdnG → formate dehydrogenase N subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,549,8160AG76.5% 20.7 / 5.0 17T806A (ACG→GCG) fdnGformate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (9/4);  total (12/5)
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.42e-01

GCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGACGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGTTTGCCATCAACAAAATGGCGGAAGGTCCGTTCCC  >  NC_000913/1549702‑1549940
                                                                                                                  |                                                                                                                            
gCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTg                                                                                                            >  2:444681/1‑133 (MQ=255)
          ggCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTATCGATATTCCTGACTTCGGCaa                                                                                                  >  2:21116/1‑133 (MQ=255)
               cTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTGCTGACTTCGACAATGCCg                                                                                             >  1:217134/1‑133 (MQ=255)
                         gcgTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTAc                                                                                   <  2:377749/133‑1 (MQ=255)
                          cgTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGCACg                                                                                  >  2:281355/1‑133 (MQ=255)
                                 tACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCg                                                                           >  2:450395/1‑133 (MQ=255)
                                    aaCCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGc                                                                        <  1:444681/133‑1 (MQ=255)
                                         tGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGttt                                                                   >  1:396990/1‑133 (MQ=255)
                                         tGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGttt                                                                   >  2:20016/1‑133 (MQ=255)
                                         tGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGACGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCggg                                                                         >  2:516004/1‑127 (MQ=255)
                                         tGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGACGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCggg                                                                         <  1:516004/127‑1 (MQ=255)
                                                   gATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCGCCGGGTACGCCAACCGGGCCGTTTATCATGCAGc                                                         <  1:72943/133‑1 (MQ=255)
                                                                        gATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGt                                    >  2:308805/1‑133 (MQ=255)
                                                                             gAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGGCGGGTATCGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGTTTGcc                               <  1:21116/133‑1 (MQ=255)
                                                                                          aTCCAGTGGAACGGCAGCAAGTGGACGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGTTTGCCATCAACAAAATgg                  >  2:406561/1‑133 (MQ=255)
                                                                                                  gAACGGCAGCAAGTGGACGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGTTTGGCATGAACAAAATGACGGAAGGt          >  1:215372/1‑133 (MQ=255)
                                                                                                          gcaAGTGGGCGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGTTTGCCATCAACAAAATGGCGGAAGGTCCGTTccc  >  2:385688/1‑133 (MQ=255)
                                                                                                                  |                                                                                                                            
GCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATATCAACGGTAAACCGTGGGATCCGAAACGGATGCTGATCCAGTGGAACGGCAGCAAGTGGACGGGTAACGATATTCCTGACTTCGGCAATGCCGCACCGGGTACGCCAACCGGGCCGTTTATCATGCAGCCGGAAGGGATGGGACGCCTGTTTGCCATCAACAAAATGGCGGAAGGTCCGTTCCC  >  NC_000913/1549702‑1549940

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