Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	66,188	A→T	18.0%	M458L (ATG→TTG)	U712_00295 →	putative membrane protein yabM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	66,188	0	A	T	18.0%	100.3 / 3.9	47	M458L (ATG→TTG)	U712_00295	putative membrane protein yabM
Reads supporting (aligned to +/- strand):  ref base A (10/28);  new base T (7/2);  total (17/30)
Fisher's exact test for biased strand distribution p-value = 6.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGTTTATTTATTGTATGATGAGAT  >  CP006881/66051‑66323                                                                                                                                          |                                                                                                                                         ggATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGAt                                                                                                                                         <  4:151422/139‑1 (MQ=255)   aTTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGt                                                                                                                                       <  3:219030/139‑1 (MQ=255)    ttGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTc                                                                                                                                      <  4:75739/139‑1 (MQ=255)       aaGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTTATTGAGATTATTTGTTGTGCTGTTTTTTTGTCaag                                                                                                                                   >  5:203373/1‑137 (MQ=255)       aaGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGGTTCCGATTATTTGATGTGGTCGTTTGTTGTGAGc                                                                                                                                   >  1:295462/1‑139 (MQ=255)            gcATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTTCGATTATTTGATGTTCTCTTTTGTTGTGAGCTt‑ttt                                                                                                                             >  6:90274/1‑139 (MQ=255)             cATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAATGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCACATTATTTGATTTGCTCTTTTTTTTTCAGCTG‑tttt                                                                                                                            >  6:165241/1‑139 (MQ=255)              aTCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGATTTCGAATTTTGTGTTGGCTGTTTTGATGTCAGCTG‑TTttg                                                                                                                           >  4:285504/1‑139 (MQ=255)               tCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTttgc                                                                                                                          <  5:289040/139‑1 (MQ=255)                 gCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTttgctt                                                                                                                        <  5:234859/139‑1 (MQ=255)                 gCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTttgctt                                                                                                                        <  7:124529/139‑1 (MQ=255)                 gCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTttgct                                                                                                                         <  8:223272/138‑1 (MQ=255)                          aGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGTTTCCTTTTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTTTAcg                                                                                                               <  4:185408/139‑1 (MQ=255)                          aGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATAcg                                                                                                               <  6:53433/139‑1 (MQ=255)                          aGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATAcg                                                                                                               <  5:19393/139‑1 (MQ=255)                          aGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATAcg                                                                                                               <  1:261722/139‑1 (MQ=255)                          aGCGGCGTCCTGTGCGCCTGTCGCCGGCCTTAACTTGTCTCAGCTGCGGCAAAATGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATAcg                                                                                                               <  5:143165/139‑1 (MQ=255)                            cGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTACGATTATTGGATGTGGTGTTTTGATGTCAGCTG‑TTTTGGTTGCGGATAcgcg                                                                                                             >  3:296667/1‑139 (MQ=255)                                             gTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGTTCTGCGCTTTTGATGTGAGCGG‑TTTTGGTTGCGTATTCGCGTCTGGGGACGTtgttgt                                                                                            >  7:110703/1‑139 (MQ=255)                                             gTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGTTTGCGATTATTGGATCTGCTCGTTTGTTGTCAGCTG‑TTTTTGTTTCGTTTTCGCGGCTCTGGACTTtgttgt                                                                                            >  3:302471/1‑139 (MQ=255)                                                     ccTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCa                                                                                    <  7:140092/139‑1 (MQ=255)                                                     ccTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCGCGATAAGCTTAGAGGCGTGTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCa                                                                                    <  6:204337/139‑1 (MQ=255)                                                              gTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGATTGCGATTATTGGATCTGCGGTTTTGATGTCAGCTG‑TTTTGCTTTCGTTTACGCCTGTCTGGACTTTGTGGTTTGCTGCAACAGGaa                                                                             >  2:228734/1‑137 (MQ=255)                                                                    gCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGATTCCGATTATTTGGTCTGGGCTTTTGTTGTCAGGGG‑TTTTTCTTTGGTTTAAGCGTCTCTGGGCTTTATTGTTTGCTGGAACAGGAATGGGAGCg                                                                     >  7:300862/1‑139 (MQ=255)                                                                     cTTCTGCAAAAGTAATGGCTCGATAAGATTAGAGGCGTTTTGATTCCTATTATTGGATCTGCTCTTTTGATCTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCgg                                                                    <  6:107522/139‑1 (MQ=255)                                                                        cGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCt                                                                  <  4:302471/138‑1 (MQ=255)                                                                                     ggCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCa                                                                                    >  7:57528/1‑107 (MQ=255)                                                                                     ggCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCa                                                                                    <  8:57528/107‑1 (MQ=255)                                                                                         cGATAAGCTTAGAGGCGTGTTTTTTACGATTATTTGATATGCTCTTTTTTTTTAAACTGTTTTTGTTTGC‑TTTAAGCGTCTCTGGGCTTTGTTGTTTTCTGCAAAAGGAATACGCGCGGCGGTGATTTAGAGtaatttt                                                >  6:109558/1‑133 (MQ=255)                                                                                          gATAAGCTTAGAGTCGTTTTTATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCTTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTc                                                 <  3:282777/137‑1 (MQ=255)                                                                                              aGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATCCGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTctgctg                                            <  8:300862/138‑1 (MQ=255)                                                                                              aGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTctgctg                                            <  6:126284/138‑1 (MQ=255)                                                                                                             tGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGgcg                             <  8:110703/138‑1 (MQ=255)                                                                                                              gATTGCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGgcgc                            <  3:200631/138‑1 (MQ=255)                                                                                                              gATTCCGATTATTGGTTCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTTGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGgcgc                            <  5:325981/138‑1 (MQ=255)                                                                                                              gATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTTGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGgcgc                            <  3:122384/138‑1 (MQ=255)                                                                                                               aTTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGcgcg                           <  3:285504/138‑1 (MQ=255)                                                                                                               aTTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGcgcg                           <  3:24796/138‑1 (MQ=255)                                                                                                                     attattGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGgtta                    >  1:252020/1‑139 (MQ=255)                                                                                                                         ttGGATCTGGTGTTTTGTTGTAAGCTG‑TTTTGCTTGGGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGtttattta                <  5:165241/139‑1 (MQ=255)                                                                                                                         ttGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGgttattta                >  1:12424/1‑139 (MQ=255)                                                                                                                          agtaTCTGGTCTTTTGTTGGCAGCTG‑TTTTGGTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGTttatttat               <  4:183682/136‑1 (MQ=255)                                                                                                                               cTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGcgcg                           >  7:49842/1‑122 (MQ=255)                                                                                                                               cTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGcgcg                           <  8:49842/122‑1 (MQ=255)                                                                                                                                 gCTCTTTTGATGTATGGTG‑TTTTGCTTGAGTATACGCGTCTCTGGACTTTGTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGTTTATTTATTGTatga        <  8:230486/139‑1 (MQ=255)                                                                                                                                  ctctTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGTTTATTTATTGTatgat       >  7:45143/1‑139 (MQ=255)                                                                                                                                       ttGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGTTTATTTTTTGTTTTATgagct  >  3:183682/1‑137 (MQ=255)                                                                                                                                          |                                                                                                                                         GGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTG‑TTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGGCGCGGTGTTTATTTATTGTATGATGAGAT  >  CP006881/66051‑66323

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