Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	849,203	A→G	Y110C (TAT→TGT)	yvcJ →	ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	849,203	0	A	G	90.3%	76.1 / 0.1	31	Y110C (TAT→TGT)	yvcJ	ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (19/9);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 4.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

GCAATTGCAATTGATTTAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTAT  >  CP000730/849067‑849320                                                                                                                                         |                                                                                                                      gcaattgcaattGATTTAAGAGGAAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAg                                                                                                                   >  2:393217/1‑141 (MQ=255)            tgATTTAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGcgtcg                                                                                                        >  2:273490/1‑141 (MQ=255)               tttAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGc                                                                                                     <  1:409065/141‑1 (MQ=255)                 tAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGcac                                                                                                   <  2:137984/141‑1 (MQ=255)                   agagGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACAt                                                                                                 >  1:208362/1‑141 (MQ=255)                   agagGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACAt                                                                                                 >  1:323382/1‑141 (MQ=255)                      ggTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCt                                                                                              <  1:321747/141‑1 (MQ=255)                                tATTTAATTCATTAGTTGCAGTAGTGGATAAAGTTAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAAGAACCGCGTCGTGCACATCCTTTGATGGAAc                                                                                    <  1:373056/141‑1 (MQ=255)                                tATTTAATTCATTAGTTGCAGTAGTGGATAAAGTTAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAAGAAACGCGTCGTGCACATCCTTTGATGGAAc                                                                                    <  1:373059/141‑1 (MQ=255)                                tATTTAATTCATTAGTTGCAGTAGTGGATAAAGTTAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAAGAAACGCGTCGTGCACATCCTTTGATGGAAc                                                                                    <  2:32831/141‑1 (MQ=255)                                  tttAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACaa                                                                                  >  2:193619/1‑141 (MQ=255)                                     aaTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGt                                                                               <  2:297434/141‑1 (MQ=255)                                      aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  2:103356/1‑141 (MQ=255)                                      aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  1:5254/1‑141 (MQ=255)                                      aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  1:98777/1‑141 (MQ=255)                                      aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  2:272650/1‑141 (MQ=255)                                                tGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGtt                                                                    >  2:157764/1‑141 (MQ=255)                                                  cAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTaa                                                                  <  1:193560/141‑1 (MQ=255)                                                                gTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAat                                                    <  2:393529/141‑1 (MQ=255)                                                                   aaaaGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAatgat                                                 >  2:25031/1‑141 (MQ=255)                                                                            agtgaCGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGagag                                        >  1:150812/1‑141 (MQ=255)                                                                            agtgaCGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGagag                                        >  1:43514/1‑141 (MQ=255)                                                                                         ttGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCaaa                           <  2:97594/141‑1 (MQ=255)                                                                                                      tttAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAgctatacttc              <  1:193619/141‑10 (MQ=255)                                                                                                            aGCAAGTACTGAAAAATTGATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTATCAAATTAGAAGTATAGCTAAttt        <  1:256762/141‑1 (MQ=255)                                                                                                             gCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAAtttt       >  2:159761/1‑141 (MQ=255)                                                                                                             gCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAAtttt       >  2:126248/1‑141 (MQ=255)                                                                                                             gCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAAtatt       >  2:52460/1‑141 (MQ=255)                                                                                                              cAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTg      >  1:129711/1‑141 (MQ=255)                                                                                                              cAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAACGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTg      >  1:184605/1‑141 (MQ=255)                                                                                                                  tACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTtat  >  2:124801/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                      GCAATTGCAATTGATTTAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTAT  >  CP000730/849067‑849320

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