Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,325,037:1	+T	11.4%	coding (282/1323 nt)	atoE →	short chain fatty acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,325,037	1	.	T	11.4%	119.8 / 12.6	35	coding (282/1323 nt)	atoE	short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base . (11/20);  new base T (2/2);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

TGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGC  >  NC_000913/2324895‑2325177                                                                                                                                                |                                                                                                                                             tGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACa                                                                                                                                         <  3:376867/148‑1 (MQ=255)   aaaaTGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAggg                                                                                                                                      <  1:235035/149‑1 (MQ=255)       tgtgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGt                                                                                                                                   <  3:164140/148‑1 (MQ=255)              tGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTgg                                                                                                                           >  2:339138/1‑149 (MQ=255)                    tttCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACt                                                                                                                      >  3:21976/1‑148 (MQ=255)                    tttCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACt                                                                                                                      <  3:207428/148‑1 (MQ=255)                        tGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttctt                                                                                                                 >  2:109111/1‑149 (MQ=255)                        tGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttct                                                                                                                  <  4:366451/148‑1 (MQ=255)                         ggAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                <  2:86462/149‑1 (MQ=255)                           aaCTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCgg                                                                                                              <  1:339138/149‑1 (MQ=255)                                     cGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGc                                                                                                     <  4:195963/148‑1 (MQ=255)                                       tttGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTctt                                                                                                                                <  3:268576/119‑3 (MQ=255)                                       tttGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTctt                                                                                                                                >  4:268576/1‑117 (MQ=255)                                                             tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttgg                                                                             <  3:136078/148‑1 (MQ=255)                                                              catcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtt                                                                           <  1:109111/149‑1 (MQ=255)                                                              catcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggt                                                                            <  3:306656/148‑1 (MQ=255)                                                                 cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtttg                                                                         <  4:21976/148‑1 (MQ=255)                                                                  gTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATttggtttggt                                                                       >  1:318178/1‑149 (MQ=255)                                                                        ggTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCg                                                                  <  3:133907/148‑1 (MQ=255)                                                                                 cTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGt                                                         >  3:324229/1‑148 (MQ=255)                                                                                  ttGACAGCTCTGCTCCGGTGAGAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTGCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGtt                                                        <  3:104562/148‑1 (MQ=255)                                                                                   tGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGttt                                                       >  3:47198/1‑148 (MQ=255)                                                                                     ccAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGcc                                                    >  1:22920/1‑149 (MQ=255)                                                                                           ctGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGa                                               <  3:147488/148‑1 (MQ=255)                                                                                            tGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAg                                             >  2:60522/1‑149 (MQ=255)                                                                                             gCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGt                                                 <  4:381147/144‑1 (MQ=255)                                                                                             gCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGt                                                 >  3:381147/1‑144 (MQ=255)                                                                                                        aaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                 >  2:1976/1‑149 (MQ=255)                                                                                                               ctgcGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCGGGCGacttcacg                           <  4:47198/148‑8 (MQ=255)                                                                                                                         ccgccTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTAt                <  1:246828/149‑1 (MQ=255)                                                                                                                          cgccTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATc               >  1:169218/1‑149 (MQ=255)                                                                                                                              tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt           <  2:98261/149‑1 (MQ=255)                                                                                                                              tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt           <  2:365285/149‑1 (MQ=255)                                                                                                                                      catgACGCCCTGTAC‑GGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTg   >  1:361515/4‑149 (MQ=255)                                                                                                                                       aaGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGc  <  2:318178/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             TGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGC  >  NC_000913/2324895‑2325177

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