Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,200,368	G→A	Q200* (CAG→TAG)	ygiF ←	inorganic triphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,200,368	0	G	A	97.1%	115.2 / NA	34	Q200* (CAG→TAG)	ygiF	inorganic triphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/16);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCCCTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTGCGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTGGTTCGC  >  NC_000913/3200290‑3200443                                                                               |                                                                            ccAGCCCCTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTacg                                                                           >  2:595034/1‑81 (MQ=255)  cAGCCCCTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTacg                                                                           <  1:594966/80‑1 (MQ=255)     ccccTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAg                                                                       <  1:594967/81‑1 (MQ=255)       ccTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAgat                                                                     <  2:594968/81‑1 (MQ=255)        cTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAgat                                                                     >  1:595020/1‑80 (MQ=255)         tGTTCCACATCCTCTTTTTCCGCCACATGCAACATGGTGGTCGGTTTGATTTCACTCGCCGGATTTCCCTACCCCAgttga                                                                   >  2:1035738/1‑81 (MQ=255)               acaTCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCa                                                             <  1:594969/81‑1 (MQ=255)                       ttttGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGcgccgc                                                     >  2:595047/1‑81 (MQ=255)                             cGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCt                                               <  1:594970/81‑1 (MQ=255)                               cAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCt                                               >  2:595026/1‑79 (MQ=255)                                aaCATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTGAg                                            >  2:595061/1‑81 (MQ=255)                                aaCATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAg                                            >  1:595050/1‑81 (MQ=255)                                aaCATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAg                                            <  1:594971/81‑1 (MQ=255)                                 aCATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAgg                                           >  1:595055/1‑81 (MQ=255)                                   aTGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCt                                         <  1:594972/81‑1 (MQ=255)                                      cAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGcc                                      <  2:594973/81‑1 (MQ=255)                                       aaaaTGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGccc                                     >  2:595054/1‑81 (MQ=255)                                       aaaaTGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGccc                                     >  2:595039/1‑81 (MQ=255)                                           tggtggTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGc                                 <  2:594975/81‑1 (MQ=255)                                           tggtggTCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGc                                 <  1:594974/81‑1 (MQ=255)                                            ggtggtCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCg                                >  2:595042/1‑81 (MQ=255)                                             gtggtCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGt                               <  2:594977/81‑1 (MQ=255)                                             gtggtCGGTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGt                               <  1:594976/81‑1 (MQ=255)                                                    gTTTGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCgg                        <  2:594978/81‑1 (MQ=255)                                                       tGATTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGtt                      <  2:594979/80‑1 (MQ=255)                                                         aTTTCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGc                   >  2:595040/1‑81 (MQ=255)                                                          tttCACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCg                  <  1:594980/81‑1 (MQ=255)                                                             cACGCGCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATa               >  2:595051/1‑81 (MQ=255)                                                               cgcgCCGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATAc              >  1:595046/1‑80 (MQ=255)                                                                    cGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTg        >  2:595053/1‑81 (MQ=255)                                                                    cGGATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTg        <  1:594981/81‑1 (MQ=255)                                                                     ggATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTgg       >  1:595057/1‑81 (MQ=255)                                                                      gATTGCCCTGCGCCAGATGATAACCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTGGt      <  1:594982/81‑1 (MQ=255)                                                                      gATTGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTGGt      >  1:595060/1‑81 (MQ=255)                                                                        ttGCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTGGTTc    <  2:594983/81‑1 (MQ=255)                                                                          gCCCTACGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTGGTTCGc  >  2:595067/1‑81 (MQ=255)                                                                               |                                                                            CCAGCCCCTGTTCCACATCGGCTTTTGCCGCAACATGCAAAATGGTGGTCGGTTTGATTTCACGCGCCGGATTGCCCTGCGCCAGATGATAGCCACGCGCCGCTTTGCTCAGGCTGCCCTGGCGTAATCCGGTTTGCGATACCAGTTGGTTCGC  >  NC_000913/3200290‑3200443

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