Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,394,838	A→C	7.3%	K86Q (AAA→CAA)	U712_12190 →	Uncharacterized protein yqzC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,394,838	0	A	C	7.3%	170.4 / 2.6	54	K86Q (AAA→CAA)	U712_12190	Uncharacterized protein yqzC
Reads supporting (aligned to +/- strand):  ref base A (21/29);  new base C (2/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACAGTAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCATTGATGCCGGCTAT  >  CP006881/2394700‑2394966                                                                                                                                           |                                                                                                                                 aCAGTAACCTAACAATACGTAAACAATAATCTTGATACTAAAAAATTGGATACTGTCAACCGAGATCAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACa                                                                                                                                  <  1:377108/139‑1 (MQ=255) aCAGTAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACa                                                                                                                                  <  7:203947/139‑1 (MQ=255) aCAGTAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACa                                                                                                                                  <  7:13328/139‑1 (MQ=255)  cAGTAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACa                                                                                                                                  >  4:393934/1‑138 (MQ=255)   aGTAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACa                                                                                                                                  >  8:188678/1‑137 (MQ=255)     tAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGATAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGt                                                                                                                              >  3:295741/1‑139 (MQ=255)     tAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGt                                                                                                                              >  1:320955/1‑139 (MQ=255)       aCCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTc                                                                                                                             >  2:312089/1‑138 (MQ=255)       aCCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCa                                                                                                                            <  3:393934/139‑1 (MQ=255)       aCCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACAGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTc                                                                                                                             >  3:264835/1‑138 (MQ=255)        cgtaaCAATACGTCAAAAATTATCTTTATTCTAAAAAATTGTTTGCTGTCACCCGTTTTTAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCaa                                                                                                                           <  1:133983/136‑1 (MQ=255)        ccGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCaa                                                                                                                           <  2:266509/139‑1 (MQ=255)        ccGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCaa                                                                                                                           <  5:76229/139‑1 (MQ=255)                  cggtaaCATTTATATTTATTATAAAAAATTGGTTTCTGTCAACGTTTATTAATATAAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAg                                                                                                                 <  2:336185/135‑1 (MQ=255)                    tcaaaaaTTATCTTGATTAAAAAAATTTGGTTTTTGCAACCCGTTATAAAAAAAAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  3:16877/134‑1 (MQ=255)                    tGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  7:66573/139‑1 (MQ=255)                     gAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  7:188678/138‑1 (MQ=255)                     gAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  1:312089/138‑1 (MQ=255)                      aaCAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               >  7:81709/1‑137 (MQ=255)                      aaCAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  8:81709/137‑1 (MQ=255)                           ttATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAACCCTACAAGCTCACCag                                                                                                         >  8:361084/1‑137 (MQ=255)                              tatttaTTCAAAAAAATTGGTTTCTGTCAACCGTTACTAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGCCAGCGGCGCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTaaa                                                                                                     <  4:347034/134‑1 (MQ=255)                               atttaTTATAAAAAATTGTTTTCTTTAAACCTTTATTAAAATCAAAAACTTCTCGATTCAAAATATAAATCATTAAAAAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGATCACCATTaaa                                                                                                     <  5:128449/134‑1 (MQ=255)                                   aTTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACgg                                                                                                >  5:184723/1‑139 (MQ=255)                                             aTTGGTTTCTGCCACCCTTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCAcg                                                                                      <  3:71648/139‑1 (MQ=255)                                             aTTGGTTGCTGCCACCCGGCATCAATATAAAAAACTTCTCGATTCAAAAGACAAATCCTTAAACAATGACCGCGGCTCAGACACGAAAACAGACAAAGTCAAAACCTACACGCTCACCATTAAAGACGGCATGAGCAcg                                                                                      <  4:318140/139‑1 (MQ=255)                                                ggTTTCTGACAACCTTTATGAAAATAAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCAcggc                                                                                    <  4:354728/138‑1 (MQ=255)                                                ggTTTCTGAAACCCGTTTTTAATAAAAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCAcggcg                                                                                   <  4:401150/139‑1 (MQ=255)                                                       gTCAACCGTTATGAAAAAAAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATTATAAAATTTCAAAAATTAAAAGATCACCATTAAAGACGGCATGAGCACGGCGGATGTa                                                                             <  2:320955/138‑1 (MQ=255)                                                                gattaATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTAt                                                                     <  5:94239/137‑1 (MQ=255)                                                                      tatCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGaa                                                              >  5:267010/1‑138 (MQ=255)                                                                         cAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  1:177645/86‑1 (MQ=255)                                                                         cAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               >  2:177645/1‑86 (MQ=255)                                                                                  ctcGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTa                                                 >  2:376243/1‑139 (MQ=255)                                                                                  ctcGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTa                                                 >  4:353468/1‑139 (MQ=255)                                                                                         cAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCtt                                           <  6:184723/138‑1 (MQ=255)                                                                                             agagAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCt                                      <  2:301207/139‑1 (MQ=255)                                                                                                 aaaTCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCaa                                                                                                                           <  2:353951/50‑1 (MQ=255)                                                                                                 aaaTCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCaa                                                                                                                           >  1:353951/1‑50 (MQ=255)                                                                                                 aaaTCCTTAAAAAATGACAGCTGCTCAGACCCGAAATCAGACAAAGTCAAACCCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTAGCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGcgc                                     >  4:161349/1‑134 (MQ=255)                                                                                                           aCAATGACCGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAACGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGa                        <  1:303250/139‑1 (MQ=255)                                                                                                           aCAATCACCGCGGCCCAGACACGAAATCACACAAAGTAAAAACCTACACGCTCACAATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGa                        <  6:267010/139‑1 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               >  7:164169/1‑48 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCt                                                                                                               <  8:164169/48‑1 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTg                                                                         >  3:109473/1‑86 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTg                                                                         <  4:109473/86‑1 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACtt                              >  7:28645/1‑129 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACtt                              <  8:28645/129‑1 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTa                            >  7:245916/1‑131 (MQ=255)                                                                                                               tGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTa                            <  8:245916/131‑1 (MQ=255)                                                                                                                     tggCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCAtt              >  2:13490/2‑139 (MQ=255)                                                                                                                     cGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCAtt              >  6:177807/1‑139 (MQ=255)                                                                                                                             acacGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCATTGATGCCgg      <  3:353468/139‑1 (MQ=255)                                                                                                                                  aaaTCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCATTGATGCCGGCTAt  >  6:145884/1‑138 (MQ=255)                                                                                                                                  aaaTCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCATTGATGCCGGCTAt  >  7:188900/1‑138 (MQ=255)                                                                                                                                     tCAGACCAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAag                                                          >  7:52883/1‑79 (MQ=255)                                                                                                                                     tCAGACCAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAag                                                          >  7:52880/1‑79 (MQ=255)                                                                                                                                     tCAGACCAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAag                                                          <  8:52880/79‑1 (MQ=255)                                                                                                                                     tCAGACCAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAag                                                          <  8:52883/79‑1 (MQ=255)                                                                                                                                           |                                                                                                                                 ACAGTAACCGAACAAGACGTGAACAATTATCTTGATTCTAAAAAATTGGTTTCTGTCAACCGTGATGAATATCAAAAACTTCTCGATTCAAAAGAGAAATCCTTAAACAATGACAGCGGCTCAGACACGAAATCAGACAAAGTCAAAACCTACAAGCTCACCATTAAAGACGGCATGAGCACGGCGGATGTATCTGCTATTCTTGAAAAAGAGGGCATTATCTCTTCGGCTCAAGACTTTAATGACTATGTCATTGATGCCGGCTAT  >  CP006881/2394700‑2394966

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