Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 466,723 G→A A197V (GCA→GTA)  ybaE ← uncharacterized protein YbaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913466,7230GA95.2% 54.7 / ‑6.1 21A197V (GCA→GTA) ybaEuncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/9);  minor base T (0/1);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CTGTGAATTGTGTTAAGCGAAAAGGACCCGTGCCGATCAGTGGGAATTGCGGATGCGCCAGATGGCTGCAATAGCTCGCCAGCCGGTGCGCAAGCCAGTAATCAGGGCGATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCC  >  NC_000913/466487‑466966
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
ctgTGAATTGTGTTAAGCGAAAAGGACCCGTGCCGATCAGTGGGAATTGCGGATGCGCCAGATGGCTGCAATAGCTCGCCAGCCGGTGCGCAAGCCAGTAATCAGGGCGATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGc                                                                                                                                                                                                                                       >  1:35394/1‑251 (MQ=255)
                                         gggAATTGCGGATGCGCCAGATGGCTGCAATAGCTCGCCAGCCGGTGCGCAAGCCAGTAAGCAGGGCGATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCGTTTGCTTCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAAcc                                                                                                                                                                                               <  1:205951/250‑1 (MQ=255)
                                                       ngcCAGATGGCTGCAATAGCTCGCCAGCCGGTGCGCAAGCCAGTAATCAGGGCGATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAAcc                                                                                                                                                                                               >  1:51545/2‑236 (MQ=255)
                                                       cgcCAGATGGCTGCAATAGCTCGCCAGCCGGTGCGCAAGCCAGTAATCAGGGCGATGTAAAAAGAAGGCCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCTATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCg                                                                                                                                                                                                    >  1:504807/1‑231 (MQ=255)
                                                                                                          gCGATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTaaa                                                                                                                                                               >  1:5376/1‑217 (MQ=255)
                                                                                                                   aaaGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGc                                                                                                                    <  1:286010/251‑1 (MQ=255)
                                                                                                                              gACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTAAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCTCTGATTATTAGTATCGAATCTGGTCAGGCCGGAAAATATc                                                                                                         <  1:158298/251‑1 (MQ=255)
                                                                                                                                     nggATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGt                                                                                                                                                                                  >  1:313857/2‑171 (MQ=255)
                                                                                                                                                cTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGAAGTTGTAACAGCATCAATAATGGATGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGcc                                                                                                                                                                                                                                      >  1:222124/1‑108 (MQ=255)
                                                                                                                                                        cGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGt                                                                                                                                                                          >  1:442934/1‑160 (MQ=255)
                                                                                                                                                        cGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTGGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACg                                                                               <  1:245575/251‑1 (MQ=255)
                                                                                                                                                          nnnnnngCNNNNNAATNATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGt                                                                                                                                                                                  <  2:313857/144‑1 (MQ=255)
                                                                                                                                                                                      gCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCAAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTc                                                 >  1:4244/1‑251 (MQ=255)
                                                                                                                                                                                       cTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGNGGCNTTTANTANATCGCCGTTATGCCANTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTNNANNNTNCCCAGTGa                                                                                                                                                                        <  2:352347/131‑1 (MQ=255)
                                                                                                                                                                                       cTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAAAGAAGCTAAAAGACCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGa                                                                                                                                                                        >  1:352347/1‑131 (MQ=255)
                                                                                                                                                                                             gTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTc                                                                                                                                                                                           <  1:417937/106‑1 (MQ=255)
                                                                                                                                                                                                    cAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCg                                                                                                                                                                                                                                                <  1:6191/46‑1 (MQ=255)
                                                                                                                                                                                                    cAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTACATCg                                                                                                                                                                                                                                                >  2:6191/1‑46 (MQ=255)
                                                                                                                                                                                                                      ntgtAAGTGTGAGGCTTTTACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCgag                                                                                              >  1:404548/2‑174 (MQ=255)
                                                                                                                                                                                                                                gNNGCTTTTNNTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTaaa                                                                                                                                                               <  2:5376/99‑1 (MQ=255)
                                                                                                                                                                                                                                      tttACTACATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTgac  <  1:465897/250‑3 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
CTGTGAATTGTGTTAAGCGAAAAGGACCCGTGCCGATCAGTGGGAATTGCGGATGCGCCAGATGGCTGCAATAGCTCGCCAGCCGGTGCGCAAGCCAGTAATCAGGGCGATGTAAAAAGAAGGTCAGACACTGCGGATGGGTGACTTCAATACGCTTCACGCTAATAAATAATTGATCCAGTGCTGGCAGTTGTAACAGCATCAATAATCGCTGGTGTAAGTGTGAGGCTTTTACTGCATCGCCGTTATGCCAGTGTAGGGTTGAACGAAGATAAAAGTCCCAGCGTAACCCGTCAGTAGAGGTTTCCCAGTGATGCGCTAAATCGCCAATCGGGCGCTGAGTATTATTATCGAAGCGGGTCAGGCCGGAAAATATCTGCCCGGCGAGATGCTGCTCGGCACGGCCGGGCAAAAAGCCTGGTTGTAGCGGTTCGAGCGGGCGATAGTAGGGAATACGCAACGTGGGTGTATCGTTTTGCC  >  NC_000913/466487‑466966

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑