Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	426,038	A→C	44.2%	intergenic (+13/+34)	rssB_2 → / ← glcB	Regulator of RpoS/Malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	426,038	0	A	C	44.2%	0.8 / 3.5	34	intergenic (+13/+34)	rssB_2/glcB	Regulator of RpoS/Malate synthase G
Reads supporting (aligned to +/- strand):  ref base A (11/5);  new base C (10/2);  total (27/7)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GACACCGATGTGGAAGAGGCAGTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTACCGGCTCGGTATAACCGTTCGGTTGCTTGCCACCTTCCACT  >  ACB122/425906‑426159                                                                                                                                      |                                                                                                                          gACACCGATGTGGAAGAGGCAGTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGccc                                                                                                                         >  3:105237/1‑135 (MQ=255)      cGATGTGGAAGAGGCAGTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCGtta                                                                                                                     >  1:48337/1‑132 (MQ=255)      cGATGTGGAAGAGGCAGTGGATGTGATGCACCTGGGGGAGCAGGACTTCCGGCTCAAGCCGGTGGATCTGGGCAAGCTGATGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGGccccga                                                                                                                     >  7:173010/1‑130 (MQ=255)        aTGTGGAAGAGGCAGTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCGTAAAAc                                                                                                                  >  3:133232/1‑135 (MQ=255)                gagGCAGTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCGTAAAACGGTccccc                                                                                                          >  3:292367/1‑135 (MQ=255)                      gTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGGCCGGTAAACGGTCGCCCTgggac                                                                                                    >  8:179754/1‑132 (MQ=255)                        ggATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGGCCCGAAAAAGGGGACCCAGgcagcc                                                                                                  >  5:264108/1‑134 (MQ=255)                                     cTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCGTACa                                                                                     >  3:291960/1‑135 (MQ=255)                                      tGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAATATTGGATGATGCGCCAGGGGGCCCGGAAACAGGGGCCCCTGGCAGCCGCtgcggcgttgcc                                                                                    >  2:172129/1‑129 (MQ=255)                                       gggggTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGATTCTGCTTACAAc                                                                                   >  4:234775/1‑135 (MQ=255)                                       gggggTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACCGTGCCCCAGGCAGCAGATGCTGCTTACAAc                                                                                   >  3:302732/1‑135 (MQ=255)                                             ggtggACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGGAAAACCGTCCCCCTGGCCGCACATGCTGCTTACAACCCGttc                                                                             >  4:417297/1‑134 (MQ=255)                                                 gACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCGGAAAACGGTCCCCCTGGCAGCAGCTGCAGCGTACAACCCCTTAcgcg                                                                         >  2:257303/1‑135 (MQ=255)                                                            cAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGGAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCCTTACGCGCCTTGAACAcg                                                              >  8:225454/1‑135 (MQ=255)                                                                 cGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCCTAAAACGGTCCCCCAGGCCGCAGATGCTGCGTACAACCCCTTTCGCGCCTTGAACTCgcggcg                                                         >  7:63633/1‑135 (MQ=255)                                                                 cGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCgcggcg                                                         >  4:138665/1‑135 (MQ=255)                                                                  gggggATCGGGCCAAGCTGTGGGCGCTGGTCAAC‑AAAGAACTTAAATTTGGATGAGGCGCCCGGGGGACCTAAAACCGGCCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGa                                                        <  4:302732/132‑1 (MQ=255)                                                                       aTCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGGCCGTAAAAAGGGGCCCCTGGCAGCAGCTGGTGGGTACAACCCGGTACGCGCCGTGAACTCGCGGCGCcggcg                                                   >  8:42025/1‑135 (MQ=255)                                                                       aTCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGCTGCTGCTTACAACCCGTTACGCGCCGTGCACTCccggcgccggcg                                                   >  5:193881/1‑135 (MQ=255)                                                                       aTCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGCACTCGCGGCGCcggcg                                                   >  4:86342/1‑135 (MQ=255)                                                                          tGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCAGGCAGCAGAAGCTGCTTACAACCCCGTACGCGCCGTGAACTCGCGGCGAccgcgggg                                                >  3:396131/1‑133 (MQ=255)                                                                            ggCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACCGTCCCCCTGGCAGCACATGCGGCTTACCACCCCTTACCCGCCATGAACTCGCGGCCACGGCGGTGCa                                              >  7:110101/1‑135 (MQ=255)                                                                                tggtgtgGGCGCTGGTAAAC‑AAATAAATAAAAATTGGATCCTCCGCCTGGGGCACCTAAAAAGGGCCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTAc                                          <  6:199410/132‑1 (MQ=255)                                                                                aggtgtgGGCGCTGG‑GAACAAAATAACTAAAATTGGGATGAGGCGCCCGGGGCACCGAAAACCGGGCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTAc                                          <  4:292367/132‑1 (MQ=255)                                                                                  ttgttgtCGCTGGTCAAC‑AACCAAGTCAAATTTGGATGATGCGCCTGGGGGCCCGTAACCCGGGCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTACCg                                        <  7:369133/134‑1 (MQ=255)                                                                                  gtGTTGGAGCTGTTCAAC‑AAACAACTCAAAATTGGATCAGGCCCCCGGGGCACCGAAAAACGGGCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTACCg                                        <  2:40890/134‑1 (MQ=255)                                                                                    gTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGGCCGTTAAAAGGTCCCCCTGGCAGCAGATTCTGCTTAAAACCCGTTACGCGCCGTGAAATCGCGGCGCCCGCGGTGCAGttaccgc                                      >  5:199410/1‑129 (MQ=255)                                                                                              gTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGCCCGTAAAAACGTGCCCCAGTCAGCAGATGCTGCTTTCAACCAGTTTCGAGCCCTGAACTCGCGGCGACGGCGGTGCCGTAACGGCTCGGTATaaa                            >  6:350927/1‑134 (MQ=255)                                                                                                gAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGGCCGGAAAAAGGGCCCCCGAGCCCCCCATGCGGCTTTCAACCCCGTACCCCCCGTGAAAACGCGGGGCCCGCGGGGCAGTAACCGCCCGGGATTAc                            >  8:398310/1‑133 (MQ=255)                                                                                                gAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCGTGAACTCGCGGCGACGGCGGTGCAGTACCGGCTCGGTATAACCg                          >  3:212908/1‑135 (MQ=255)                                                                                                gAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGGCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGCCGGCGGGGCCGTACCGGCACGGTATAACCg                          >  8:346864/1‑135 (MQ=255)                                                                                                  cc‑AAACAACTCAAAATGGCATCAGGCCCCCGGGGGACCTTAAACCGGCCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTACCGGCTCGGTATAACCGtt                        <  3:234775/134‑1 (MQ=255)                                                                                                          aCTGAAAATTGGATGATGCGCCAGGGGGCCCCTAAAAACGGGCCCCAGGCAGCAGATGCTGCGTACAAACCCGTAAGCGCCTTGAAATAGCCGCGCCGGCGGGGCATGACCGGCGCGGTATAACCGTTCGGttgc                 >  2:138709/1‑135 (MQ=255)                                                                                                                 aTTGGATGATGCGCCAGGGGGCCCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCGTGAACTCGCGGCGCCGGCGGGGCAGTACCGGCGCGGTATAACCGTTCGGTTGCTTGGCAc          >  3:112682/1‑135 (MQ=255)                                                                                                                         atgCCCCTGGGGCCCCGTAACCGGGCCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTACCGGCTCGGTATAACCGTTCGGTTGCTTGCCACCTTCCACt  <  6:114087/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                          GACACCGATGTGGAAGAGGCAGTGGATGTGATGCACCTGGGGGTGGTGGACTTCCTGCTCAAGCCGGTGGATCTGGGCAAGCTGTTGGAGCTGGTGAAC‑AAAGAACTGAAAATTGGATGATGCGCCAGGGGGACCGTAAAACGGTCCCCCTGGCAGCAGATGCTGCTTACAACCCGTTACGCGCCTTGAACTCGCGGCGACGGCGGTGCAGTACCGGCTCGGTATAACCGTTCGGTTGCTTGCCACCTTCCACT  >  ACB122/425906‑426159

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