Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,336,411	2 bp→AT	coding (83‑84/1260 nt)	murA ←	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,336,411	0	G	A	89.3%	73.9 / 2.4	28	I28N (ATC→AAT)	murA	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (13/12);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-02
*	NC_000913	3,336,412	0	A	T	89.3%	82.5 / 3.3	28	I28N (ATC→AAT)	murA	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base T (13/12);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

TTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATA  >  NC_000913/3336171‑3336638                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                   gtGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCGTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCagag                                                                                                                                                                                                                           <  1:54241/250‑1 (MQ=255)      ccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAgcagca                                                                                                                                                                                                                      <  1:162690/251‑1 (MQ=255)                 cgccttacccacggtcccATCGCCCATATATATGCACGCATGGTGTTACCCCGATTGTAATGTGCGCATAAGTCATTAACGTCGCTGGCATCAATATGCACAGAACCATTACTTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATTTACTTCTTTCAGTGTCGGGACGTTCTGTATCTCTACCGGTTCTTCCGCCCGTAGTGCGGCAAAAAGATTAGGCAGATCAGCATTTTTAGCGcc                                                                                                                                                                                                           <  2:183968/236‑1 (MQ=255)                                 ccAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTcc                                                                                                                                                                                           >  2:118664/1‑250 (MQ=255)                                 ccAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCg                                                                                                                                                                                           >  1:226123/1‑251 (MQ=255)                                   aGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGcc                                                                                                                                                                                         <  1:172880/251‑1 (MQ=255)                                          aGATAGAAGCTCGCTTGGTTTTAACCCGATCGGAAGGTGCCCAGAATACATTAACGGCTCGGGTATCAATATGCACAGAACCATGACGGTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATCGATGTATCGACGTCTTTCATGTTCGGGACGATCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGTCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCTCCCTGGAGc                                                                                                                                                                                  <  2:68300/251‑1 (MQ=255)                                            aaaGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGATTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGAGCCGGAAATTGTGAATTCGCCCTGGAGCtt                                                                                                                                                                                >  1:139373/3‑251 (MQ=255)                                                aaGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGtt                                                                                                                                                                            >  1:259660/1‑251 (MQ=255)                                                aaGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGtt                                                                                                                                                                            >  1:293376/1‑251 (MQ=255)                                                 aGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCATGGAGCTTCGTTg                                                                                                                                                                           >  1:49809/1‑251 (MQ=255)                                                           gTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCt                                                                                                                                                                                       <  1:71208/229‑1 (MQ=255)                                                           gTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGAAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCtagatcgacagagcgccgtgtgg                                                                                                                                                                 >  2:71208/1‑229 (MQ=255)                                                               tAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACg                                                                                                                                                              <  1:306299/250‑1 (MQ=255)                                                                    agaacGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAAttt                                                                                                                                                        <  2:176976/247‑1 (MQ=255)                                                                           aaGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCGGGCTTAGCAGCTTCATTGATGTATCGACGTCTGTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAtt                                                                                                                                                 <  2:293376/251‑1 (MQ=255)                                                                                 cgcAGAATACATTACCGTCGCGGGCATAATTCTTCACAGAACTCGTACGTTCTACTTGCGTGCCCAGCTGGCTTAGCAGCTTCATTGATGTAGCGACGTCTTTCAGTTTCGGGACGTTCGGTATTTATTCCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCAGTTTTAGCGCCGGAAATTGTGACTTCGCCCGGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGttt                                                                                                                                           <  2:58005/251‑1 (MQ=255)                                                                                          cATTAACTTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGt                                                                                                                                  <  1:279207/251‑1 (MQ=255)                                                                                             tAACGGCGCGGGCATCAATATGCACAGAACCATTACGGTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATTTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTATGCAGAGCAGCATTTTTATCGCCGGAAATTGTGACTTCGCCCTGGAGCGTTGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTa                                                                                                                                <  1:18453/250‑1 (MQ=255)                                                                                                                         aCCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGATCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCTCCGGAAGTTGTGACTTCGCCCTGGCGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGaaacg                                                                                                    >  2:188772/1‑247 (MQ=255)                                                                                                                                             cACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTc                                                                               <  1:92210/251‑1 (MQ=255)                                                                                                                                                      ggCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGc                                                                        <  1:188772/249‑1 (MQ=255)                                                                                                                                                             cagcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATcgcgcg                                                               >  1:329788/1‑251 (MQ=255)                                                                                                                                                                 ttCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGcccc                                                           >  1:201200/1‑250 (MQ=255)                                                                                                                                                                 ttCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGccc                                                            >  1:121651/1‑250 (MQ=255)                                                                                                                                                                       gATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTAATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTcagc                                                     >  1:288756/1‑248 (MQ=255)                                                                                                                                                                                               gggACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGacac                             >  1:185207/1‑251 (MQ=255)                                                                                                                                                                                                                           ccGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAatata  <  1:37870/250‑1 (MQ=255)                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                   TTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATA  >  NC_000913/3336171‑3336638

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑