Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 587,362 A→G V207A (GTC→GCC)  ybcH ← DUF4434 domain‑containing protein YbcH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913587,3620AG96.9% 99.3 / ‑4.3 32V207A (GTC→GCC) ybcHDUF4434 domain‑containing protein YbcH
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (15/16);  total (16/16)
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.92e-01

GTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAAC  >  NC_000913/587246‑587496
                                                                                                                    |                                                                                                                                      
gTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACg                                                                                                                 <  1:133223/140‑1 (MQ=255)
   cATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACgtgt                                                                                                              >  1:44690/1‑140 (MQ=255)
        acaacGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCcag                                                                                                         <  1:386089/140‑1 (MQ=255)
           acGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGTGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCcagcag                                                                                                      >  2:753730/1‑140 (MQ=255)
           acGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCcagcag                                                                                                      <  2:44690/140‑1 (MQ=255)
            cgccgcTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGt                                                                                                     <  1:635514/140‑1 (MQ=255)
             gccgcTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGtt                                                                                                    >  2:307745/1‑140 (MQ=255)
                 cTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCg                                                                                                >  2:646142/1‑140 (MQ=255)
                        gggCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCGGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCa                                                                                         >  2:632294/1‑140 (MQ=255)
                         ggCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCAt                                                                                        >  2:630000/1‑140 (MQ=255)
                         ggCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCAt                                                                                        <  2:378531/140‑1 (MQ=255)
                         ggCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCAt                                                                                        <  2:335283/140‑1 (MQ=255)
                         ggCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCAt                                                                                        <  2:325561/140‑1 (MQ=255)
                             ggACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCGGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCggg                                                                                    <  1:632294/140‑1 (MQ=255)
                                   tGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTATACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACAt                                                                              >  2:513168/1‑140 (MQ=255)
                                            gcgcTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGc                                                                     <  1:513168/140‑1 (MQ=255)
                                                         aaTAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTg                                                        <  2:505987/140‑1 (MQ=255)
                                                                 tCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAAc                                                >  2:620740/1‑140 (MQ=255)
                                                                     gCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGt                                            >  2:180736/1‑140 (MQ=255)
                                                                     gCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGt                                            <  1:180736/140‑1 (MQ=255)
                                                                         ttCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTg                                        >  2:27068/1‑140 (MQ=255)
                                                                                  cAGTTTATCCACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACAt                               <  1:630000/140‑1 (MQ=255)
                                                                                          ccACGCCGCTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATc                       >  2:565879/1‑140 (MQ=255)
                                                                                             cgccgcTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGc                    >  2:44052/1‑140 (MQ=255)
                                                                                              gccgcTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGaaa                                  >  2:319821/1‑125 (MQ=255)
                                                                                              gccgcTGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGaaa                                  <  1:319821/125‑1 (MQ=255)
                                                                                                   tGCCATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCTAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTgc              >  2:291624/1‑140 (MQ=255)
                                                                                                      cATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGcgcg           <  1:307745/140‑1 (MQ=255)
                                                                                                      cATCCTGTACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGcgcg           <  1:481779/140‑1 (MQ=255)
                                                                                                             tACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTa    >  2:193417/1‑140 (MQ=255)
                                                                                                             tACCCAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTa    <  1:193417/140‑1 (MQ=255)
                                                                                                               cccAGGCATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAAc  >  2:688838/1‑140 (MQ=255)
                                                                                                                    |                                                                                                                                      
GTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAAC  >  NC_000913/587246‑587496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑