Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	825,618	G→A	6.3%	A31A (GCG→GCA)	U712_04000 →	Uncharacterized protein yfkI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	825,618	0	G	A	6.3%	204.0 / 2.4	63	A31A (GCG→GCA)	U712_04000	Uncharacterized protein yfkI
Reads supporting (aligned to +/- strand):  ref base G (22/37);  new base A (2/2);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-01

CATCCATTCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGAACGCTTTCAAAGGAAGTCGA  >  CP006881/825480‑825752                                                                                                                                           |                                                                                                                                       cATCCATTCTAAAAGAAAAACAATTGTGAGTGATAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCg                                                                                                                                        <  5:34476/139‑1 (MQ=255)   tCCATTCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTg                                                                                                                                      <  1:30646/139‑1 (MQ=255)    ccATTCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGc                                                                                                                                     >  1:20171/1‑139 (MQ=255)    ccATTCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGc                                                                                                                                     >  2:45756/1‑139 (MQ=255)     cATTCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCAGACAAGACAGGCGTGc                                                                                                                                     <  1:45756/138‑1 (MQ=255)       ttCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTg                                                                                                                                      >  8:265188/1‑135 (MQ=255)       ttCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTg                                                                                                                                      <  7:265188/135‑1 (MQ=255)       ttCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCggg                                                                                                                                  >  1:354760/1‑139 (MQ=255)       ttCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCgg                                                                                                                                   >  8:17143/1‑138 (MQ=255)                        ttGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATg                                                                                                                 <  1:357061/139‑1 (MQ=255)                          gtgagtgaAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGc                                                                                                                <  7:17143/138‑1 (MQ=255)                          gtgagtgaAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGc                                                                                                                <  5:215077/138‑1 (MQ=255)                           tgagggaAAGTTGAATATGCATGAACAACGACCGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTTCTCCATAAGCCGACAAGACAGGCGAGCGGGATGAGATGGCTGACATGCa                                                                                                               <  5:131310/138‑1 (MQ=255)                                gcaaGGAGAATATGCATGACCAACGAACGGATAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATa                                                                                                         <  4:330420/137‑1 (MQ=255)                                        aaTTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 <  5:282153/139‑1 (MQ=255)                                        aaTTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 <  1:334779/139‑1 (MQ=255)                                         aTTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTg                                                                                                                                      <  8:145203/101‑1 (MQ=255)                                         aTTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 <  5:277977/138‑1 (MQ=255)                                         aTTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCAc                                                                                                <  1:238973/139‑1 (MQ=255)                                         aTTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGAAATGCGTg                                                                                                                                      >  7:145203/1‑101 (MQ=255)                                                   acaacGAACGTTTAATGCTGTAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      <  5:1616/139‑1 (MQ=255)                                                   acaacGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      <  8:13802/139‑1 (MQ=255)                                                   acaacGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      <  1:170981/139‑1 (MQ=255)                                                   acaacGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGAGGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      <  4:102500/139‑1 (MQ=255)                                                   acaacGAACGGTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      <  2:339946/139‑1 (MQ=255)                                                             tttAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACATTCACTGTACAAAAGCAACCCagag                                                                            <  7:379261/139‑1 (MQ=255)                                                                   gCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGtt                                                                       <  2:354760/138‑1 (MQ=255)                                                                      gAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTaaa                                                                    >  2:76829/1‑138 (MQ=255)                                                                       aaaGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTaaaa                                                                   >  4:397285/1‑138 (MQ=255)                                                                       aaaGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTGATCACAGCTCCATAACCCGACAAGACAGGCGTGCTGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACATAAGCACACCATAGCTGTTaaaa                                                                   >  5:22875/1‑138 (MQ=255)                                                                        aaGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAAcc                                                                                 <  6:282815/123‑1 (MQ=255)                                                                        aaGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAAcc                                                                                 >  5:282815/1‑123 (MQ=255)                                                                          gggATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCGGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAaca                                                               <  8:85395/139‑1 (MQ=255)                                                                              tatTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCACCCCAGAGCTGTTAAAAAACACTGt                                                           >  2:11897/1‑139 (MQ=255)                                                                              tatTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCATGCGGGATGAGATGGCTGACATg                                                                                                                 <  6:26310/85‑1 (MQ=255)                                                                              tatTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCATGCGGGATGAGATGGCTGACATg                                                                                                                 <  6:20760/85‑1 (MQ=255)                                                                              tatTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCATGCGGGATGAGATGGCTGACATg                                                                                                                 >  5:26310/1‑85 (MQ=255)                                                                              tatTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCATGCGGGATGAGATGGCTGACATg                                                                                                                 >  5:20760/1‑85 (MQ=255)                                                                                     cGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGCCAGGCGTGCGGGATGAGATGGCTTACATGCAAGCATAAACTTTCACTGTACAAAAGAAACCCAGAGCTGTTAAAAAACCATGTAATcac                                                     >  6:227254/1‑138 (MQ=255)                                                                                      ggAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGCGGGCTGACATGCAAGCATAAACTTTCACTGTACAACAGCCACCCAGAGCTGTTAAAAAACCCTGTCATcac                                                     <  4:39554/137‑1 (MQ=255)                                                                                      ggAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATg                                                                                                                 >  3:2402/1‑77 (MQ=255)                                                                                      ggAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATg                                                                                                                 <  4:2402/77‑1 (MQ=255)                                                                                      ggAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGtcatca                                                      <  2:20171/136‑1 (MQ=255)                                                                                      ggAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATcac                                                     <  7:201469/137‑1 (MQ=255)                                                                                      ggAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATcac                                                     <  1:48753/137‑1 (MQ=255)                                                                                               gCAGGCGCGGCGTTATCACTGCTCCCTAAGCCGACAAGACCGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATg                                          <  8:398000/139‑1 (MQ=255)                                                                                                    cgcgGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      >  8:356366/1‑90 (MQ=255)                                                                                                    cgcgGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAg                                                                                      <  7:356366/90‑1 (MQ=255)                                                                                                                   gCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTcgcgcg                       >  1:320732/1‑138 (MQ=255)                                                                                                                   gCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCg                       >  7:233322/1‑138 (MQ=255)                                                                                                                    cTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACaaaa                                                 <  6:329531/111‑1 (MQ=255)                                                                                                                    cTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACaaaa                                                 >  5:329531/1‑111 (MQ=255)                                                                                                                    cTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGaa                     <  1:19226/139‑1 (MQ=255)                                                                                                                      ccATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 <  5:334446/61‑1 (MQ=255)                                                                                                                      ccATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 >  6:334446/1‑61 (MQ=255)                                                                                                                         tAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 >  4:128070/1‑58 (MQ=255)                                                                                                                         tAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 <  3:128070/58‑1 (MQ=255)                                                                                                                         tAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGaa                     >  6:311262/1‑134 (MQ=255)                                                                                                                         tAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGaa                     <  5:311262/134‑1 (MQ=255)                                                                                                                             ccGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGAACGCTTTCa             >  8:216158/1‑138 (MQ=255)                                                                                                                              ctaCAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGAACGCTTTCaaa           <  7:125752/137‑1 (MQ=255)                                                                                                                              cGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 <  5:198695/53‑1 (MQ=255)                                                                                                                              cGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCa                                                                                                 >  6:198695/1‑53 (MQ=255)                                                                                                                                       aGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGAACGCTTTCAAAGGAAGTCGa  <  1:76829/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                       CATCCATTCAAAAAGAAAAACAATTGTGAGTGAAAGGAGAATTTGCATGAACAACGAACGTTTAATGCTGAAAGGGATATTTCTCGGAGCTGCGGCAGGCGCGGCGTTATCACTGCTCCATAAGCCGACAAGACAGGCGTGCGGGATGAGATGGCTGACATGCAAGCATAAACTTTCACTGTACAAAAGCAACCCAGAGCTGTTAAAAAACACTGTCATCACAAAAGTGGATGAGGCCAAAAAGCTCGCCCGAACGCTTTCAAAGGAAGTCGA  >  CP006881/825480‑825752

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: ‑