Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,709,919	G→T	32.9%	Y336* (TAC→TAA)	oatA ←	possible acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,709,919	0	G	T	32.9%	44.5 / 14.4	32	Y336* (TAC→TAA)	oatA	possible acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base T (2/9);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 6.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACTGTTCTAGCAAATTGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGAA  >  CP000730/2709780‑2710057                                                                                                                                            |                                                                                                                                           acacTTTTCTAGCAAATTGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAAGAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACg                                                                                                                                            <  1:100502/140‑1 (MQ=255)     ctttCTAACAAATGGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTtata                                                                                                                                        <  1:161410/138‑1 (MQ=255)                 ttGCCCCTACTTTTTAAGTAAAAAAGCAAAACCTTAAAATCCTTTTTTACTTTTATGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTc                                                                                                                            <  1:9925/140‑1 (MQ=255)                      ccTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGt                                                                                                                       <  2:144304/140‑1 (MQ=255)                             tttAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGtaat                                                                                                                >  2:301540/1‑140 (MQ=255)                                     aaaaTGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTc                                                                                                        >  1:292953/1‑140 (MQ=255)                                               aGCTTTAAATCCTTTTTTACGTATAGGTTTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACaa                                                                                              <  1:307178/140‑1 (MQ=255)                                                  tttAAATCCTTTTTTACCAAACGGTTTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGa                                                                                           <  2:324589/140‑1 (MQ=255)                                                     aaaTCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATACCTTTTCACAAAAACAATGAtcg                                                                                        >  2:190600/1‑138 (MQ=255)                                                      aaTCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATACAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAgg                                                                                       <  1:53189/140‑1 (MQ=255)                                                          cTTTTTTATGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATaa                                                                                   <  1:267318/140‑1 (MQ=255)                                                              tttACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTCACAAAATGTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTATTAACTGTTAACAAAAACAATGATAGGATAATGcc                                                                               >  2:24611/1‑140 (MQ=255)                                                                tACGTATAGTTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACTGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCAt                                                                             >  2:62756/1‑140 (MQ=255)                                                                tACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCAt                                                                             <  1:203259/140‑1 (MQ=255)                                                                tACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGATAACAAAAACAATGATAGGATAATGCCAt                                                                             <  2:177479/140‑1 (MQ=255)                                                                     ataGGTGTTTGAATAAGGGGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATAGCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGta                                                                        <  2:25449/140‑1 (MQ=255)                                                                           gtTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCt                                                                  >  2:95933/1‑140 (MQ=255)                                                                               cAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCAATAAGTATAAGCTaaat                                                              >  2:123171/1‑140 (MQ=255)                                                                                 aTAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGCTAATGCCATAAGTATAAGCTATATGa                                                            >  1:21655/1‑140 (MQ=255)                                                                                 aTAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGa                                                            <  1:42153/140‑1 (MQ=255)                                                                                    aaGCGATACGAAATTTCAGCCATTAACGCTTTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCg                                                         <  2:150287/140‑1 (MQ=255)                                                                                    aaGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCg                                                         <  1:301540/140‑1 (MQ=255)                                                                                       cGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGttt                                                      <  2:292953/140‑1 (MQ=255)                                                                                          tACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATACCTGTTCACAAAAACAATGATAGGATAATGCCATAATTATCAGCTATATGATCGTTTAcc                                                   >  2:13011/1‑140 (MQ=255)                                                                                           accaaaTTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACTTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCt                                                  <  1:93348/137‑1 (MQ=255)                                                                                           aCGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCt                                                  >  1:113722/1‑140 (MQ=255)                                                                                                    catctatcatcgtggTTAACAAACTTTCTATAATATACACTTATACCGGCATTTGTCCTTGTACGTAATAACTTTTATCATATATAAGGGTAGTATAATGCCATAAGTGTAAGCTATATGATCGTTTATCTATAATTAGt                                         <  2:75623/126‑1 (MQ=255)                                                                                                     aGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTa                                        <  2:5150/140‑1 (MQ=255)                                                                                                             aCGCTGTTAACAAAATTTCTCTAACATAAACTTATACCGGTATTTGTCCTTGGACCTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTc                                <  2:57168/140‑1 (MQ=255)                                                                                                                              tCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAATTTTAGc               >  2:238268/1‑140 (MQ=255)                                                                                                                                      ataaACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCTTTTACCTATAATTAGTAAATGTTTCATACTTAAAAATTTAGCAAATAAAc       >  1:49914/1‑140 (MQ=255)                                                                                                                                           cGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGaa  >  2:107119/1‑140 (MQ=255)                                                                                                                                           cGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGaa  >  1:335593/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                           ACACTGTTCTAGCAAATTGCCCCTTCTTTTTAGGTAAAAATGCAAAAGCTTTAAATCCTTTTTTACGTATAGGTGTTTCAATAAAGCGATACGAAATTTCAGCCATTAACGCTGTTAACAAAATTTCTATAATATAAACGTATACCGGTATTTGTCCTTGTACGTAATAACTGTTCACAAAAACAATGATAGGATAATGCCATAAGTATAAGCTATATGATCGTTTACCTATAATTAGTAAAGGTTTCATACTTAAAAATTTAGCAAATAAACTAGAA  >  CP000730/2709780‑2710057

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