Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,014,967 Δ1 bp coding (91/771 nt) yqeC ← uncharacterized protein YqeC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,014,9670T.96.7% 125.2 / ‑3.6 30coding (91/771 nt)yqeCuncharacterized protein YqeC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (13/16);  total (13/17)
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.91e-01

CGTTGCACCAGTGCATCAATCGCTTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTTTCCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGCCATCGACTGTATAATTACCGCTGCCGGATTATCATCAAGAATGGGGCAA  >  NC_000913/3014723‑3015216
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
cGTTGCACCAGTGCATCAATCGCTTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCa                                                                                                                                                                                                                                                    <  1:238644/251‑1 (MQ=255)
                 aatcccgttgggCGTAAATCCCTGCACGTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAAGGAAATGGGGGATGTAAGAGGCGCATAAGGAAGCATGGCGGTATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTTTCCCCC‑AGCACCAACAACATAAAt                                                                                                                                                                                                                                   <  2:216069/241‑1 (MQ=255)
                     gCTTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATAGGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATcaca                                                                                                                                                                                                                               <  2:234843/251‑1 (MQ=255)
                                        gCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAAGGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCAc                                                                                                                                                                                                            <  2:102102/251‑1 (MQ=255)
                                         cACTTTTCCCTGGTTAGATTTCCAGCTGTGAAAACAACATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGCTCACGAAAGAAAACCACGGGCCAGTGAGACATTGGCATAAACATATGTGTCGTCGTCGTGATTAACACCCGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCACAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCAcc                                                                                                                                                                                                           >  1:187413/1‑251 (MQ=255)
                                          aCTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGCGACGTTGGCCTAAACATATGTGTCATCGGCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCATCACCAACAACAGAAATCACAGTGGGACGTTTCTgaccaccc                                                                                                                                                                                                          >  2:180296/1‑244 (MQ=255)
                                                cccTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATc                                                                                                                                                                                                     >  2:186138/1‑250 (MQ=255)
                                                 ccTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGAAAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCa                                                                                                                                                                                                    >  1:362778/1‑250 (MQ=255)
                                                        cGCTTTCCAGATGGGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGGGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCCCCTAAATCAATGACTa                                                                                                                                                                                             <  2:250222/250‑1 (MQ=255)
                                                        cGCTTACCAGCTGTAAAAACAAAAGGAAAGGGGGGATGTAAGAGATGCATGAGGATTCATGGCGGGATCACGACAGAAACCCACGTGCCAGTAAGACTTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGTATCAGTGCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCAAAGTTCGACGTTTCTGCGCACCTAAATCATTGGCTaa                                                                                                                                                                                            <  2:223354/251‑1 (MQ=255)
                                                                           cAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACtat                                                                                                                                                                          >  2:109558/1‑250 (MQ=255)
                                                                                   aatgGGGGATGTAAGAGACGCATGAGAAAGCATGGCGGGATCACGTCAGAAACCCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCAc                                                                                                                                                                 <  2:288893/251‑1 (MQ=255)
                                                                                        gggATGTAAGAGACGCATGAGGAAGCATGGGTGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGGCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCTATGGGGCAACTAGGCTTGTCACGAGTc                                                                                                                                                            <  2:149534/251‑1 (MQ=255)
                                                                                           aTGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAACAACACGGGCCAGTGAGACGATGGCATAAACCTATGTGTCGTCGTCGTGATTAACACAAGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCCACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCttt                                                                                                                                                         >  1:375730/1‑251 (MQ=255)
                                                                                                  gggcgcgtGAGGAAGAAGGGCGGGATCACGAAAGAAAACCACGGGCCAGTGAGCCGTGGGCAGAAACAGTGGTGTCGTTGTAGTGATTAACACCCCCCTGCCGCTAGCCGGGAGCTGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGGGGGACGTTTCTGCGCACCTAAATCAAGGACTAATGCCGATGGGTCAACTATACTGGTCACTAGTCTGTATGACct                                                                                                                                                  <  2:207390/244‑1 (MQ=255)
                                                                                                    aCGCATGAGGAAGCATGGCGGGATCACGCCAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctct                                                                                                                                                >  1:309303/1‑251 (MQ=255)
                                                                                                     cGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACctctc                                                                                                                                               <  1:109558/251‑1 (MQ=255)
                                                                                                                     gagggATCGCGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATGTGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCGGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAAGCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTGTAGGACCTCTCTGGGTTAAATGATccc                                                                                                                               <  2:208676/249‑1 (MQ=255)
                                                                                                                     gCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACTTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATccc                                                                                                                               <  2:176993/251‑1 (MQ=255)
                                                                                                                                     aaaCCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATca                                                                                                                                                                                                                                 >  2:167065/1‑137 (MQ=255)
                                                                                                                                     aaaCCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATca                                                                                                                                                                                                                                 <  1:167065/137‑1 (MQ=255)
                                                                                                                                              gCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACg                                                                                                      <  1:44615/251‑1 (MQ=255)
                                                                                                                                                                           gacgtagtcggGATTAAAACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAAAAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACATCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCtatc                                                                          >  2:221800/7‑249 (MQ=255)
                                                                                                                                                                                                 gCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGc                                                   <  2:89624/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTa                                                                                                         <  2:24865/192‑1 (MQ=255)
                                                                                                                                                                                                      ccGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTa                                                                                                         >  1:24865/1‑192 (MQ=255)
                                                                                                                                                                                                        gCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTAt                                                                                               <  1:359796/200‑1 (MQ=255)
                                                                                                                                                                                                        gCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTAt                                                                                               >  2:359796/1‑200 (MQ=255)
                                                                                                                                                                                                                               gccaAAAAAGCAAGCTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGCCATCGACTGTATAATTACCGCTGCCGGAtt                     >  1:118864/1‑251 (MQ=255)
                                                                                                                                                                                                                                         caagcTGGTTT‑CCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGCCATCGACTGTATAATTACCGCTGCCGGATTAtcatca              >  2:65415/1‑248 (MQ=255)
                                                                                                                                                                                                                                                ggtttCCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGCCATCGACTGTATAATTACCGCTGCCGGATTATCATCAAGAATgggg     >  1:324/3‑251 (MQ=255)
                                                                                                                                                                                                                                                ggtttCCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGCCATCGACTGTATAATTACCGCTGCCGGATTATCATCAAGAATggg      <  2:118864/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    tCCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAAGTAGCAGGCGTTATTAGCCCTTAGAAATAATGGAAGCACGATGTCAGCCATCGACTGTATAATGACCGCTGCCGGATTattttgaagattgtggcaa  <  2:188543/250‑18 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
CGTTGCACCAGTGCATCAATCGCTTCTGGCGTAAATCCCTGCACTTTTCCCTGGTTCGCTTTCCAGCTGTGAAAACAAAATGAAATGGGGGATGTAAGAGACGCATGAGGAAGCATGGCGGGATCACGACAGAAAACCACGGGCCAGTGAGACGTTGGCATAAACATATGTGTCGTCGTCGTGATTAACACACGCCTGCCGCTAGCCTGGAGCAGTTCTGCCAGCCAAAAAAGCAAGCTGGTTTTCCCCCCAGCACCAACAACAGAAATCACAGTGGGACGTTTCTGCGCACCTAAATCAATGACTAATGCCGATGGGTCAACTATACTTTTCACGAGTCTTTATGACCTCTCTGGGATAAATTATCCCCAACTTAATCCATCAGGAAGTAACGCAATTATCAGGCGTTATTAGCCCCTATAAATAATGGAACCACTATGTCAGCCATCGACTGTATAATTACCGCTGCCGGATTATCATCAAGAATGGGGCAA  >  NC_000913/3014723‑3015216

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