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%	145.3 / NA	37	Q68* (CAA→TAA)	ilvN	acetohydroxy acid synthase I subunit IlvN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/19);  total (18/19)

ACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCA  >  NC_000913/3850846‑3850936                                              |                                              aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                            >  1:517073/1‑49 (MQ=255) aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                            >  1:2767557/1‑49 (MQ=255)  cTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTa                                          <  1:1591214/50‑1 (MQ=255)   tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                         >  1:1475346/1‑50 (MQ=255)     aTTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatc                                        <  1:1452668/49‑1 (MQ=255)     aTTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatc                                        <  1:1943607/49‑1 (MQ=255)      ttACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcat                                      <  1:468516/50‑1 (MQ=255)         cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCt                                    <  1:1108726/49‑1 (MQ=255)         cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCt                                    <  1:125409/49‑1 (MQ=255)          gCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGc                                  >  1:2263154/1‑50 (MQ=255)              cACTTTCACGACATCTCCCAGCTTATCGATTTAGCTTATCATCTGCTc                                >  1:2426770/1‑48 (MQ=255)                cTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                            <  1:2220105/50‑1 (MQ=255)                    cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGc                         <  1:425417/49‑1 (MQ=255)                    cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                        >  1:940889/1‑50 (MQ=255)                    cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                        >  1:681564/1‑50 (MQ=255)                    cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                        >  1:509891/1‑50 (MQ=255)                        aCATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtc                    <  1:1366941/50‑1 (MQ=255)                        aCATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtc                    >  1:3261/1‑50 (MQ=255)                            cTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCa                >  1:1676431/1‑50 (MQ=255)                            cTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCa                >  1:2327600/1‑50 (MQ=255)                              tCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCAtt              >  1:1217525/1‑50 (MQ=255)                               ccAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCAt               <  1:1488112/48‑1 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTg             <  1:2076663/49‑1 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTg             <  1:600988/49‑1 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGa            <  1:2072534/50‑1 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGa            >  1:232553/1‑50 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGa            <  1:2249458/50‑1 (MQ=255)                                  gCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAcc          <  1:938090/50‑1 (MQ=255)                                   cTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc           >  1:2467635/1‑48 (MQ=255)                                      aTCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg        >  1:1522181/1‑48 (MQ=255)                                      aTCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg        >  1:583294/1‑48 (MQ=255)                                      aTCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGt       <  1:2026843/49‑1 (MQ=255)                                       tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTa      <  1:2138149/49‑1 (MQ=255)                                       tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAg     <  1:443231/50‑1 (MQ=255)                                        cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc    >  1:2760763/1‑50 (MQ=255)                                          aTTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCa  <  1:2556399/50‑1 (MQ=255)                                           tttAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc    >  1:2023830/1‑47 (MQ=255)                                              |                                              ACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCA  >  NC_000913/3850846‑3850936

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: ‑