Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,057,150	C→T	A295A (GCG→GCA)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,057,150	0	C	T	95.0%	49.1 / ‑5.9	20	A295A (GCG→GCA)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/13);  minor base G (1/0);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CAGCAAAGCACAGGTACGGGTTAGCTGCCGGATCCGGGAAACGTACTTCGATACGACGTGCTTTCGGAGAAGAAACCACCGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTTGTGCACAACATATTTGTAGATCTGAATTTCGTCAGCTTTTTTG  >  NC_000913/4056915‑4057345                                                                                                                                                                                                                                            |                                                                                                                                                                                                    nagcaAAGCACAGGTACGGGTTAGCTGCCGGATCCGGGAAACGTACTTCGATACGACGTGCTTTCGGAGAAGAAACCACCGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTcagac                                                                                                                                                                                         >  1:359211/2‑248 (MQ=255)                      tAGCTGCCGGATCCGGGAAACGTACTTCGATACGACGTGCTTTCGGAGAAGAAACCACCGGAATACGGATAGACGCAGAACGGGTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGTCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGc                                                                                                                                                                  <  1:172195/250‑1 (MQ=255)                           gCCGGATCCGGGAAACGTACTTCGATACGACGTGCTTTCGGAGAAGAAACCACCGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGaan                                                                                                                                                              <  1:265815/249‑2 (MQ=255)                                      gAAACGTACTTCGATGCGACGTGCTTTCGGAGAAGAAACCACCGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTACGCAAGCCTTACAGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGTCGTTACTCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGt                                                                                                                                                 >  1:197571/1‑251 (MQ=255)                                           gTACTTCGATACGACGTGCTTTCGGAGAAGAAACCGCCGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTGTAGCGAGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACg                                                                                                                                                     <  1:314208/242‑1 (MQ=255)                                                                             caccGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCgn                                                                                                          <  1:271815/251‑2 (MQ=255)                                                                                 ggAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATGACGCCGCCAATGTAGTACAGTGCCTGCTCATACAGACCTGCGTATTTGTCTCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATc                                                                                                      <  1:412219/251‑1 (MQ=255)                                                                                               cGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATACGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCCATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGtt                                                                                        >  1:195428/1‑251 (MQ=255)                                                                                                        gTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACACGAAGCTTATAAGCGTTGGTGGTCGGGTTTGCAAGGGCGTTCATCGCTTTAGCGTGTTTGATTAAGCCGCCACTGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCggcg                                                                                                              >  1:410071/1‑218 (MQ=255)                                                                                                                       aGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTTTAAGAGTTGGTGGTCGGGTTTGCCATGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAAGGTAGTACAGTGCCTGCTCATACAGACCTGCGTATTTGTCTCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCgg                                                                      <  1:428104/245‑1 (MQ=255)                                                                                                                                          tCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCGGCTCAGACAGACCTGCGTATTGGTCGCCTGCGAACAGGTTAACGCCGTTTTGAGACAGAGATATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTcg                                                                         <  1:69485/223‑1 (MQ=255)                                                                                                                                             tAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGGACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGTTTAACGCCGTTTTTAGACATAGACATGGGGCATTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCgg                                                                      <  1:259385/223‑1 (MQ=255)                                                                                                                                                              cTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGACGCCAATGAAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTAcc                                                               >  1:416137/1‑213 (MQ=255)                                                                                                                                                                                  ggTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTa                                                                 >  1:155705/1‑191 (MQ=255)                                                                                                                                                                                      tGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTTGTGCACAACATATTTGTAGATCTGAATTTCGTCAGCTTTTTTg  <  1:344349/250‑1 (MQ=255)                                                                                                                                                                                        ccAGGGCGTTAATCGCTTTAGCGTGTTTGATTACTCCGCCAAGGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGc                                                                                                                                                    <  1:300335/102‑1 (MQ=255)                                                                                                                                                                                                                gtTTGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACa                                                                                                                                                             <  2:11526/69‑1 (MQ=255)                                                                                                                                                                                                                gtTTGATTACGCCGCCAATGTAGTACAGGGCCTGCTCAGACAGACCTGCGTCTTTGTCGCCTGCGAACa                                                                                                                                                             >  1:11526/1‑69 (MQ=255)                                                                                                                                                                                                                 tttGATTACGCCGCCAATGTAGTACAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCg                                                                                           <  1:138844/134‑1 (MQ=255)                                                                                                                                                                                                                                   gtagtaCAGTGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTGTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCAGAAAGGTCTCGGTTTTACCGAAGCGGTGCGCTACGTTGTGCACAACATATTTGTAGATCTGAATTTCGt             <  1:233336/194‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                    CAGCAAAGCACAGGTACGGGTTAGCTGCCGGATCCGGGAAACGTACTTCGATACGACGTGCTTTCGGAGAAGAAACCACCGGAATACGGATAGACGCAGAACGGTTACGCGCAGAGTAAGCCAGCATTACCGGTGCTTCATAGCCCGGGACCAGACGCTTATAAGAGTTGGTGGTCGGGTTTGCCAGGGCGTTAATCGCTTTAGCGTGTTTGATTACGCCGCCAATGTAGTACAGCGCCTGCTCAGACAGACCTGCGTATTTGTCGCCTGCGAACAGGTTAACGCCGTTTTTAGACAGAGACATGTGGCAGTGCATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTTGTGCACAACATATTTGTAGATCTGAATTTCGTCAGCTTTTTTG  >  NC_000913/4056915‑4057345

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