Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,112,209 C→T D449D (GAC→GAT)  opgH → osmoregulated periplasmic glucans biosynthesis protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,112,2090CT94.4% 47.5 / ‑6.0 18D449D (GAC→GAT) opgHosmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/4);  minor base G (1/0);  total (14/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGACTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGT  >  NC_000913/1111967‑1112444
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
gCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTgg                                                                                                                                                                                                                                      >  1:65572/1‑250 (MQ=255)
  cGGGATCATTAAGTTGTCGCCGAAAGCGTCCGGTATGGATACGCGGTATGCGAGCGGTCAGCAGTTCGCGATCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAAACGTTTATCGAGCACTGCGCACTGGCTCCGCTGGCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAg                                                                                                                                                                                                                                   <  1:125857/251‑1 (MQ=255)
                                cgttgtggttACGCTGTATGCGCGCTGTAAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTAGCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTgggggg                                                                                                                                                                                                     <  2:145587/242‑1 (MQ=255)
                                 ggTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGGTggggggg                                                                                                                                                                                                    >  1:372503/1‑250 (MQ=255)
                                                tATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACg                                                                                                                                                                                      >  1:161925/1‑250 (MQ=255)
                                                 aTGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGAt                                                                                                                                                                                    >  1:275230/1‑251 (MQ=255)
                                                                               tgtATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCt                                                                                                                                                      >  1:132063/1‑251 (MQ=255)
                                                                                gtATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCTCACTACTGGGGACATAACGCGATTATCCGCGGGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTa                                                                                                                                                     <  1:345557/251‑1 (MQ=255)
                                                                                                      ggtgtGCACTTCTGGCATCTTGTCGATTCGCACTGCTGGGGACATAACTCGATTATCCGCGTGAAAGCGTTTCTCTATCACTGCGCGCTGGCTCCGCTGCCGTGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTTGATTGCTTACGATCGCCCGGGTTCTTATGAAGTATTTCCGCCTAACTTGCTTGATGAGCTAAAACg                                                                                                                               <  2:275081/247‑1 (MQ=255)
                                                                                                                      aacgtGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGGTGCGCCGTGCAGGTTGGGGGGTCTGGATTTCTTACGATCTCCCGGGGTCTTATGAAGAATTGCCGCCTAACGTGCTTGATGAGCTAAACCATGAccgccgcaggtgc                                                                                                               >  2:301245/5‑244 (MQ=255)
                                                                                                                                     actataggagACATAACGCGAGTATCCGCGTGAAACCGGTTATCGCGCACCGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCAGTTCACTCCTGACACATGAGTTCCTGGAAGCGGCATTGGTGCACCGTGCAGGTTGGGGGGTCTGGGTGGGTAACGATCTCCAGGGTTCTTATTGAGAATTGCGGCTTAACTTGCTTGATGCCCTAAAGGGAGACCGCCGATGGTGCCGCGGTAAcgagatg                                                                                                >  2:107727/10‑245 (MQ=255)
                                                                                                                                             gtaCATAACGCGATTATCAGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCTGGCGAAGGTGCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGACGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGATGACTTGCCGCCTATCTTGCTTGATGCACTAAAACGTGACCGTCGATGGTGACACGGTAACCTGAGGAACTTcgg                                                                                        >  2:59543/3‑249 (MQ=255)
                                                                                                                                                                                cGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCAcc                                                     >  1:215481/1‑251 (MQ=255)
                                                                                                                                                                                         cgcACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGaa                                                                         >  2:225077/1‑222 (MQ=255)
                                                                                                                                                                                                                               gTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTTTTATCTCTCCGCTCCGCtg      >  1:301844/1‑251 (MQ=255)
                                                                                                                                                                                                                               gTTCAATCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGACGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGAGGTGTTCCTGACGGGCGGGATGTTTAATCTCTccgctccgct       >  2:266361/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaTCCTGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTTTTATCTCTCCGCTCCGCTGTGGt  >  1:347432/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tGTCACATGATTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTTTTATCTctccg              >  1:109500/1‑234 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
GCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCGGTTCAATCCTGTCACATGACTTCGTGGAAGCGGCGTTGATGCGCCGTGCAGGTTGGGGGGTCTGGATTGCTTACGATCTCCCGGGTTCTTATGAAGAATTGCCGCCTAACTTGCTTGATGAGCTAAAACGTGACCGCCGATGGTGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGT  >  NC_000913/1111967‑1112444

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