Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,568,885	T→A	9.1%	L245* (TTA→TAA)	fmtA →	FmtA protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,568,885	0	T	A	9.1%	138.3 / 5.9	44	L245* (TTA→TAA)	fmtA	FmtA protein
Reads supporting (aligned to +/- strand):  ref base T (15/25);  new base A (2/2);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AGGTTTAATTATCCAAAACGTTACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCTGATAAATACAAATCAT  >  CP000730/2568745‑2569022                                                                                                                                             |                                                                                                                                          aGGTTTAATTATCCAAAACGTTACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAAtt                                                                                                                                           <  1:16999/141‑1 (MQ=255)           aTCCAAAACGTTACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACCTCATTCAAACAAACCAATTACAAATCCAAACCTGATGCTATTGGCTATGAATTACAAGGTTCg                                                                                                                                 <  2:459218/141‑1 (MQ=255)             ccAAAACGTTACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGac                                                                                                                               >  2:73475/1‑141 (MQ=255)                     ttACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTgtcg                                                                                                                       >  1:284787/1‑141 (MQ=255)                     ttACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTgtcg                                                                                                                       >  1:337270/1‑141 (MQ=255)                                tCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACCAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTAGACCCCTGTCCTCTCTAAACCt                                                                                                            <  2:211151/141‑1 (MQ=255)                                    ataAAAAATATATTACAAATTCATGGCTCAAGCCATTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAAt                                                                                                        <  2:425339/141‑1 (MQ=255)                                     taataattatatTTCTAATTCATGTCTCAAGCCTTTGCATATGGCACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTATAAACCTGaact                                                                                                       <  2:16712/134‑3 (MQ=255)                                     taCAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAAtt                                                                                                       <  1:93313/141‑1 (MQ=255)                                     taCAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAAtt                                                                                                       <  1:457632/141‑1 (MQ=255)                                       cAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTa                                                                                                     >  2:334442/1‑141 (MQ=255)                                       cAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTa                                                                                                     >  2:424330/1‑141 (MQ=255)                                              tatTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATAACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTg                                                                                              >  1:288599/1‑141 (MQ=255)                                              tatTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATAACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTg                                                                                              >  1:226790/1‑141 (MQ=255)                                                  atagaTTCATGGCTCAAGCATTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGAt                                                                                          <  2:147001/137‑1 (MQ=255)                                                  aCAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGAt                                                                                          <  2:238905/141‑1 (MQ=255)                                                  aCAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGAt                                                                                          <  2:341037/141‑1 (MQ=255)                                                  aCAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGAt                                                                                          <  2:340421/141‑1 (MQ=255)                                                    aaaTTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATac                                                                                        <  1:402035/141‑1 (MQ=255)                                                      agtCATGGCTCAAGCCTTTGCATATTACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATTCTATTGGCTATGAATTACAATGTTCTACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATacac                                                                                      <  2:437540/139‑1 (MQ=255)                                                      aTTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATacac                                                                                      <  2:146012/141‑1 (MQ=255)                                                      aTTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATacac                                                                                      <  2:398975/141‑1 (MQ=255)                                                      aTTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATacac                                                                                      <  1:380878/141‑1 (MQ=255)                                                      aTTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATacac                                                                                      <  1:155934/141‑1 (MQ=255)                                                                                         tcatTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTg                                                   >  1:321595/1‑141 (MQ=255)                                                                                          catTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATAACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTgg                                                  <  2:288599/141‑1 (MQ=255)                                                                                          catTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATAACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTgg                                                  <  2:226790/141‑1 (MQ=255)                                                                                               aaacaaacCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACAt                                             >  1:411165/1‑141 (MQ=255)                                                                                               aaacaaacCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACAt                                             >  1:296688/1‑141 (MQ=255)                                                                                               aaacaaacCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACAt                                             >  1:445774/1‑141 (MQ=255)                                                                                                  caaacCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTgg                                          <  1:390560/141‑1 (MQ=255)                                                                                                  caaacCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGACTATAAACATGCATTTAACCTTTGGGATAAACAATCAGAATATATGATGACATCAACTGAAGATTTGGAACATTgg                                          <  1:67464/141‑1 (MQ=255)                                                                                                      ccAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATaa                                      <  1:186952/141‑1 (MQ=255)                                                                                                      ccAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATaa                                      <  1:440745/141‑1 (MQ=255)                                                                                                       cAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATaaa                                     >  2:166601/1‑141 (MQ=255)                                                                                                                 tCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACtt                           <  1:166601/141‑1 (MQ=255)                                                                                                                      aCATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATcc                      <  2:321595/141‑1 (MQ=255)                                                                                                                         tgatgCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCAcc                   >  2:226175/1‑141 (MQ=255)                                                                                                                         tgatgCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCAcc                   <  2:219123/141‑1 (MQ=255)                                                                                                                          gatgCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGCCAGAATTCCAACTTAATCCACCt                  >  1:214971/1‑141 (MQ=255)                                                                                                                          gatgCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCt                  >  2:271668/1‑141 (MQ=255)                                                                                                                          gatgCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCt                  >  2:6225/1‑141 (MQ=255)                                                                                                                                 ttGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCTGATAAAt           <  1:280166/141‑1 (MQ=255)                                                                                                                                         gAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCTGATAAATACAAATCa   <  1:106428/141‑1 (MQ=255)                                                                                                                                          aaTTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCTGATAAATACAAATCAt  >  1:85053/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                          AGGTTTAATTATCCAAAACGTTACGAAGCAATCCTATACAAAATATATTACAAATTCATGGCTCAAGCCTTTGCATATGACACATACATCATTCAAACAAACCAATTACAAATCAAAACATGATGCTATTGGCTATGAATTACAAGGTTCGACACCTGTCGTCTCTAAACCTGAATTTAACCTTTGGGATACACCATCAGCATATATGATGACATCAACTGAAGATTTGGAACATTGGATAAAATTCCAACTTAATCCACCTGATAAATACAAATCAT  >  CP000730/2568745‑2569022

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