Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,850,891	G→A	Q68* (CAA→TAA)	ilvN ←	acetohydroxy acid synthase I subunit IlvN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,850,891	0	G	A	100.0%	101.7 / NA	29	Q68* (CAA→TAA)	ilvN	acetohydroxy acid synthase I subunit IlvN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/7);  total (22/7)

GGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCAG  >  NC_000913/3850844‑3850937                                                |                                               ggaCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGc                                              >  1:305491/1‑50 (MQ=255)  gaCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                             >  1:2031407/1‑50 (MQ=255)   aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCtt                                            <  1:1834336/50‑1 (MQ=255)   aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                             >  1:2323880/1‑49 (MQ=255)      gATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                          >  1:178200/1‑49 (MQ=255)         tACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcatc                                      >  1:1641104/1‑50 (MQ=255)            gCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGc                                   >  1:40256/1‑50 (MQ=255)             cTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCt                                  <  1:1808348/50‑1 (MQ=255)              tGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTc                                 >  1:1229942/1‑50 (MQ=255)                cACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTc                                 >  1:467768/1‑48 (MQ=255)                cACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCa                               <  1:794402/50‑1 (MQ=255)                 aCTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTc                                 >  1:233473/1‑47 (MQ=255)                 aCTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAg                              <  1:1848854/50‑1 (MQ=255)                  cTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                             <  1:2121071/50‑1 (MQ=255)                  cTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                             >  1:431077/1‑50 (MQ=255)                  cTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                             >  1:1043035/1‑50 (MQ=255)                   tttCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGAc                            >  1:738213/1‑50 (MQ=255)                      cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                         >  1:976919/1‑50 (MQ=255)                       aCGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTg                        >  1:2148584/1‑50 (MQ=255)                        cGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTgg                       >  1:266028/1‑50 (MQ=255)                         gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgt                      >  1:2199155/1‑50 (MQ=255)                          aCATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtc                     >  1:584838/1‑50 (MQ=255)                           cATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcg                    >  1:121162/1‑50 (MQ=255)                           cATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcg                    <  1:1615719/50‑1 (MQ=255)                                    gCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc            >  1:2212144/1‑49 (MQ=255)                                     cTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc            >  1:224286/1‑48 (MQ=255)                                        aTCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg         >  1:1847931/1‑48 (MQ=255)                                          cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc     >  1:1538727/1‑50 (MQ=255)                                             tttAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCAg  <  1:216049/50‑1 (MQ=255)                                                |                                               GGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCAG  >  NC_000913/3850844‑3850937

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