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

GGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGC  >  NC_000913/3850844‑3850934                                                |                                            ggaCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGc                                           <  1:2622076/50‑1 (MQ=255)   aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                          >  1:5624800/1‑49 (MQ=255)   aCTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCt                                          >  1:4013154/1‑49 (MQ=255)    cTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTa                                        >  1:3395270/1‑50 (MQ=255)    cTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTa                                        >  1:3930942/1‑50 (MQ=255)    cTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTa                                        >  1:1880837/1‑50 (MQ=255)     tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:4564621/1‑50 (MQ=255)     tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:5363295/1‑50 (MQ=255)     tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:3497519/1‑50 (MQ=255)     tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:2268803/1‑50 (MQ=255)     tGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:1956316/1‑50 (MQ=255)      gATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:4036356/1‑49 (MQ=255)      gATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTat                                       >  1:3382942/1‑49 (MQ=255)       aTTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatca                                     <  1:673373/50‑1 (MQ=255)       aTTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatca                                     <  1:740774/50‑1 (MQ=255)       aTTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatc                                      <  1:6022521/49‑1 (MQ=255)        ttACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcat                                    >  1:2135524/1‑50 (MQ=255)        ttACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcat                                    <  1:4937502/50‑1 (MQ=255)         tACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTatcatc                                   >  1:2974536/1‑50 (MQ=255)           cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCt                                  <  1:248338/49‑1 (MQ=255)           cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCt                                  <  1:1560355/49‑1 (MQ=255)           cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCt                                  <  1:5003749/49‑1 (MQ=255)           cGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTg                                 <  1:2521309/50‑1 (MQ=255)             cTGCACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCt                               >  1:1154417/1‑50 (MQ=255)                cACTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTc                              >  1:2228995/1‑48 (MQ=255)                 aCTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTc                              >  1:3836046/1‑47 (MQ=255)                  cTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                          <  1:1285922/50‑1 (MQ=255)                  cTTTCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGa                          >  1:3134989/1‑50 (MQ=255)                   tttCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGAc                         >  1:5723724/1‑50 (MQ=255)                     tCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGAc                         <  1:1546517/48‑1 (MQ=255)                     tCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGAc                         <  1:1725989/48‑1 (MQ=255)                     tCACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGc                       >  1:4613793/1‑50 (MQ=255)                      cACGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                      >  1:4653360/1‑50 (MQ=255)                       aCGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                      <  1:5812873/49‑1 (MQ=255)                       aCGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCt                      <  1:3635752/49‑1 (MQ=255)                       aCGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTg                     <  1:315662/50‑1 (MQ=255)                       aCGACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTg                     <  1:1625334/50‑1 (MQ=255)                         gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgt                   >  1:4368025/1‑50 (MQ=255)                         gACATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgt                   >  1:2121665/1‑50 (MQ=255)                          aCATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtc                  <  1:5878191/50‑1 (MQ=255)                           cATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcg                 >  1:3767606/1‑50 (MQ=255)                           cATCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcg                 >  1:3111241/1‑50 (MQ=255)                            aTCTTCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGgtcgt                >  1:409157/1‑50 (MQ=255)                               ttCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCAt             >  1:952130/1‑50 (MQ=255)                                tCCAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCAtt            >  1:4256132/1‑50 (MQ=255)                                  cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTg           <  1:2785927/49‑1 (MQ=255)                                  cAGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGa          <  1:1661374/50‑1 (MQ=255)                                   aGCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc         >  1:3653654/1‑50 (MQ=255)                                    gCTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc         >  1:2631662/1‑49 (MQ=255)                                     cTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc         >  1:1618467/1‑48 (MQ=255)                                     cTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGAc         >  1:3572334/1‑48 (MQ=255)                                     cTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCa       <  1:3955307/50‑1 (MQ=255)                                     cTTATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCa       <  1:2193908/50‑1 (MQ=255)                                      ttATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg      >  1:892305/1‑50 (MQ=255)                                      ttATCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg      >  1:619251/1‑50 (MQ=255)                                        aTCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAg      >  1:5419703/1‑48 (MQ=255)                                         tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTa    <  1:1484626/49‑1 (MQ=255)                                         tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTa    <  1:2814452/49‑1 (MQ=255)                                         tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAg   <  1:630433/50‑1 (MQ=255)                                         tCGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAg   <  1:456691/50‑1 (MQ=255)                                          cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  >  1:3254988/1‑50 (MQ=255)                                          cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  >  1:5127273/1‑50 (MQ=255)                                          cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  >  1:4143162/1‑50 (MQ=255)                                          cGATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  >  1:1543622/1‑50 (MQ=255)                                           gATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  <  1:5593454/49‑1 (MQ=255)                                           gATTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  >  1:5986528/1‑49 (MQ=255)                                            aTTTAGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGc  >  1:4539828/1‑48 (MQ=255)                                                |                                            GGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCTGGTCGTCATTGACCAGTAGC  >  NC_000913/3850844‑3850934

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