Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,163,829	T→C	V72A (GTG→GCG)	ysnB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,163,829	0	T	C	97.1%	120.4 / NA	35	V72A (GTG→GCG)	ysnB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/17);  minor base G (1/0);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGATGCAGATGTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGTGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACATGTGGCAAAATATGAGTATATTAATGGTG  >  CP000730/1163686‑1163969                                                                                                                                                |                                                                                                                                             cccAGATGCAGATGTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGaaa                                                                                                                                          >  1:77158/1‑148 (MQ=255)  ccAGATGCAGATGTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAg                                                                                                                                         >  1:141503/1‑148 (MQ=255)  ccAGATGCAGATGTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAg                                                                                                                                         <  2:197531/148‑1 (MQ=255)     gatgcagatgTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCAt                                                                                                                                      <  1:102999/148‑1 (MQ=255)     gatgcagatgTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCAt                                                                                                                                      <  2:273820/148‑1 (MQ=255)     aatgcagatgTATATTTACATTTATTAGATTCAGAATTCGCGTATGAGGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCAt                                                                                                                                      <  2:36414/147‑1 (MQ=255)                      aCATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACAt                                                                                                                     >  1:47611/1‑148 (MQ=255)                         tttAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTta                                                                                                                  >  1:361455/1‑148 (MQ=255)                         tttAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTta                                                                                                                  >  2:266548/1‑148 (MQ=255)                             ggAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATc                                                                                                              >  2:374480/1‑148 (MQ=255)                               agaTTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCaa                                                                                                            >  2:383080/1‑148 (MQ=255)                                gaTTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAg                                                                                                           >  2:304103/1‑148 (MQ=255)                                gaTTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAg                                                                                                           >  1:299419/1‑148 (MQ=255)                                gaTTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAg                                                                                                           >  1:299417/1‑148 (MQ=255)                                            gcgTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACaa                                                                                               <  1:291430/148‑1 (MQ=255)                                              gTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAaga                                                                                             >  2:387367/1‑148 (MQ=255)                                              gTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAaga                                                                                             >  2:121366/1‑148 (MQ=255)                                              gTATGATGATACGGAACTTAGCTTATTAAATAGAGTAAAAGGCAATTGTGAATTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCACTTTATACTCATGGACATTTATATCAAGTCAATCGAACaaca                                                                                             >  2:147021/1‑146 (MQ=255)                                                                  gCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAgg                                                                         <  2:71527/148‑1 (MQ=255)                                                                            tAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGGAAAAGGACTTGAATTa                                                               >  1:335448/1‑148 (MQ=255)                                                                            tAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTa                                                               >  1:397956/1‑148 (MQ=255)                                                                                gTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGgtt                                                           >  2:112580/1‑148 (MQ=255)                                                                                      ggCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGttgtttgt                                                     <  1:362121/148‑1 (MQ=255)                                                                                      ggCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGttgtttgt                                                     <  1:19985/148‑1 (MQ=255)                                                                                      ggCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGttgtttgt                                                     <  1:234426/148‑1 (MQ=255)                                                                                      ggCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCAGAAAAAGGACTTGAATTAGGttgtttgt                                                     <  2:135111/148‑1 (MQ=255)                                                                                      ggCAATTGTGATTATTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATATAAAAAAAAAAAAAGGACTTGAATTAGGttgtttgt                                                     <  2:160358/148‑1 (MQ=255)                                                                                                   tttACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATgg                                        <  1:101299/148‑1 (MQ=255)                                                                                                   tttACCCAGAATTTGAAAAAGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATgg                                        <  1:101313/148‑1 (MQ=255)                                                                                                          agggttagAAAGTGAAGCGGATGCAAAATGTAATCACGGAGAAACAATCAAGACACTTGAAAAATAATATCAAGTCAATAAAACAAGAAATGTAAGAGCTGAAAAAGGACTTGAATTAGGCTGTTTGGTTGCATTTTCTGTAGATacc                                 >  2:200444/8‑147 (MQ=255)                                                                                                              tttGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACAtg                             <  1:240869/148‑1 (MQ=255)                                                                                                              tttGAAAATGAAGCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACAtg                             <  2:143160/148‑1 (MQ=255)                                                                                                                          gCGGTCGCAAAATATAATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACATGTGGCAAAATATg                 <  1:383080/148‑1 (MQ=255)                                                                                                                                      tataATGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACATGTGGCAAAATATGAGTATATTAATg     >  1:51246/1‑148 (MQ=255)                                                                                                                                         aaTGACGCGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAATAAGAGATAAAAAAAAAAAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACATGTGGCAAAATATGAGTATATTAATGgtg  <  2:51158/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             CCCAGATGCAGATGTATATTTACATTTAGGAGATTCAGAATTCGCGTATGATGATACGGAACTTAGCTTATTTAATAGAGTAAAAGGCAATTGTGATTTTTACCCAGAATTTGAAAATGAAGCGGTCGCAAAATATAATGACGTGAAAGCATTTTATACTCATGGACATTTATATCAAGTCAATCGAACAAGAGATTTATTAGCTGAAAAAGGACTTGAATTAGGTTGTTTGTTTGCATTTTATGGACATACACATGTGGCAAAATATGAGTATATTAATGGTG  >  CP000730/1163686‑1163969

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