Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,927,532	A→C	39.9%	intergenic (+111/‑104)	PP_RS08755 → / → ygfA	tRNA‑Asp/5‑formyltetrahydrofolate cyclo‑ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,927,532	0	A	C	39.9%	43.3 / 6.9	45	intergenic (+111/‑104)	PP_RS08755/ygfA	tRNA‑Asp/5‑formyltetrahydrofolate cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base A (22/3);  new base C (2/13);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 3.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCGGTGTGCATGTCCGCCCGACGGAAAGCCTTAACGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGACCGAAAACGTGCCGCTCACCC  >  ACB122/5927403‑5927660                                                                                                                                  |                                                                                                                                  ggCCGGTGTGCATGTCCGCCCGACGGAAAGCCTTAACGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATTGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAC‑GGCggg                                                                                                                            >  4:311803/1‑136 (MQ=255)           cATGTCCGCCCGACGGAAAGCCTTAACGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGttt                                                                                                                  >  3:26715/1‑136 (MQ=255)                    ccGACGGAAAGCCTTAACGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCgcc                                                                                                          >  2:133137/1‑133 (MQ=255)                      gACGGAAAGCCTTAACGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACtt                                                                                                        >  3:151578/1‑135 (MQ=255)                          gAAACCCTAAACGCCCCCGGAAATCACCACTTGGAACTTTCGGGTTCATGGGAAACACCTATAGGTTTCTTTTCGGGTCGATTGAATTTTGTTCCCCCCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTc                                                                                                   <  8:11869/136‑1 (MQ=255)                           aaccccTTACCGCCCCCTGAAATCTCCACTTTTAACTTTGGGGTTAATGGGCTACACCTATAGGGTTGTTAGGGGGGAGCCTTATTTTTCTTGCCCGCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTc                                                                                                   <  6:117416/131‑1 (MQ=255)                            cggCCTAACCGCCCCCTCCAATCCCCACATTTAACTTTGGGGTTAATGGGCAACACCGATAGCGGTCTTATCGGGGAGCTTTATTTTTCTTGCCCCCCCCCAC‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCaa                                                                                                 <  4:104314/134‑1 (MQ=255)                                   aaCGCCCCCTGCAATCCCCACCTTGAACTTTCGGGTTCACGGGCTACACCGATCGCGGTCTTAGCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGt                                                                                          <  7:208629/136‑1 (MQ=255)                                     cGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTac                                                                                        >  7:123587/1‑136 (MQ=255)                                      ccccccTCCAAGCACCACCTTCAAGTTTGGGTTTAATGGGCAACCCGAATAGGGGTTTTTGGGGGGCGCCTGATTTTTCTGCCCCCCCCCCAG‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTacg                                                                                       <  7:396742/135‑1 (MQ=255)                                            tGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTAcgacgact                                                                                 >  3:104314/1‑134 (MQ=255)                                            tGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGGTTTGGCCGACTTTATTCAAGCTTTGTACGACGCTc                                                                                 >  3:366379/1‑136 (MQ=255)                                            tGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATAGGGGAGCCTGAATTTTCTTGCCATCACCAAA‑GGAGGGCAATTCGTTTTGGCCGACTTTAATCAAGCTTTGTACGAAGAta                                                                                 >  7:286945/1‑135 (MQ=255)                                              cAATCTCCCCCTTGAACTTGCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTTAATTTTCTTGCCCCCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATc                                                                                 <  7:300955/134‑1 (MQ=255)                                                      aCCTTGAACTTTCGGGTTCACGGGCTACCCCGATAGCGGTCTTAGGGGGGAGCCTGAATTTTCTTGCCCCCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAg                                                                       <  4:151578/136‑1 (MQ=255)                                                      aCCTTGAACTTGCGGGTTCATGGGCTACCCGGATAGCGGTCTTATCGGGGCGCCTTAATTTTCTTCCCCGCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCa                                                                        <  2:71019/135‑1 (MQ=255)                                                      aCCTTGAAATTTCGGTTTCACGGGCTACACCTATAGCGTTCTTATCGGGGAGCCTTAATTTTCTTGCCCGCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACcg                                                                            <  2:329078/131‑1 (MQ=255)                                                                ttCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGttag                                                             >  7:247613/1‑134 (MQ=255)                                                                 ccGGGTCCACGGGATACACCCATCGGGGTCTTAGCGGGGAGCCTCAATTTTCTTCCCCGCCCCCAC‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGt                                                            <  3:116206/135‑1 (MQ=255)                                                                          tgggCTACACCGATCGGGGTATTATGGGGTCGCCTTATTTTTCTTGCCGGCCCCCAA‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGc                                                    <  2:148055/134‑1 (MQ=255)                                                                          ggggCAACCCCTATCGGGGTCTTCGCGGGGACCCTTATTTTTC‑TGCCCCCCCCCAAGGGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGc                                                    <  6:359261/134‑1 (MQ=255)                                                                          aGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAT‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGcc                                                   <  8:123587/136‑1 (MQ=255)                                                                           gggCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCATATCGTATTGGCCGACTTTATTCAAGCATTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCa                                                  >  4:380819/1‑136 (MQ=255)                                                                                   ccGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAC‑GGCGGGCAATTCGTTTAGGCCGACTTTATTCAAGCTTTTTACGACGATCCCCCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCt                                          >  8:204894/1‑136 (MQ=255)                                                                                     gATAGCGGTTTTTGGGGGGAGCATTAATTTTCTTGCCCGCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCAACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTgg                                        <  3:171626/136‑1 (MQ=255)                                                                                     cttAGCGGTCTTTTGGGGTAGCGTTATTTTTCTTGCCCGCCCCCAG‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTgg                                        <  4:26715/134‑1 (MQ=255)                                                                                      aTAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTgg                                        >  8:243229/1‑135 (MQ=255)                                                                                       tAGCGGTCTTATCGGGGAGCCTTATTTTTCTTCCCCGCCCCCAC‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTgg                                        <  8:243016/134‑1 (MQ=255)                                                                                         gCGGTCTTATCGGGGAGCCTTAATTTTCTGGCCCGCCCCCAC‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATc                                    <  8:108953/136‑1 (MQ=255)                                                                                         gCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCCA‑GGCGGGCAATTAGTTTTGGCCGCCTTTATTCAAGCATTGTACGACGAGCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACAACTGGGATc                                    >  8:144088/1‑136 (MQ=255)                                                                                         gCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAATCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATc                                    >  7:1584/1‑136 (MQ=255)                                                                                         gCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGATTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATc                                    >  8:17286/1‑136 (MQ=255)                                                                                              cTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGccc                                >  6:177410/1‑135 (MQ=255)                                                                                               ttATCGGGGAGCCTGAATTTTCTTGCCTGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAg                                                                       <  1:312377/95‑1 (MQ=255)                                                                                               ttATCGGGGAGCCTGAATTTTCTTGCCTGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAg                                                                       >  2:312377/1‑95 (MQ=255)                                                                                               ttATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GCCGGGCAATTCGTTTTGGCCGACTTTATTCAACCTTTGTAAGACGATCCACCGCCCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGt                              >  4:51759/1‑136 (MQ=255)                                                                                                        aGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGTCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGAcc                     >  4:354774/1‑136 (MQ=255)                                                                                                        aGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAACCCACCACATCTGGGATCGCCCGTGCCATGAcc                     >  4:360806/1‑136 (MQ=255)                                                                                                          ccTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCCCCGCGCAGGCACTGGTCGTACTAAAGCCAGCACATCTGGGCTCGCCCGTGCCATGACCg                    >  3:61539/1‑135 (MQ=255)                                                                                                          ccTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTGTATTCAAGCTTTGTCCGCCGATCCCCCGCGCAGGCACTGGTCGTACTAAAGCCAGCACAACGCGGCTCGCCCGTGGCCTGACCg                    >  3:53277/1‑135 (MQ=255)                                                                                                           cTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCCCCGCGCAGTCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGGCATGACCg                    >  6:85693/1‑134 (MQ=255)                                                                                                            tGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTTGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGACCGaa                  >  8:168146/1‑135 (MQ=255)                                                                                                                        cccgcccCAAG‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCTCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGACCGAAAACGTGCCGCTc      <  5:273130/135‑1 (MQ=255)                                                                                                                          cccccccAC‑GGGGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGACCGAAAACGTGCCGCTCAcc   <  4:366379/134‑1 (MQ=255)                                                                                                                            cccccAG‑CGGGGCCAATTCGTTTTTGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGACCGAAAACGTGCCGCTCAccc  <  6:60921/135‑1 (MQ=255)                                                                                                                                  |                                                                                                                                  GGCCGGTGTGCATGTCCGCCCGACGGAAAGCCTTAACGCCCCCTGCAATCTCCACCTTGAACTTTCGGGTTCAAGGGCTACACCGATAGCGGTCTTATCGGGGAGCCTGAATTTTCTTGCCCGCCCCCAA‑GGCGGGCAATTCGTTTTGGCCGACTTTATTCAAGCTTTGTACGACGATCCACCGCGCAGGCACTGTTCGTACTAAAGCCAGCACATCTGGGATCGCCCGTGCCATGACCGAAAACGTGCCGCTCACCC  >  ACB122/5927403‑5927660

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: ‑