Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,188,087 T→C T110A (ACG→GCG)  yehA ← putative fimbrial adhesin YehA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,0870TC93.3% 41.8 / ‑5.9 15T110A (ACG→GCG) yehAputative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (6/8);  minor base G (0/1);  total (6/9)
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.63e-01

AGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAAC  >  NC_000913/2187874‑2188315
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
aGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGTATAATGTGGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCGTCAGGGGTCACCGGAACTGAATAACCACTGGTGAAACTGGTCTGATTAttt                                                                                                                                                                                                 <  1:245891/251‑1 (MQ=255)
aGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTAttt                                                                                                                                                                                                 >  2:201181/1‑251 (MQ=255)
          attttgagagtatgaaTTTCTGTAGATATTTACACTTGCACCCCATTACTTGCAACCTAATAAACACAGTAGAAATTGTCCTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAGCATCATTAATCATCACTGAGCTGCCGACCCCCTCAGAAACATTCACATTGGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGCAACTGGTCTGATTATTTGACAACGTGa                                                                                                                                                                                       <  2:142465/236‑1 (MQ=255)
                  ccAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGttt                                                                                                                                                                               <  1:204014/251‑1 (MQ=255)
                               gTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCa                                                                                                                                                                    <  1:19276/249‑1 (MQ=255)
                                     aTTTAAACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCATCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCaaaa                                                                                                                                                             <  2:132438/250‑1 (MQ=255)
                                                 aTCCCATGACTTGCCACCTAATAAACACAGTAGAACTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCCACTGATTAACCACTGGTTAAACTGGTCTGATTATTTGACAACGTTATGTTGTTTTTAGTGGTCCATTCAAAATGTAaatttaatt                                                                                                                                                >  2:208328/1‑251 (MQ=255)
                                                       tGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGa                                                                                                                                                                                               >  1:176502/1‑198 (MQ=255)
                                                       tGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGa                                                                                                                                                                                               <  2:176502/198‑1 (MQ=255)
                                                                                                           gtCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGt                                                                                      >  2:94000/1‑251 (MQ=255)
                                                                                                                                                       caCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTTTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAAGTTAATTTTATTGTATCGTTATAAAGTCCAACAATAATATCCGACTGGCTCATAGGGGTGATAGTATTGGTACTTACGCATGACGAATCGTTAAAATTATGCTCatat                                           <  2:163714/250‑1 (MQ=255)
                                                                                                                                                                 ccgccCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTtg                                <  1:113872/251‑1 (MQ=255)
                                                                                                                                                                                               ttACTGGCCGCAGGCGTCACCGCAACTGAATCACCACTCGTGAAACTGGTCTGAATATTTGAAAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATAGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGGGCAAAATACCCCTGaac  >  2:241850/1‑251 (MQ=255)
                                                                                                                                                                                                          aGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTCTTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc        >  1:247796/1‑234 (MQ=255)
                                                                                                                                                                                                          aGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTTTACTTACGCATGACGAATCATTAAAATTATGCTCAGATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc        <  2:247796/234‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
AGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAAC  >  NC_000913/2187874‑2188315

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