Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 807,020 C→T G88S (GGC→AGC)  ybhC ← outer membrane lipoprotein YbhC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913807,0200CT88.9% 17.1 / ‑5.4 9G88S (GGC→AGC) ybhCouter membrane lipoprotein YbhC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/3);  minor base G (0/1);  total (5/4)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GGTTGACGTCGTGACGCCAGTCGGCAGGGCTCATGCCACCATCAAGGGAAAGCCCAATCTTCACATCAATCGGTTTTTCACCTGTACCGTACAGAGTAATTCCACCCGGAGCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAG  >  NC_000913/806780‑807227
                                                                                                                                                                                                                                                |                                                                                                                                                                                                               
ggTTGACGTCGTGACGCCAGTCGGCAGGGCTCATGCCACCATCAAGGGAAAGCCCAATCTTCACATCAATCGGTTTTTCACATGTACCGTACAGAGTAATTACACACGGAGCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCTCGCCGGGcc                                                                                                                                                                                                        >  1:194754/1‑250 (MQ=255)
                       gCAGGGCTCATGCCACCATCAAGGGAAAGCCCAATCTTCACATCAATCGGTTTTTCACCTGTACCGTACAGAGTAATTCCACCCGGAGCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCGCCTGAATCGTGGGATGCGTTACACCTTGAGTGCTCGCCGGGCCGACAACAAAGTCAGGTTGCgcagg                                                                                                                                                                                >  1:50439/1‑251 (MQ=255)
                                   ccaccaTCAAGGGAAAGCCCAATCTTCACATCAATCGGTTTTTCACCTGTACCGTACAGAGTAATTCCACCCGGAGCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTTAGTGCgcgccg                                                                                                                                                                                                            <  1:174207/211‑1 (MQ=255)
                                                                                                     ccACCCGGAGCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTAAGCTTGATAATTGCCGCATCTAACGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCTCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTACACGCTGCAGCACCt                                                                                                                                           >  1:199680/1‑210 (MQ=255)
                                                                                                              gCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGAATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCCTCTACCCCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCTCGCCGGGCCGACAACCAAGTCAGGTTTCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCttt                                                                                         >  1:203415/1‑251 (MQ=255)
                                                                                                                      gcnatancnngNNCCCNGATACTCACCAGGCATCACGGNNATATACTGGNGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCGCCTGAATCGTGGTANGCGTTACACNNTGNGTGCGNGNNNGNNNGANAACAAAGTCAGGTTGCGCANGCNNGGTAATNGGGGAAGGATTCCACGCTGCAGCAcc                                                                                                                                            <  2:50439/182‑1 (MQ=255)
                                                                                                                          tataCCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCTCCTGAATCGTTGGATGCGGTACACCTTGAGTGCTCGCCGGGCCGACAACAAAGTCAGGTTGCGcaggcaggn                                                                                                                                                                           <  1:21435/157‑2 (MQ=255)
                                                                                                                                               ccAGGCATCACGGCAATATACAGGCGCTTGTTGGTACGCTTTATAATTGCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCTCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCATGGTAATCGGGGAAGGATTCCTCGCTGCAGCACCTGTTGTCAGGGATTCAAAATAGTGTCGAGCATCTAAATTCTGCGCTGCTTTTGCCGACAGAATCGGGCGAGAAGAGGTCCCAgg                                                        <  1:200190/251‑1 (MQ=255)
                                                                                                                                                                                               gCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCTCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGc        >  1:181092/1‑251 (MQ=255)
                                                                                                                                                                                                        aCCGCCGCCTGAAGCGTGGTATGCGTTACTCCTTGAGTGCTCGCCGCGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCTCGCTGCAGCACCTGGTGTCATGTATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCTGTTTGATCAGAAGGTGGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAg  <  1:92287/248‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                               
GGTTGACGTCGTGACGCCAGTCGGCAGGGCTCATGCCACCATCAAGGGAAAGCCCAATCTTCACATCAATCGGTTTTTCACCTGTACCGTACAGAGTAATTCCACCCGGAGCGGCAGGGACATATACCGTTCCCTGATACTCACCAGGCATCACGGCAATATACTGGCGCTTGTTGGTACGCTTGATAATTGCCGCATCTACCGCCGCCTGAATCGTGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAG  >  NC_000913/806780‑807227

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