Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,000,980 A→T 13.9% L26H (CTT→CAT)  yigE ← DUF2233 domain‑containing protein YigE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,000,9800AT13.9% 69.3 / 2.1 29L26H (CTT→CAT) yigEDUF2233 domain‑containing protein YigE
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base T (2/2);  total (13/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.88e-01

TGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACAC  >  NC_000913/4000744‑4001227
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
tGTACAAACCGAGCGGCTCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCCCCTGACCCTGACTATTAATATCCGCCAGCAGTTCATGTAACCTTCCCCACGCTTCGACATTTGGTTGTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGAGTAACGGTATACGCCTGTACGGTCATCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAg                                                                                                                                                                                                                   <  1:393060/275‑1 (MQ=255)
                                                               gaccctgacCCCTAATATCCGCCAGCAGCGCATGTAACCTTCACCACGCTTCGCCATTGGCTCTGTGCCAGCACATTTTCACCCGCTCCGTTAGCAGCTAACCGGTATTCGCCTGTACGGTCAGCGTCGCACCTGAACGTGCGCAATCATCTGCGGCAACAGCCAACAAGGGACCTAAGGTGAGGGCGCGACAACTCCGTCTGAGATACAAGGTGATCATTCCTTTACCAATGAGTAACTGATGCGACATGAGAGGTCCTGGGGGTGCGAgtttt                                                                                                                                                    >  2:73916/1‑275 (MQ=255)
                                                                      acTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGAtt                                                                                                                                                                                                                       >  1:70128/1‑201 (MQ=255)
                                                                          ttAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGATATTCAAg                                                                                                                                                                                                                   <  2:70128/201‑1 (MQ=255)
                                                                           tAATATCCGCCAGCAGAGCAGGTAACGTTCCCCACGCTTCGACATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGGGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTc                                                                                                                                                                                                         <  1:20951/210‑1 (MQ=255)
                                                                                            gCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAg                                                                                                                                                                                                                                                        <  2:173527/146‑1 (MQ=255)
                                                                                                       tCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGAtt                                                                                                                                                                                                                       >  2:156916/1‑168 (MQ=255)
                                                                                                           cACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAg                                                                                                                                                                                                                   <  1:156917/168‑1 (MQ=255)
                                                                                                                 tCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGGTTGAGGATTAACGGTATACGCCTGTACGGGCAGCGTCGGATCTGGGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGGGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCACCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTa                                                                                                  <  1:73916/275‑1 (MQ=255)
                                                                                                                                                                   aaCGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGATGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGa                                                                                                                                                                                     >  2:146464/1‑142 (MQ=255)
                                                                                                                                                                   aaCGGTATACGCCTGTACGGTCAGCGGCGGGTGTGGGAGTGCGCAATCATCCGCATCAACAGCAAACAACGGATGTAAGGGGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGAACATGCCTTTACCAATGAGTAGCTGa                                                                                                                                                                                     >  2:146453/1‑142 (MQ=255)
                                                                                                                                                                       gTATACGTCTGTACGGGCCGCGTCGGATCTTAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGATGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCCAGGTGATCATTCCTTTACCAATGAGTAGCTGATgcg                                                                                                                                                                                 <  1:146454/142‑1 (MQ=255)
                                                                                                                                                                       gTATACGCCTGTACGGGCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGATGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATgcg                                                                                                                                                                                 <  1:146465/142‑1 (MQ=255)
                                                                                                                                                                               cTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGCCAACAGCAAACAAGGGAAGTAAGGTGGGGGCGAGAAACATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTAAGGCGCCATCATAGATCCTGGATGAGGGCTTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAATTTCa                                                                                                      >  1:143136/1‑209 (MQ=255)
                                                                                                                                                                                  tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACACGGGAAGTAAGGTGAGGGCGAGAACAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTGGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTTTTGCCGACAAACTTTTGCTGGCGTAACATGCGCGCACGAt                                 >  1:239788/1‑275 (MQ=255)
                                                                                                                                                                                  tACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACATTGACGAGTAACAAAAAGCGCGAAGTTCAACTAt                                                                                                 >  1:165126/1‑211 (MQ=255)
                                                                                                                                                                                  tACGGTAAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCTTGTACCAATGAGTAGCTGATGCGCCATTATAGGTCTTGGATGTGGGATTTTTTTATCCTGTTAGCGAACTTGACGAGTACCAAAAAACGCGAAGTTCa                                                                                                      >  1:345051/1‑206 (MQ=255)
                                                                                                                                                                                   aCGGTAGGCGTCGGATCCGGGGGTGCGCAATCATCAGCGGCAACAGCAAAGAAGAGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGGTGCGCCATTATAGGTCCTGGATGTGGGATTGTTTTATCCTTTTAGCGACCTTGACTAGTACCAAAAGGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGAt                                 <  1:148096/274‑1 (MQ=255)
                                                                                                                                                                                      gtaaGCGTCGGATCAGAGAGTGCGCAATAATCAGCGTGTACAGCAAACAAGTGAAGTAAGGTGAGGGCGAGNAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAACGAGTAGCTCATGAGCCATTATAGGTCCTGGATGTGGGATTTGTTTATCCTGTTATCGACCTAACCGGGTACAAAAAAGCNCGAATTTCAACTa                                                                                                  <  2:345050/203‑1 (MQ=255)
                                                                                                                                                                                      gTCAGCGTCGGATCTGAGAGTTCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGAGCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGAt                                 <  1:189329/271‑1 (MQ=255)
                                                                                                                                                                                      gTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTGACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGNCGGCAAAATTTTGCTGGCGTAACATGCGCGCACGAt                                 <  1:150164/271‑1 (MQ=255)
                                                                                                                                                                                           cGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGTAAAATCCGTTTGAGATGCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATGTTTGTAGCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAGTTTTg                                                        <  1:57148/243‑1 (MQ=255)
                                                                                                                                                                                                            cgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAgg                                                                                                                                                                                                                  >  1:170334/1‑72 (MQ=255)
                                                                                                                                                                                                            cgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTACTGGATGTGGGATTTTTTTATCCTGTTAGCGACCt                                                                                                                                  >  1:107476/1‑152 (MQ=255)
                                                                                                                                                                                                             gcAATCATCAGCGGCAACAGCAAAAAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGGGATTGAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATTTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGac      <  1:370063/275‑1 (MQ=255)
                                                                                                                                                                                                                atcatcAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGAc                                                                                                                              <  2:107476/152‑1 (MQ=255)
                                                                                                                                                                                                                  agatggCGGCAAAAGCAAACAAATGAAGTAAGGTGGAGGCGGGAAAAATCCGGTTGGGATTCAAGGTTAGCATTCCTTTACCAAGTAGTAGCAGATCGGGCATTATAGGTACTGGATGTGGGATTTTTTTATCCTATACGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTGGACGGAAAAATTTTGCGGGCGTAATATGCGCGCACGATCACTCTAAGAGGACATTCTCCTTGGacacac  <  2:64502/269‑1 (MQ=255)
                                                                                                                                                                                                                             aCAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCt                                                                                                                                                                                       >  1:104100/1‑82 (MQ=255)
                                                                                                                                                                                                                               aGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGTGCGAAGTTCACCTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGATATTCGCCTTGGac      >  1:64502/1‑257 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
TGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACAC  >  NC_000913/4000744‑4001227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: