Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	225,710	C→G	T266R (ACA→AGA)	murQ →	possible sugar phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	225,710	0	C	G	94.1%	49.2 / ‑2.4	17	T266R (ACA→AGA)	murQ	possible sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (8/8);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GTTAAAGCAACCAATCAAAAACTGATCGACCGTTCAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGACAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCTTAGGAGGGATTTAAATGACCAAAGAA  >  CP000730/225586‑225833                                                                                                                             |                                                                                                                            gTTAAAGCAACCAATCAAAAACTGATCGACCGTTCAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTAtgtgt                                                                                                             >  2:184281/1‑141 (MQ=255)     aaGCAACCAATCAAAAACTGATCGACCGTTCAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCa                                                                                                         <  1:265454/141‑1 (MQ=255)               tCAAAAACTGATCGACCGTTCAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGga                                                                                               >  2:121318/1‑141 (MQ=255)                    aaCTGATCGACCGTTCAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGaagaag                                                                                          <  2:175737/141‑1 (MQ=255)                                   cAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGACAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGttat                                                                           <  1:297035/141‑1 (MQ=255)                                                    gAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACAtt                                                          <  2:171308/141‑1 (MQ=255)                                                     aaaTATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGATTATTAAACAATGGTGACATTg                                                         >  1:253000/1‑141 (MQ=255)                                                          tgtgCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTaaa                                                    <  2:90267/141‑1 (MQ=255)                                                                   acaTATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGACGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATc                                           <  1:233623/141‑1 (MQ=255)                                                                      tatGATGAAGCAATGGCGATATATCAGGTATGGGCGCGTTATGTGAAAGTTGTGAGAGTTATGGGTATGTGCAGCATGTCTAAGGAATAAGAAACAAGACGGGTAATAAACATAGTGGACATTGTTAAACGATCGATCata                                        >  2:91882/1‑139 (MQ=255)                                                                           tgaAGCCATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAg                                   <  2:66453/141‑1 (MQ=255)                                                                                  aTGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCt                            <  2:218384/141‑1 (MQ=255)                                                                                     gCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCTTag                         >  2:77562/1‑141 (MQ=255)                                                                                              cAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCTTaggtttatggg                >  1:43143/1‑133 (MQ=255)                                                                                              cAGGTATCTGAGCATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCTTAGGAGGGAttt                >  1:171515/1‑141 (MQ=255)                                                                                                          cATGATGTGAAAGTTGCGAGAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCTTAGGAGGGATTTAAATGACCAAAg    >  1:275503/1‑141 (MQ=255)                                                                                                            tgatgTGAAAGTTGCGAGAGTTATGGGTATGTGTGGGATTTCTAAGGAAGAAGCAACTAGATGGTTATTAAACAATGGTCACATTGTTAAACGAGCAGTCAGAGATAGACAACCTTAGGAGAGATTTAAATGACCTAAGaa  <  2:149070/141‑1 (MQ=255)                                                                                                                             |                                                                                                                            GTTAAAGCAACCAATCAAAAACTGATCGACCGTTCAGTGCGTATTATTCAAGAAATATGTGCTATCACATATGATGAAGCAATGGCGTTATATCAGGTATCTGAGCATGATGTGAAAGTTGCGACAGTTATGGGTATGTGTGGCATTTCTAAGGAAGAAGCAACAAGACGGTTATTAAACAATGGTGACATTGTTAAACGAGCAATCAGAGATAGACAACCTTAGGAGGGATTTAAATGACCAAAGAA  >  CP000730/225586‑225833

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