Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,029,851	A→G	8.7%	T105A (ACC→GCC)	yedA →	amino acid exporter for phenylalanine, threonine

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,029,851	0	A	G	8.7%	122.8 / 4.1	47	T105A (ACC→GCC)	yedA	amino acid exporter for phenylalanine, threonine
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGCGGCAATCCGTGGGGCGCGATTC  >  NC_000913/2029709‑2029998                                                                                                                                               |                                                                                                                                                    gcgcGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg                                                                                                                                                >  7:118/1‑148 (MQ=255) gcgcGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg                                                                                                                                                >  6:303726/1‑148 (MQ=255) gcgcGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg                                                                                                                                                >  4:172951/1‑148 (MQ=255)   gcgGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc                                                                                                                                              <  5:342528/148‑1 (MQ=255)            aaCTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTAc                                                                                                                                     <  2:260915/148‑1 (MQ=255)            aaCTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTAc                                                                                                                                     >  6:302937/1‑148 (MQ=255)                 cccccGCTTCGTCCGCTGCTCAATGCCGCGGTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCtg                                                                                                                                 <  1:278143/147‑1 (MQ=255)                  ccccGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcccc                                                                                                                                            >  4:88886/1‑135 (MQ=255)                  ccccGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcccc                                                                                                                                            <  3:88886/135‑1 (MQ=255)                  ccccGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCtgtg                                                                                                                               <  3:222461/148‑1 (MQ=255)                         aCGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAgc                                                                                                                        >  6:117239/1‑148 (MQ=255)                                  gctCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGttt                                                                                                               >  6:177243/1‑148 (MQ=255)                                          ccGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTaa                                                                                                       <  1:243353/148‑1 (MQ=255)                                            gcgcTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTaaaa                                                                                                     <  5:302937/148‑1 (MQ=255)                                              gcggATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAACCGCCTGTTTGGCATTAAAAcg                                                                                                   <  8:289253/145‑1 (MQ=255)                                               cTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAAcgc                                                                                                  >  3:250634/1‑148 (MQ=255)                                                   ttGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCaa                                                                                              <  3:172951/148‑1 (MQ=255)                                                   ttGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCaa                                                                                              <  1:162552/148‑1 (MQ=255)                                                      gCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACt                                                                                           >  6:127407/1‑148 (MQ=255)                                                      gCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACt                                                                                           >  4:47534/1‑148 (MQ=255)                                                      gCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTAACCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACt                                                                                           >  6:71428/1‑148 (MQ=255)                                                       ccTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTg                                                                                          >  3:317511/1‑148 (MQ=255)                                                          gTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGaa                                                                                       >  4:90101/1‑148 (MQ=255)                                                           ttATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGaa                                                                                       <  4:94597/147‑1 (MQ=255)                                                           ttATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGaa                                                                                       <  5:303726/147‑1 (MQ=255)                                                           ttATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGa                                                                                        <  3:173669/146‑1 (MQ=255)                                                              ttGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAAtgg                                                                                    <  5:177243/147‑1 (MQ=255)                                                                  tggctgTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGttt                                                                                                               >  6:155323/1‑116 (MQ=255)                                                                  tggctgTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGttt                                                                                                               <  5:155323/116‑1 (MQ=255)                                                                   ggctgTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATgggtgggt                                                                              <  8:118/148‑1 (MQ=255)                                                                    gctgTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTa                                                                             >  1:217229/1‑148 (MQ=255)                                                                          gggaaTGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCa                                                                       <  6:217358/145‑1 (MQ=255)                                                                            tAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTGGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCAtt                                                                     <  3:335563/148‑1 (MQ=255)                                                                               tGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCATCCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGGATTGCCTTTggg                                                                  <  2:122369/148‑1 (MQ=255)                                                                                         aCGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGaa                                                        <  3:47534/148‑1 (MQ=255)                                                                                                  gAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGc                                               >  6:284670/1‑148 (MQ=255)                                                                                                      aTCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCaa                                           >  7:108573/1‑148 (MQ=255)                                                                                                         aaaaTGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTCCCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGGGGGTATTGCCATTGGGCTTGCCGGAATCATCATGATCAATAg                                        >  6:163545/1‑148 (MQ=255)                                                                                                             tGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAGCCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGt                                    >  2:4795/1‑148 (MQ=255)                                                                                                             tGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAGCCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGt                                    >  4:25663/1‑148 (MQ=255)                                                                                                                       ggCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTaa                          <  8:108573/148‑1 (MQ=255)                                                                                                                        gCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAg                         <  5:284670/148‑1 (MQ=255)                                                                                                                         cATCGCCGCCGTAGTGGTTGCAGCCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGc                        <  3:25663/148‑1 (MQ=255)                                                                                                                         cATCGCCGCCGTAGTGGTTGCAGCCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGc                        <  1:4795/148‑1 (MQ=255)                                                                                                                            cgccgcCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGCgg                      >  7:227032/1‑147 (MQ=255)                                                                                                                            cgccgcCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGCgg                      >  3:49469/1‑147 (MQ=255)                                                                                                                            cgccgcCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGCgg                      >  8:210392/1‑147 (MQ=255)                                                                                                                                               aCCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGCGGCAATCCGTGGGGCGCGATTc  >  7:360233/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                    GCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCCCCTCTTTACCCTGTGCTTCAGCCGCCTGTTTGGCATTAAAACGCGCAAACTGGAATGGGTGGGTATTGCCATTGGGCTTGCCGGAATCATCATGCTCAATAGCGGTGGAAATTTAAGCGGCAATCCGTGGGGCGCGATTC  >  NC_000913/2029709‑2029998

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