Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,316,131	T→C	11.9%	E358G (GAG→GGG)	trpE ←	component I of anthranilate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,316,131	0	T	C	11.9%	91.5 / 4.0	34	E358G (GAG→GGG)	trpE	component I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base C (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACGGGCA  >  CP009273/1315871‑1316384                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                              gacgacCTCCCGCCTCGACAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATGCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTa                                                                                                                                                                                                                                                 <  1:100375/275‑1 (MQ=255)                                    ggTACTTGCGGCTCACCGCTTTACGTCCTCATATTCATACAGGCGCGATAAGGGTGCAGGGCGTCAAGATCGTGACGGAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGGGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATctct                                                                                                                                                                                                              <  1:251375/274‑1 (MQ=255)                                                                          aCAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCgtg                                                                                                                                                                        >  1:109706/1‑275 (MQ=255)                                                                                 cgATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCCCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCgg                                                                                                                                                                  >  2:311964/1‑274 (MQ=255)                                                                                                                                              gCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCCCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGAc                                                                                                    <  1:311969/275‑1 (MQ=255)                                                                                                                                                              aCGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAAt                                                                                                                                                                                                                                >  1:581172/1‑135 (MQ=255)                                                                                                                                                               cGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTcc                                                                                                                                                                                                          >  2:317094/1‑156 (MQ=255)                                                                                                                                                                  gcAACGTGGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACgg                                                                                                                                                                                                                            <  2:581159/134‑1 (MQ=255)                                                                                                                                                                   cAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTc                                                                                                                                                                                                                                   >  2:453876/1‑127 (MQ=255)                                                                                                                                                                   cAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTg                                                                                                                                                                                                      <  1:317099/156‑1 (MQ=255)                                                                                                                                                                   cAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtggtg                                                                                                                                  >  1:195442/1‑224 (MQ=255)                                                                                                                                                                    aaCTGTGGTGAGATCGGCGACGTAGCAGCTGCCGGGGGTGCAAATACGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGGTCAGACAGCTATTTATGATCGGTACGCATTGCCAGTTCAATACGACGATCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCCGCAAGCTGTTAGCACTCAAAGTGGCGGCTGGTGGACACAGACTTGATCGAGTTTACCGGCGGCGAGCCCAATCGGGTGaccagat                                                                              >  2:230199/1‑269 (MQ=255)                                                                                                                                                                       tttGGTGAGATCGGCGACTTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat                                                                                                                              <  2:195440/224‑1 (MQ=255)                                                                                                                                                                       tttGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATa                                                                                                                                                                                                                               <  1:453883/127‑1 (MQ=255)                                                                                                                                                                                 tgggggaaggaGCGGCTGGCGGGGGTGTAAATGCGTGCCAGATCATTGCGGGCGAGATAAACAAGCATCAGATGTGCAGACAGCTATTTATGATCGGTACGCATTGCAAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATaa                                                                  <  1:185796/264‑1 (MQ=255)                                                                                                                                                                                   ggCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGa                                                                                                     >  2:723344/1‑237 (MQ=255)                                                                                                                                                                                       aCGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTTCGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGAcgcg                                                                                                 <  1:723358/237‑1 (MQ=255)                                                                                                                                                                                                  gCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGtt                                                                                                                                                                      >  2:195727/1‑157 (MQ=255)                                                                                                                                                                                                      ggggTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCgg                                                                                                                                                                  <  1:195729/157‑1 (MQ=255)                                                                                                                                                                                                            cAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCCCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAg                                                                                         >  1:203036/1‑224 (MQ=255)                                                                                                                                                                                                                tGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCCCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTg                                                                                     <  2:203034/224‑1 (MQ=255)                                                                                                                                                                                                                  cGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGcc                                                                                               >  1:727094/1‑212 (MQ=255)                                                                                                                                                                                                                    tGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATc                                                                                                                                                                                                                                                     >  1:140375/1‑60 (MQ=255)                                                                                                                                                                                                                    tGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATaaa                                                                 >  2:332400/1‑240 (MQ=255)                                                                                                                                                                                                                      ccAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtg                                                                                                                                     >  2:302149/1‑170 (MQ=255)                                                                                                                                                                                                                        aGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTa                                                                                                                                                                                                                                                 <  2:140374/60‑1 (MQ=255)                                                                                                                                                                                                                        aGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAAc                                                             <  1:332406/240‑1 (MQ=255)                                                                                                                                                                                                                          aTCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtggtgg                                                                                                                                 <  1:302154/170‑1 (MQ=255)                                                                                                                                                                                                                               tACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGAtctc                                                                                                                                                                                                               >  2:4235/1‑87 (MQ=255)                                                                                                                                                                                                                               tACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATc                                                                                                                                                             >  2:120822/1‑137 (MQ=255)                                                                                                                                                                                                                                   ggCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTc                                                                                                                                                                                                           <  1:4235/87‑1 (MQ=255)                                                                                                                                                                                                                                   ggCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGt                                                                                                                                                         <  1:120823/137‑1 (MQ=255)                                                                                                                                                                                                                                                tggttaaggnggtCAGGCAGTGCTGTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACgggca  <  2:272739/263‑1 (MQ=255)                                                                                                                                                                                                                                                           ttCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACg      >  1:272743/1‑260 (MQ=255)                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                              GACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAGGCCGCCAGCGGTGACGGGCA  >  CP009273/1315871‑1316384

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