Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,244,908 C→T L61L (CTG→TTG)  ymgE → PF04226 family protein YmgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,244,9080CT95.0% 49.8 / ‑4.1 21L61L (CTG→TTG) ymgEPF04226 family protein YmgE
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (5/14);  total (5/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 = 4.69e-01

TACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATCTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCC  >  NC_000913/1244668‑1245129
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                               
taCTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCtt                                                                                                                                                                                                                       >  1:122863/1‑251 (MQ=255)
       agCGTTGATGTAAGCATGTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCtggtg                                                                                                                                                                                                                <  1:267588/251‑1 (MQ=255)
          gTTGATGTAAGCACTGTGTGCGTGCTATGCATCGAAACGAGGGCTTGAACATGGGAATTATGGCCTTGAGTATTTGTGTCCTGAGAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTGTCCTGACCTGTATTCTCGGGATAGGCGGTGCGGTGGTCGGCTGCTGGATGGCGACCATGTTTGGCATTGGCGGGTCCATCATTGGTTTTAATTTGCACAGCTTCCTGGTGGCg                                                                                                                                                                                                             <  2:186005/251‑1 (MQ=255)
                  aagtattgctggcttCCGATTAATCGAAACGCGTGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGAGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGtggtgg                                                                                                                                                                                                      <  2:122863/236‑1 (MQ=255)
                            tgctgccGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCGGATAGCCGGCATTATCTCCAAGCTAATCATGCCGGGTCTTGATGGTGGGGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGGTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATtc                                                                                                                                                                                             <  2:320641/244‑1 (MQ=255)
                                            aaaCGAGGGCCCGAATATGGTAATGATTGCCTGGATTATTTTTGACCTGATAGCCGGCATGATAGTCAAGCTTATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGTATAGTAGGGGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTAGTTTTAATGTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTcc                                                                                                                                                                          <  2:78306/251‑1 (MQ=255)
                                                                  attattGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTACTGGTGGCGGGGGGGGGAACTATTCTTGTTCTGGGCATATTCCGCCTCCTGCGGAGAGGATa                                                                                                                                                       >  1:293425/1‑248 (MQ=255)
                                                                    tattGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGAt                                                                                                                                                   <  1:323673/250‑1 (MQ=255)
                                                                       tGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTCATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGCGGTTTTAATTTGCACAGCTTCCTGGTGGCGGGGGTGGGAGCCTTTTTCTGTTTGGGCCTATTCCGCCTCCTGGGGAGAGAATAAGCttttt                                                                                                                                               >  1:127858/1‑250 (MQ=255)
                                                                             gATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAgg                                                                                                                                         <  1:108126/251‑1 (MQ=255)
                                                                                    tttGACCTGATAGCCGGCATTATCGTCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATa                                                                                                                                   <  1:221893/250‑1 (MQ=255)
                                                                                     ttGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCGGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAg                                                                                                                                  <  1:366514/250‑1 (MQ=255)
                                                                                                gCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGc                                                                                                                      <  1:266912/251‑1 (MQ=255)
                                                                                                                                           gtggAGTTTTCCTGACCGGTATGCTGTGGATAGTCTGTGCGGTGGTCGGCGGCGGGCTGGCGACCATGTTTGGCATTGGCGGCGCCATCAGTGGTTTTAATCTGCACAGCTTCCTGGTGTCGGTGGTGGGAGCTATTCTCGTGCTGGGCATATTCCGCCTTCTGCGAAGAGAATAAGATTTTCAGAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCtt                                                                           <  1:201299/251‑1 (MQ=255)
                                                                                                                                                  tttCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGTTTTTAATTTGCACAGCGTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGAGGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACggg                                                                     <  1:43730/250‑1 (MQ=255)
                                                                                                                                                   ttCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGAGGCTTTAATTTGCACAGCTTCCTGGTGGCGGTGGGGGGAGCTATTCTCGTTCTGGGCCTATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAAGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCCGGACGTATACTCGACACCTGCTTTAAggggg                                                                   >  1:332082/1‑249 (MQ=255)
                                                                                                                                                                                                tggctggcGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGGGCTATTCTCGTTCTTGGCATATTCCGCCTCCCGCGGAGAGAATAAGATTTTCATAAGGCGGATAGCGATAAAGATGCCGCTATCCGCTTTCACATCTGAACGTTTACTCGACACCCGCTTTACAGGTGAAAAAAATCAATATAGTCACTCTATATCCGAAAAAATCGGtgtcc                     >  2:131047/1‑247 (MQ=255)
                                                                                                                                                                                                 ggttgtctaCCATGTTTGGCATGGTCGGCTACGTAAGTGTTTTTAATTTGCAAAGCTGCCTGGTGGCTGGGGTGGGATATATGTTCGTTCGGGGCATATTCCGCCTCCTGCGAATAGTATATGATTTTCATAAGGCTGTTAGGGATACAGATGCCGCTATACGCGTTCACATCATAACGTCTACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATCAAGCCGGTTcc                      <  2:146534/241‑1 (MQ=255)
                                                                                                                                                                                                 ggctggcGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTGCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTccc                     <  1:179880/250‑1 (MQ=255)
                                                                                                                                                                                                              gTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAacaacaaca        <  1:163933/250‑1 (MQ=255)
                                                                                                                                                                                                                   gCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATcc  <  1:78187/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                               
TACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATCTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCC  >  NC_000913/1244668‑1245129

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