Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,430,615 C→T S251F (TCC→TTC)  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,430,6150CT95.0% 56.1 / NA 20S251F (TCC→TTC) cycAD‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/8);  minor base G (0/1);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTG  >  NC_000913/4430371‑4430853
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
gtcgCCATTGTGGCGCTGATTGTCTTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGGCATTCGCGCATTTGTGGAATGACGTCGGCTGGTTCCCGAAAGGTTTAATTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTGCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTGCATTcc                                                                                                                                                                                                                                          <  1:247354/251‑1 (MQ=255)
                             ccTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGGCATTCGCGCATTTGTTGAATGACGGCGGCTGGTTCCCGATAGGTTTAAGTGGCTTCTTTGCCGTATTCCATATAGCGGTTTTCGCTTTCGTTGGGATTGAGCTTGTAGGTACAACAGCTGCGTAAACCAAATATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGt                                                                                                                                                                                                                         <  1:303246/239‑1 (MQ=255)
                                             aTGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCGCTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTcg                                                                                                                                                                                                             <  1:46512/235‑1 (MQ=255)
                                                                                                     cggcgcatGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTAGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAAACAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGt                                                                                                                                     >  1:75072/8‑251 (MQ=255)
                                                                                                       ncggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGt                                                                                                                                     >  1:65531/2‑249 (MQ=255)
                                                                                                                  ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGTTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTGCGCGCGGAGTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGc                                                                                                                                                     <  1:167302/222‑1 (MQ=255)
                                                                                                                                  gCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAACTCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGc                                                                                                         >  1:211390/1‑250 (MQ=255)
                                                                                                                                                  ncAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTcc                                                                                                                                                                >  1:192372/2‑179 (MQ=255)
                                                                                                                                                   cAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCACGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTAGCGCTGATTGTGTTTACGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTtgt                                                                                       >  1:316362/1‑251 (MQ=255)
                                                                                                                                                         gaggTTTTCGCTTTCGTGGGGATTGANCTGGTANNTANAACAGCTGCGGAAACCAAANATNNNGNNNANTCACTNCCNNGCGCGATTANCTTCATTCCGATCCGTCTCATTATGTTCTACGTCTTCGCGCTGATTGTGATTANGTCCGTGACNNCGTGGAGTTCGGTAGTCCCGGNGAANNGNNCNTttntt                                                                                                                                            >  2:69802/3‑192 (MQ=255)
                                                                                                                                                                           nnnnnnnnnntGNNNNNTACNACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTcc                                                                                                                                                                <  2:192372/144‑1 (MQ=255)
                                                                                                                                                                                             acaacaGCTGCGGAAACCAAAGATCCAGAGAAATCACTGTCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGc                                                                                                                                                     >  1:25164/1‑147 (MQ=255)
                                                                                                                                                                                                  aGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCAACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGGGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTctct                                         >  1:256150/1‑250 (MQ=255)
                                                                                                                                                                                                        ggAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGCTTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGACAAGCCCGTTTGTTGAATTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGcc                                  >  1:316231/1‑251 (MQ=255)
                                                                                                                                                                                                        ggAAACCAAAGATCCAGAGAAATCACCGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCAGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGcc                                  >  1:142856/1‑251 (MQ=255)
                                                                                                                                                                                                              cAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATNc                            <  1:255153/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcgATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGtgtg  <  1:43137/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gcgATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGTCGTCTTCTCTACCAGCCGTATTCTGTGTGGTCt                  <  1:259229/233‑1 (MQ=255)
                                                                                                                                                                                                                                           cgATTAACTTCNTTCCGATCCGTATCATTATGTTCTACGtn                                                                                                                                                                                                                 >  2:43388/1‑40 (MQ=255)
                                                                                                                                                                                                                                           cgATTAACTTCATTCCGATCCGTATCATTATGTTCTACGTc                                                                                                                                                                                                                 <  1:43388/41‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
GTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTG  >  NC_000913/4430371‑4430853

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