Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,356,932 G→A T60I (ACC→ATC)  sapA ← putative periplasmic binding protein SapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,9320GA90.9% 52.6 / NA 22T60I (ACC→ATC) sapAputative periplasmic binding protein SapA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/10);  minor base C (2/0);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGC  >  NC_000913/1356725‑1357138
                                                                                                                                                                                                               |                                                                                                                                                                                                              
cGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAgg                                                                                                                                                                      <  1:79480/250‑1 (MQ=255)
    ggAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGgtgt                                                                                                                                                                   <  1:907/249‑1 (MQ=255)
                           ncattcattTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCCGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGAc                                                                                                                                          >  1:146220/2‑251 (MQ=255)
                            cattcattTTACGAGTGGGAGTAAACCAGTCGGATTTTTGAAACGGAACCTCGCGCCGCAGGAGGAAGCGATAGATCGCGCAGTTGTCGAGTACTTCCAAGCTTTCGGCAAGTTCCGGCATCAGGCGCTAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGCTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGCCGCAATAGAc                                                                                                                                          >  1:382764/1‑250 (MQ=255)
                                     tACGAGTGGGCGTAAACCAGTCGGTTTTGTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGGGGTCGAGTACTTCCCGGCTGTCGGCAAGTTCCGGCATCAGGCTATAGGTATAGGGATCGACATCCATCAGTCTATCATAAAACTGGGTGGCAAGGAGATCGACAATTAACCCACGGCTCGCttc                                                                                                                                                                                    <  1:156554/199‑2 (MQ=255)
                                             ngAGTAAACCAGTCGGTTTTTTGAAACTGAACATTGCGCCGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCCGCCTTGTGCAAGGTCCGGCCTCAGGCGCTAGGTATCGGGATCGACCTCCAGCCGGCGATCATAAAACTTGGCGGCAAGGATATCGACAATTAACCCGCTGCTCGCTTTGGCTGGGTTAAAGGGGTTGACTTTCCCGCTGACGCAATAGGCAAAACCGCTGTCGCGGat                                                                                                                        >  1:90013/2‑251 (MQ=255)
                                                                                  naNGCNNGTGGNAGCGATAGGTCGCGCCGTTGTCGANTACTTCCAAGCGTTCGGCAAGTTCCGGCATCNGNCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGNCAAGGATATCNACANNTAACNCNCNGCTCGCTTTGGATGGGTTAAAGGTGTTGANTTGCCCGCTGACGCAATAGACAANNCCNNNGTNNCGGat                                                                                                                        <  2:382764/213‑1 (MQ=255)
                                                                                  gACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGcc                                                                                   <  1:1280/251‑1 (MQ=255)
                                                                                                                 gTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCTGCAAGGATATCg                                                                                                                                                                                                           <  1:411213/100‑1 (MQ=255)
                                                                                                                         cTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATACGTATAGGGATAGACATCCAGCAGTAGCTCATCAAACTGGGCGGAAAGGCTATCGACCATTAACCCAGTGCTCGCTTTGGATGGGt                                                                                                                                                                            >  1:478887/1‑123 (MQ=255)
                                                                                                                           tCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCg                                                                                                                                                                                        <  1:359984/109‑1 (MQ=255)
                                                                                                                           tCCCAGCTNTNNNCNNGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTNGATCATAAAACTGGGNGGNAAGGNTATNGACAATTAACCCACTGCTCg                                                                                                                                                                                        >  2:359984/1‑109 (MQ=255)
                                                                                                                                 cggtctgCAAGTGCCGGCATCATGCTATAGGTAGAGGGATCGGCATCCATCAGGCGATCATAAAACTGGGCGGCAAGGATATCGTCAATTAACTCACAGCGCGCTGGGGATggg                                                                                                                                                                             <  1:410786/108‑1 (MQ=255)
                                                                                                                                     nggCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGgcgc                                >  1:131935/2‑251 (MQ=255)
                                                                                                                                      ggCAAGTTCCGGCATCAGGCGATCGGTATAGGGCCCGAAATGCAGCAGTCCATCATAATACTGGGCGGCCAGGATATCGACCATTAACCCCCTGCTCGCTTTGGACGGGTTAACGGTGTTGACTTGACCGCTGACGCAACAGCCAAAACCGCTGTCGCGGATATCAGCATGc                                                                                                              >  1:436378/1‑172 (MQ=255)
                                                                                                                                                ggCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGGCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGGCGCAATAGACAAAACCGCTGTCGCgg                                                                                                                          >  1:376213/1‑150 (MQ=255)
                                                                                                                                                                tataGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACa                                                                                                                                         >  1:219399/1‑119 (MQ=255)
                                                                                                                                                                   aGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTTACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAATTCCAGCAATCACCAAAATAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGc  <  1:434264/251‑1 (MQ=255)
                                                                                                                                                                        tCGACATCCAGCAGTCGATCATAAAACAGGGAGGCAAGGATATCGACAATTAACCCAATGCTCGATTTGGATGGGTTAACGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCgcggc                                                                                          >  1:155635/1‑158 (MQ=255)
                                                                                                                                                                                        gATCATAAATCTGGACGGCAGGGATANCAACATNNAANACACTGCTCGATTTGCATGNGTNNNANNNGNTGAAANGCNNGCTGACGCANTAGACAAAACCGCTGTCGCGGATATCAGCATGAGGGCGAGATTCCGGCGCGGCNATTGCAAGANNACTCCCAAGTCCAGCAATCACNCAANNANNCNatantan                                       >  2:112041/1‑192 (MQ=255)
                                                                                                                                                                                        gATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGcc                                                                                   >  1:26882/1‑149 (MQ=255)
                                                                                                                                                                                                    gggCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGAGGTGATAAAGGTGTTGACTTGCCCTCTGACGCAATAGACAAAACCGCTGGCGCGGATAGCAGCATGCGGTGGAGATTCAGGCGCGGCGATTGCCTGaac                                                                              <  1:484529/142‑3 (MQ=255)
                                                                                                                                                                                                          cAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTcaca                                                                       >  1:338745/1‑143 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                              
CGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGC  >  NC_000913/1356725‑1357138

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