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

ACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCA  >  NC_000913/3850846‑3850936                                              |                                              aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                            <  1:1225079/49‑1 (MQ=255)  cTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCtt                                           <  1:2092253/49‑1 (MQ=255)   tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTa                                          >  1:1161905/1‑49 (MQ=255)    gATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                         >  1:744143/1‑49 (MQ=255)    gATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                         >  1:475090/1‑49 (MQ=255)      ttACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatca                                       >  1:818138/1‑49 (MQ=255)       tACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcat                                      >  1:2219766/1‑49 (MQ=255)         cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcatc                                     <  1:158415/48‑1 (MQ=255)            tGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCt                                 >  1:661219/1‑49 (MQ=255)                   tCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                            <  1:28720/47‑1 (MQ=255)                   tCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                            <  1:2618574/47‑1 (MQ=255)                   tCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGAc                           <  1:521547/48‑1 (MQ=255)                    cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGc                         >  1:929907/1‑49 (MQ=255)                    cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGc                         >  1:307463/1‑49 (MQ=255)                     aCGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGc                         <  1:2113711/48‑1 (MQ=255)                      cGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTg                       >  1:1579429/1‑49 (MQ=255)                       gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTgg                      >  1:54874/1‑49 (MQ=255)                       gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTgg                      >  1:2238918/1‑49 (MQ=255)                       gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTgg                      >  1:413790/1‑49 (MQ=255)                       gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTgg                      >  1:2866397/1‑49 (MQ=255)                        aCATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgt                     <  1:1240643/49‑1 (MQ=255)                          aTCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcg                   >  1:1377456/1‑49 (MQ=255)                            cTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcgt                  >  1:2501280/1‑48 (MQ=255)                               ccAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcgt                  <  1:1228221/45‑1 (MQ=255)                               ccAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcgt                  <  1:796301/45‑1 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCAtt              >  1:2860148/1‑48 (MQ=255)                                cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCAtt              <  1:2226220/48‑1 (MQ=255)                                      aTCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg        >  1:357206/1‑48 (MQ=255)                                       tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTa      <  1:1062532/49‑1 (MQ=255)                                        cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAg     >  1:1580986/1‑49 (MQ=255)                                         gATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAg     <  1:905424/48‑1 (MQ=255)                                          aTTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc    >  1:671825/1‑48 (MQ=255)                                          aTTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc    >  1:114155/1‑48 (MQ=255)                                           tttAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc    >  1:2062886/1‑47 (MQ=255)                                           tttAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc    >  1:1972619/1‑47 (MQ=255)                                           tttAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCa  <  1:1089634/49‑1 (MQ=255)                                              |                                              ACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGCCA  >  NC_000913/3850846‑3850936

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