Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,440 C→T T421T (ACC→ACT)  creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,4400CT86.7% 23.5 / NA 17T421T (ACC→ACT) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (4/9);  minor base G (0/2);  total (5/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTT  >  NC_000913/4639201‑4639648
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                
gcgcTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGa                                                                                                                                                                                                       <  1:218435/251‑1 (MQ=255)
                   gTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGAAGGCGGTATTGAAAGGCTTGCGCGACAGCACGTTGTTTAACCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAACGCGCGGGTTTTGGGAACCAGTGTGCTGGGGGTGGCGCTGGCCGGCATGGTGTTTTGGACTCGTAATATCGACTGGTATGGGTTTTCACTg                                                                                                                                                                                    >  1:103121/1‑251 (MQ=255)
                       tacgggtgCATGTATAAACGCTAGTCTGATAGGGGCGTTAATTAACGGTATTCATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGTTTGTTGTTTACCCGCGCGCGTTTGTTGCTGGATGGTTTGATGTGGGGATTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGTTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACGCGTAATATCGACTTGTATGCGTTTTTACTGCCg                                                                                                                                                                                 <  2:343827/243‑1 (MQ=255)
                                          gCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAaaga                                                                                                                                                             <  1:111018/251‑1 (MQ=255)
                                                  taTTGTGGCGTTAATGAACGGTATTTATTTGTAGGCGGTATGGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGTATGGTGTGATGTGGTGTCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATTATGGTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTATAGAAGTTa                                                                                                                                                        <  1:42803/247‑1 (MQ=255)
                                                   agttgggcattAATGACCGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTCCAACAGCATGTTGTTTACCCTCGCGCTGTGGTTGCTGGATGGTGTGATGTGGGGATTGCTCATCTCCTCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTTGCGCTGGCCGGCATGATGTTTGTGACCCGTAATAACGACTGTTATGCGTTTTCCCTGCCGAAAATGAAAGCCAGTAAAGAAGATACaa                                                                                                                                                     <  2:93819/240‑1 (MQ=255)
                                                      ggggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGTCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGTTGTGATGTGGGTACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATAACGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGa                                                                                                                                                 <  1:222128/251‑1 (MQ=255)
                                                       gggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGAc                                                                                                                                                <  1:74858/251‑1 (MQ=255)
                                                                 gcaCGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAAAAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGAGATGTGGGGACTGCACACCTCTGCCGATAGCGCGCTGTTGTTGGGAAACAGTGTGCTGGGGGTGGCGATGGCCGGCATGATGATTGGGACCCGTAATATCGACGGGTGTGCGTTTTAACTTCCGAACTCGAAAGCCAGAAAAGAAGGTAAAACGGaaggttctgtgc                                                                                                                                      >  2:101945/3‑240 (MQ=255)
                                                                                           tGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGGGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCTGCATGATGTTTGTGACTCGTAATATCGACTGGGATGCGTTTTCACTGCCGAAAATGAAAGCCAGGAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAAt                                                                                                             >  1:260180/1‑250 (MQ=255)
                                                                                               aggctggcgCAACTGCTTGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGGTTGTGATTTGGTGACTTCTCAACTCTGCCGATTGCGCGCTTTTGTTGGGAACCAGTGTGCTGGTGGTGGTGCTGGCCGGCATTATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTATAGAATTTACAATGGACGATCAAGTACGTATCTGGAAATAAGGTTGAAAAATaaaaa                                                                                                        <  1:221951/243‑1 (MQ=255)
                                                                                                                                                                    aaaCTCGGCCGATAGCGGGCTGTTGTTGGGTACCAGTGTGCTGGTGGTGGCGCTGGCCGGCAGGATGTTTGTGACTCTTAAAATCGACTGGTATGCTTTTTCCCTCCCGAAAATGAAAGCAAGTAAAAAAGTTACAATGGCCGATCAATTACGTTTCTGGAAATAAGTTTGAAAAATAAAAACGGCGTTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGa                                      <  1:129224/247‑1 (MQ=255)
                                                                                                                                                                           gCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACGCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATaa                              <  1:299231/249‑1 (MQ=255)
                                                                                                                                                                                tAGCGCGCTGTTGTTGGGAACCAGTGAGCTGGTGGTTGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGCTCAATTACGTATCTGGAAATAAGGTTGAAAATTAAAAACGGCGCTAAAAAGAGCCGTTTTTTTTGACGGTGGAAAAGCCGAATAATCTTCCACATCACAACAGAAGCGATAACATTCAc                       >  2:191426/1‑251 (MQ=255)
                                                                                                                                                                                          gttatgGGAACCAGTGTGATGGTGGTGGGGCTGTCCGGCATGATGTTTGTGACTCGTAAATTCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGAtggt             <  2:89709/247‑1 (MQ=255)
                                                                                                                                                                                                ggAACCAGTGTGCTGGAGGTGGCGCCGGCCGGCATGATGTTTGGGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTAAAATGGACGATCAATTACGAATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAACAAGCGCCGTTTTTTTTGTCGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAAACTTCACCGTGGATGGGGGCgat       >  1:308714/1‑251 (MQ=255)
                                                                                                                                                                                                     cAGTGTGCGTGTGGTGGCGCTGGCCGGCATGATGTTTTTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGAttt  <  1:356645/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                
GCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTT  >  NC_000913/4639201‑4639648

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