Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 462,097 C→T S61S (AGC→AGT ppiD → periplasmic folding chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913462,0970CT93.8% 32.6 / ‑3.3 16S61S (AGC→AGTppiDperiplasmic folding chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (10/5);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

GTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCAAATTTGATAACAGCCG  >  NC_000913/461870‑462324
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   
gTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAAc                                                                                                                                                                                                               <  1:248462/250‑1 (MQ=255)
                                                            cGCACGCCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTc                                                                                                                                                   >  1:140197/1‑250 (MQ=255)
                                                                      naaaCAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTcc                                                                                                                                                                                           >  1:306458/2‑200 (MQ=255)
                                                                                   gCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAAATACGCCGCAAAAGTGAATGACCAGGCAATAAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTc                                                                                                                                 >  1:338268/1‑245 (MQ=255)
                                                                                                                    nnnnnnnnnnaCNNNNNTGANTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTcc                                                                                                                                                                                           <  2:306458/144‑1 (MQ=255)
                                                                                                                                    gTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCtat                                                                                                                                                                   >  1:25984/1‑162 (MQ=255)
                                                                                                                                                             gCAATAACTACGCCGCAAAAGTGAATNACCAGGNNATNAGCCGTGGGCAATTCGAGANCGNNNTNNNCNGTGAGNGTNNTCGCATGCANCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGCANACCCTGCGTNNACAGGTGCTGAATCGTCTGATNGACNNGNNGNTGCNggn                                                                                                           >  2:113908/1‑192 (MQ=255)
                                                                                                                                                                   aCTACGCCGCAAAAGTGAATGACCAGTAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTcgcg                                           <  1:420197/251‑1 (MQ=255)
                                                                                                                                                                       cgccgcAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGAccc                                       >  1:13386/1‑251 (MQ=255)
                                                                                                                                                                                 gTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCAGGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCa                                                                                                                                                  <  1:237634/134‑1 (MQ=255)
                                                                                                                                                                                        aCCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTAAACagtgag                                                                                                                                                                                                                                  >  1:225889/1‑47 (MQ=255)
                                                                                                                                                                                            ggAAATCAGCCGTGGGCAATTCGAGAAAGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGc                   >  1:428739/1‑250 (MQ=255)
                                                                                                                                                                                                          ggCAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACCGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGACGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACCGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCAAATTTGATAACAGc    >  1:446883/1‑251 (MQ=255)
                                                                                                                                                                                                           gCAATTCGNNACNNCNTNCAACANTGAGCGTAATCGCATGCAGCAACAGCTggg                                                                                                                                                                                                        >  2:342657/1‑54 (MQ=255)
                                                                                                                                                                                                           gCAATTCGAGAACGCCTGCAACAGTGAGCGTAATCGCATGCAGCAACAGCTggg                                                                                                                                                                                                        <  1:342657/54‑1 (MQ=255)
                                                                                                                                                                                                            cAATTCGAGAACGCCTTCAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCTGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCAAATTTGATAACAGCCg  >  1:504101/1‑251 (MQ=255)
                                                                                                                                                                                                                             cAACAGTGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATcgacga                                                                                                                       <  1:430844/117‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   
GTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCAAATTTGATAACAGCCG  >  NC_000913/461870‑462324

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