Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,098,087	A→C	24.1%	A34A (GCT→GCG)	U712_10435 ←	Hypothetical Protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,098,087	0	A	C	24.1%	94.0 / 8.8	59	A34A (GCT→GCG)	U712_10435	Hypothetical Protein
Reads supporting (aligned to +/- strand):  ref base A (34/10);  new base C (2/13);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 1.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTACCATATTGTATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGTTTTCATCCTTTCTTATCACGATATACG  >  CP006881/2097964‑2098220                                                                                                                            |                                                                                                                                       cTTTTACCATTTTGTAAGACTGTCCAAAAGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAATGCACTATTGAAATTACAGCATCGAGCAAAAAAAA‑GCGGCTCCAACaat                                                                                                                         >  1:198524/1‑138 (MQ=255) cTTTTACCATATTGTATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACaat                                                                                                                         >  1:342213/1‑138 (MQ=255)   tttACCATATTGTATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACaataat                                                                                                                      >  6:277335/1‑139 (MQ=255)          tatTGTATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATtaataa                                                                                                                >  2:94750/1‑138 (MQ=255)           atTGTATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATaaaa                                                                                                              >  5:193553/1‑139 (MQ=255)            ttGTAAGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGTAAGGCTCTTTTGAATTTTCAGCATCTAGCAAATAAAC‑GCGGCTCCAACAATAATTAATAAAAt                                                                                                             >  4:309519/1‑139 (MQ=255)               tATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGaa                                                                                                          >  1:352269/1‑139 (MQ=255)               tATGAGTGTCCAAATGGCGGCATGACGTTCCACGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGaa                                                                                                          >  6:69153/1‑139 (MQ=255)                 tGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAAc                                                                                                         >  8:307993/1‑138 (MQ=255)                   aGTGTCAAAATGGCGTCATGACGTACCCCGAAAAAACTTTTTCAAATAAAAAGAGCATTCATAATAAGTATGGCTATTTTTAAATTTCAGAATCTCCAAAATAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACaa                                                                                                       <  7:307993/138‑1 (MQ=255)                    gtgtCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGa                                                                                                     >  7:322958/1‑139 (MQ=255)                      gtCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAAAAATTAATAAAATGAACAAGac                                                                                                    >  1:79452/1‑138 (MQ=255)                       tcaaattGGCGGCATCACGTACCCCCAAAAAAATTTTTAAAAAAAAAAGAGCCTTCCTAATAAGCAAGGCTCTTTTTAAATTTCAGCATCTACAAAATAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGaca                                                                                                   <  1:249153/132‑1 (MQ=255)                       tCCAAATGGCGGGATGACGTTGCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATATAGCAATTAACA‑GCGGCTCCAACAATTATTAATAAAATGAACATGacaa                                                                                                  >  5:157920/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAATTTTCAGCATCTAGCAAATATAA‑GCGGCTCCAACAAAAATTAATAAAATGAACAAGccac                                                                                                  >  8:160323/1‑138 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCTGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  6:232993/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGccaa                                                                                                  >  2:285344/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  6:108356/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  6:183556/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  8:316490/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  6:405802/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  7:329334/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGacaa                                                                                                  >  7:70332/1‑139 (MQ=255)                       tCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGaca                                                                                                   >  1:323984/1‑138 (MQ=255)                               gcggcATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGGTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCTGCTCGAACAATATTTAATAAAATGATCAAGACAACAATTac                                                                                           >  6:229389/1‑137 (MQ=255)                                 ggcATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGACCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAaca                                                                                         >  7:270962/1‑138 (MQ=255)                                           tCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTAttt                                                                              >  6:108561/1‑139 (MQ=255)                                               cGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTcc                                                                          >  3:277117/1‑139 (MQ=255)                                               cGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTcc                                                                          >  4:9016/1‑139 (MQ=255)                                                  aaaaaCTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACACTTAACACAAATGTATTTGTCCCg                                                                        >  5:45631/1‑138 (MQ=255)                                                      aatttttAAAATAAAACGACCAATCATAATAAGGAAGGCTTTTTTTAAATTTCAGCATAAACAAAAAAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGcc                                                                   <  1:94750/137‑1 (MQ=255)                                                      aCTTTTTCACATAAAAAGAGCCTTCATAAAACGGATGGCTCTTTTTAAATTTCAGCATCTAGAAAATAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGcc                                                                   <  6:193553/139‑1 (MQ=255)                                                      aCTTTTTCACATAAAAACAGCCTTCATAATAAGGAAGGCTATTTTTAAATTTCAGCATATAGCAAAAAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGcc                                                                   <  7:202544/139‑1 (MQ=255)                                                      aCTTTTTCAAAAAAAAAGACCCTACATAAAAAGTATGGCTATTTTTAAATTTCAGCATCTACAAAAAAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGcc                                                                   <  1:99982/139‑1 (MQ=255)                                                      aCTTTATCACATAAAACCACCCTGCATAATATGGCCGGATTTTTTAAAATTCCCGCATAAACCAAAAAACC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGcc                                                                   <  5:5849/139‑1 (MQ=255)                                                      aCTTATACACAAAAACAGAGCCTCCAAAAAACGGATGGCTCTTTTGAAATTTCCCACACACGCAAAAAACC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGcc                                                                   <  4:255475/139‑1 (MQ=255)                                                        tttttAAAAAAAAAACAGCATCCCAAATAGGTATGGCTCTTTTTAAATTTCCGAAACTTGAAACAAAAC‑GCGGCTCAAACATTAATTAATAAAATGAACAAGACAACCATTAACACAAATGTATTTGTCCCGCAGcgat                                                                 <  6:357044/139‑4 (MQ=255)                                                             aaaaTAAAAAGCCCCTTCATAATATGGATGGTTTTTTTTAAATTTCAGAATATACAAAAAAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCg                                                            <  6:365199/136‑1 (MQ=255)                                                                 tAAAAATAGCATTCCTAATAAGGAAGGATATTTTTAAATTTCAGCATATAGAAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAAc                                                        <  8:322958/139‑1 (MQ=255)                                                                 tAAAAAGAGCCTTCATAATAAGTATGGCTTTTTTTAAATTTCAGAATCAAGAAAAAAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAAc                                                        <  4:95257/139‑1 (MQ=255)                                                                  aaaaaGAGCCTTCATAAAAGGTACGGCTCTTTTTAAATTTCAGCATCTAGCAAATAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAAcc                                                       <  8:330546/139‑1 (MQ=255)                                                                          ccTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACTCAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCag                                               >  6:317909/1‑139 (MQ=255)                                                                          ccTGCCTAATAAGGAAGGCTCTTTTGAAATGTCCGCATCGAGCAAATAACA‑GCGGCGCCAACAATAATTAATACCATGAACAAGACAACAATGAACACAAATGTATTTGTCCCGCAGCCAACTCCGTAACCGTACCCag                                               >  1:238/1‑139 (MQ=255)                                                                           ctccataaaaaGGATGGTTTTTTTTAATTTTCATCATCTAGAAAATAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCaga                                              <  7:229005/134‑1 (MQ=255)                                                                           ctaacTAATATGTAAGGCTTTTTTTAAATTACAGCATCTACAAAAAAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCaga                                              <  4:310833/135‑1 (MQ=255)                                                                             tCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATTAAATGAACAAGACAACAATTAACACAAATGTATTTGTACCGCAGCCAAATCCGTAACCTTACCCagag                                             >  1:122312/1‑138 (MQ=255)                                                                              cataataaTGATGGCTATTTTTAAATTTCACCATCTAGAAAAAAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCagag                                             <  4:277117/135‑1 (MQ=255)                                                                              cataataaGTATGGCTATTTTTAATTTTCACAATAAAACAAAAAAAC‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTa                                           <  2:325131/137‑1 (MQ=255)                                                                              cataataaGTAAGGCTCTTTTTAAATTTCACATTCTAGAAAAAAACA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTa                                           <  5:77876/137‑1 (MQ=255)                                                                                    aaggaagcCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAATTGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTAACCATa                                     >  4:187014/1‑139 (MQ=255)                                                                                       gaaggCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACATGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACAt                                  >  2:78874/1‑139 (MQ=255)                                                                                         tttgttttttTAAATTTCAGATTTTAGAAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATg                                 <  8:70332/132‑1 (MQ=255)                                                                                         cggCTATTTTTAAATTTCAGC‑TAAAACAAAAAAAACGCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATg                                 <  7:316490/137‑1 (MQ=255)                                                                                               ttttGAAATTTCAGCATCTAGCAAATAAAA‑GCTGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGtttt                          >  4:78955/1‑139 (MQ=255)                                                                                                                ctacaaaaaaaCA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGTTTTCATCCTTTCTTATCACg         <  5:229389/134‑1 (MQ=255)                                                                                                                aaaaaaaaaaaaa‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGTTTTCATCCTTTCTTATCACg         <  2:352269/134‑1 (MQ=255)                                                                                                                  aGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGTTTTCATCCTTTCTTATCACGa        >  6:232592/1‑138 (MQ=255)                                                                                                                       aaaaaa‑GATGCTCACACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCAGCAGCCAAATCAGTAACAGTACCCAGAGTATCCATACATTCTTTTTTCATCATTTCTTATCACGATATACg  <  5:372732/137‑1 (MQ=255)                                                                                                                       aTAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCag                                               <  6:184941/94‑1 (MQ=255)                                                                                                                       aTAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCag                                               >  5:184941/1‑94 (MQ=255)                                                                                                                       aTAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGTTTTCATCCTTTCTTATCACGATATACg  <  2:323984/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                       CTTTTACCATATTGTATGAGTGTCCAAATGGCGGCATGACGTTCCCCGCAAAAACTTTTTCACATAAAAAGAGCCTTCCTAATAAGGAAGGCTCTTTTGAAATTTCAGCATCTAGCAAATAAAA‑GCGGCTCCAACAATAATTAATAAAATGAACAAGACAACAATTAACACAAATGTATTTGTCCCGCAGCCAAATCCGTAACCGTACCCAGAGTATCCATACATGCTGTTTTCATCCTTTCTTATCACGATATACG  >  CP006881/2097964‑2098220

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