Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,274,965	G→C	45.0%	E106D (GAG→GAC)	bglG →	transcriptional antiterminator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,274,965	0	G	C	45.0%	7.8 / 23.3	20	E106D (GAG→GAC)	bglG	transcriptional antiterminator
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base C (4/5);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAGGGCTACGCATTGCTGGCACAGATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGAAAGCGTGGCGAAGGCATTTA  >  CP000730/2274829‑2275099                                                                                                                                         |                                                                                                                                       aaaaaGGGCTACGCATTGCTGGCACAGATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAg                                                                                                                                    >  2:36586/1‑141 (MQ=255)            cGCATGGCTGGCACAGATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTATACAatat                                                                                                                         <  2:121251/141‑1 (MQ=255)                        acaGATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCaa                                                                                                             >  1:95505/1‑141 (MQ=255)                           gATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGaa                                                                                                          >  1:99993/1‑141 (MQ=255)                           gATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGaa                                                                                                          >  1:243720/1‑141 (MQ=255)                                      aaaCGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGCTTc                                                                                               >  2:96495/1‑141 (MQ=255)                                          gATTTGAAGCTATCGTGTGCGCGACATCAAACCATTAACTTATCTGTTGAAGTGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTc                                                                                           <  2:153244/141‑1 (MQ=255)                                                                  caTCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATTTTAgat                                                                   >  1:63883/1‑141 (MQ=255)                                                                            ttGAATTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCTAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGc                                                         >  1:190886/1‑141 (MQ=255)                                                                                    tCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTa                                                 >  1:237328/1‑141 (MQ=255)                                                                                       gTTGAAGAGCAGAAAGTAATTATTATATACGGTGTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTGAGATGATTTAGAGCTTGATTTAAAt                                              <  2:9559/141‑1 (MQ=255)                                                                                           aaGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGt                                          <  2:167513/141‑1 (MQ=255)                                                                                                   aaaGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACAAACTa                                  >  1:12220/1‑141 (MQ=255)                                                                                                        aaTTATTATATACGCTTTGATACAAGCCAAGGACCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCttt                             <  1:114546/141‑1 (MQ=255)                                                                                                               atacaCGCTTTGATACAAGCCAAGGACCCCGTTAAACAACAGAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGa                      <  1:144902/141‑1 (MQ=255)                                                                                                                    cGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGAAAGCg                 >  2:103560/1‑141 (MQ=255)                                                                                                                          gATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGAAAGCGTGGCGa           >  1:257270/1‑141 (MQ=255)                                                                                                                            tACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGAAAGCGTGGCGAAg         >  1:75669/1‑141 (MQ=255)                                                                                                                              caggCCAAGGAGCCAGTTGAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCGACTATCTTTATTTCGAAAGCGTGGCGAAGGc       >  2:259854/4‑141 (MQ=255)                                                                                                                                   cAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGAAAGCGTGGCGAAGGCATTTa  <  1:240976/141‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       AAAAAGGGCTACGCATTGCTGGCACAGATTCTGATTTAAACGATTTGAAGCAATCGATTGCACAACATCAAACCATTGACTTATCTGTTGAAGAGCAGAAAGTAATTATTATATACGCTTTGATACAAGCCAAGGAGCCAGTTAAACAATATAGTTTAGCGCAAGAAATTGGCGTTTCTGTCCAGACTTTAGCAAAGATGTTAGATGATTTAGAGCTTGATTTAAATAAGTACCAACTATCTTTATCTCGAAAGCGTGGCGAAGGCATTTA  >  CP000730/2274829‑2275099

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