Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	849,263	C→A	A130E (GCA→GAA)	yvcJ →	ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	849,263	0	C	A	95.7%	68.3 / ‑3.6	23	A130E (GCA→GAA)	yvcJ	ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (15/7);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATTCGTCGATACTATGAAGATGAAGAGTTTGAAACTTT  >  CP000730/849132‑849398                                                                                                                                    |                                                                                                                                        cAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAatga                                                                                                                                >  2:199575/1‑141 (MQ=255)       tgaaagtgaCGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGCGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAatga                                                                                                                                <  2:115441/135‑1 (MQ=255)       tgaaagtgaCGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGCGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAatga                                                                                                                                >  1:115441/1‑135 (MQ=255)                     tcatTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGtctc                                                                                                            <  1:296783/141‑1 (MQ=255)                      catTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCa                                                                                                           >  1:232128/1‑141 (MQ=255)                                atgttTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGt                                                                                                 >  2:226934/1‑141 (MQ=255)                                   tttttAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGtata                                                                                              >  1:176011/1‑141 (MQ=255)                                   tttttAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGtata                                                                                              >  2:236355/1‑141 (MQ=255)                                    ttttAGAAGCAAGTACTGAAAAATTAAGTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAg                                                                                             >  2:228873/1‑141 (MQ=255)                                      ttAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCt                                                                                           <  1:236355/141‑1 (MQ=255)                                      ttAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCt                                                                                           <  2:191271/141‑1 (MQ=255)                                              aaGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGt                                                                                   <  1:58105/141‑1 (MQ=255)                                                       aaaaTTAATTTCAAGATATAAGGAAACGTGTTGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAATTTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGAtaa                                                                          >  2:208748/1‑140 (MQ=255)                                                          attgattTTAATATTTAAGAAAACACGTAGTCCACATCCGTTGATGCAACAAGCTAAAAGATTGTTAATCAACGAAATTAATGAGGTGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGAtactac                                                                       >  2:42411/1‑141 (MQ=255)                                                             aattTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACaaa                                                                    >  1:121769/1‑141 (MQ=255)                                                                    aGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCa                                                             >  2:260570/1‑141 (MQ=255)                                                                              aaaCGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAAt                                                   >  2:207829/1‑141 (MQ=255)                                                                                      cgtGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGaa                                           >  2:139958/1‑141 (MQ=255)                                                                                               ccTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATtcgt                                  <  1:124887/141‑1 (MQ=255)                                                                                                     aTGGAACAAGGTAAAAGATCGCTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATTCGTCGATAc                            >  2:211063/1‑141 (MQ=255)                                                                                                          aCAAGGTAAAAGATCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATTCGTCGATACTatga                       <  1:138666/141‑1 (MQ=255)                                                                                                                  aaaGACCGTTAATCAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATTCGTCGATACTaagaagatgaag               >  2:18602/1‑141 (MQ=255)                                                                                                                               cAATGAAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATTCGTCGATACTATGAAGATGAAGAGTTTGAAACttt  <  2:232128/141‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        CAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTATAGATACTACAAAGTTATCACCTAAAGAATTAAAAGAACGCATTCGTCGATACTATGAAGATGAAGAGTTTGAAACTTT  >  CP000730/849132‑849398

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