Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,429,997 C→T S45F (TCT→TTT)  cycA → D‑serine/alanine/glycine/:H(+)symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,429,9970CT94.1% 33.9 / NA 18S45F (TCT→TTT) cycAD‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/9);  minor base G (0/1);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGT‑TGTTTAT‑GGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTG  >  NC_000913/4429775‑4430223
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                  
gATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGG‑TGTTTAT‑GGGGTTTTGCAAAACGATTAGcagatcgggaga                                                                                                                                                                                                         >  1:118488/1‑239 (MQ=255)
                                               gTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCAGTGGTACGGGGT‑TGTTTAT‑GGGGTGTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGt                                                                                                                                                         <  1:140318/251‑1 (MQ=255)
                                                             acacAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCAGTGGCGGTGCCATTGGTACGGGGT‑TGTTTAT‑GTGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATgg                                                                                                                                           <  1:163033/251‑1 (MQ=255)
                                                                      ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATGGGTACGGGGT‑TGTTTAT‑GGGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGTGCAATGGGGGAATtgc                                                                                                                                  <  1:210742/251‑1 (MQ=255)
                                                                      ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGAGCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCCTTGGCGGTGCCATTGGTACGGGGT‑TGTTTAT‑GGTGTTTGGCAAAACGGTTAGCCTTGCCGGGCCGTCGATCATTTTCGTGTATATGATCATTGGTTTGATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtgc                                                                                                                                  <  1:294926/251‑1 (MQ=255)
                                                                                                     gTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGCGGCATTGGGACGGGGC‑GGTTTATTGGGGT‑TTGAAAAACGANTAgnnnnntt                                                                                                                                                                                                              >  2:191534/1‑137 (MQ=255)
                                                                                                             cGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGT‑TGTTTAT‑GGGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTcn                                                                                           <  1:191534/251‑2 (MQ=255)
                                                                                                                 gCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGG‑TGTTTAT‑GGGGTTTGGCAAAAAGATTAGCCTTGCCGGGACGGCGATCCTTTTCGTTTATATGATTATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTTGTGCTTTTGAATCTTGAATACAAATCTTTTAGTGACTTTGCttt                                                                                       >  1:15889/1‑250 (MQ=255)
                                                                                                                                        aCAGTCGCTACGGCGCACTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGAGCGCTTGGTTAGGGGTGTGTTT‑T‑TGGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGAACATTTTTGTTTTTCTGTTCATTGGGTTTATGCTTTTTTTCGTTGTGCGGGTAATGGGGGGATTGGTGTGTTCGGCTATGGGATACAAATCTTTTTGTGACTTCGGTTACCATTTTCTCGGGGCGTGGGc                                                                  >  1:210378/1‑249 (MQ=255)
                                                                                                                                              gCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGACATTGGTAAGGGGG‑TGTTTTT‑GGGGGTTGGCaaaaag                                                                                                                                                                                                                          >  1:305508/1‑89 (MQ=255)
                                                                                                                                                       cAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTAAGGGGG‑TGGTTAT‑GGGGGTTGGAAAAAAGGTTAGCCTTGCCGGGGCGTCGATCATTTTCGTTTATAAGATTATTGGGTTTATGCTCTTTTTTGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATAAAAATCTTTTAGTGACTTCGCTTCCGATTTAATCGGGGCGGGGGCAGGAAAttttaacgggt                                                 >  1:42011/1‑243 (MQ=255)
                                                                                                                                                                       cATATTCAGCTTATTGCCATTGGCGGTGGCCTTCGTACGGGGT‑TGGTTAT‑GGGGGTTGGCAAAACGATTTGCCTTGCCGGGCCGGCGATCATTTTTGGTTTTATGATCATTGGTTTTa                                                                                                                                                                      >  1:92335/1‑118 (MQ=255)
                                                                                                                                                                          tctgcctatcttGCCAGTGGCGGTGCCAGTGGTACGGGTT‑TGTTTAT‑GTGGTTTGGCAAAACGATTAGCCTTGCCGGGTCGGTGATCATTTTCGTTGATAGAATCATTTGTTTTATGCGCTTTTTCGTGATGCGTGCAATGGGGGAATTGCTGCCTGCGAATCTGGAATTCAAATCTGTGATTGACGTCTCTTCCGAGTTACGCGTGCCTGGGGCAGGATATTTCACCGGCTGGACTTACTGGTTctgtt                               <  1:316128/240‑3 (MQ=255)
                                                                                                                                                                                                    cATTGGTACGGGGT‑TGTTTAT‑GGGGTTTGGCCAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTCTATGATCATTGGTTTTCTGCTCTTTTTCGt                                                                                                                                                         >  1:22910/1‑102 (MQ=255)
                                                                                                                                                                                                       tNNTACGGGNN‑TGTTTAT‑GGGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGt                                                                                                                                                         <  2:22910/99‑1 (MQ=255)
                                                                                                                                                                                                       tGGTACGGGGT‑TGTTTAT‑GGGGTTTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGAAGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGTTGGGTTGTAACCGGTATGGCAGACgtggtg  <  1:323443/250‑1 (MQ=255)
                                                                                                                                                                                                              ggtt‑tgttTAT‑GGGGTTTGGCAAAGCGATTAGCCTGGCCGGGTCGTCGATCATTTTCGTTTATATGATCATTGGGTTTATGCTCTTTTTCTTGATGCGGGCAATGGGGGAATTGCTGCTTTCg                                                                                                                          <  1:55957/120‑1 (MQ=255)
                                                                                                                                                                                                                           ggggTTTGGCAAAACGAGTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTGCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGGGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGTTGGGTTGTAACCGGTATg              <  1:46320/220‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                  
GATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGT‑TGTTTAT‑GGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTG  >  NC_000913/4429775‑4430223

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑