Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	593,444	G→C	intergenic (+315/+102)	unknown → / ← unknown	hypothetical protein/UPF0758 protein YicR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	593,444	0	G	C	85.0%	40.8 / 3.5	20	intergenic (+315/+102)	unknown/unknown	hypothetical protein/UPF0758 protein YicR
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base C (6/11);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

CTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTT  >  1/593213‑593563                                                                                                                                                                                                                                         |                                                                                                                        cTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAg                                                                                                      >  1:793873/1‑251 (MQ=255)         cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:126024/251‑5 (MQ=255)         cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:622206/251‑5 (MQ=255)              tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         <  1:821071/251‑10 (MQ=255)              tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         <  1:783339/251‑10 (MQ=255)                    aGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTTTTTATGTTTATATGATGATTCTGCGTTAAATTTTAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAgtgt                                                                                                   >  1:322279/1‑235 (MQ=255)                      caaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattc                                                                                 <  1:543917/251‑18 (MQ=255)                       aaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCATCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTTAGGGGGGTTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              >  1:158898/1‑233 (MQ=255)                          caaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattcatac                                                                             <  1:859192/251‑22 (MQ=255)                                      gctgAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaatta                                                                                    <  1:482745/232‑15 (MQ=255)                                         gAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTTTGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              >  1:826309/1‑215 (MQ=255)                                                            tGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgtt                                                                                               >  1:775368/1‑196 (MQ=255)                                                                                     gCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         >  1:440256/1‑171 (MQ=255)                                                                                                        aaaTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:258670/156‑5 (MQ=255)                                                                                                          aTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat                                                                                      <  1:284319/162‑13 (MQ=255)                                                                                                                     gctCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat                                                                                      <  1:366774/149‑13 (MQ=255)                                                                                                                                                      aaaCGGCGGGTTAATACTCCGCCGTTTTTCTTTTCATGTACATGATGATTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaatgt  <  1:82718/203‑12 (MQ=255)                                                                                                                                                                           tcgtTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:72476/90‑5 (MQ=255)                                                                                                                                                                                                tgatgaTTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAacac              >  1:770647/1‑149 (MQ=255)                                                                                                                                                                                                          gctttAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaa     >  1:450714/4‑140 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                        CTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTT  >  1/593213‑593563

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