Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,568,848	C→T	A408T (GCC→ACC)	gadB ←	glutamate decarboxylase isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,568,848	0	C	T	100.0%	154.9 / NA	65	A408T (GCC→ACC)	gadB	glutamate decarboxylase isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/40);  total (25/40)

ATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAA  >  NC_000913/1568799‑1568936                                                  |                                                                                         aTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCa                                               >  1:4729503/1‑93 (MQ=14) aTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCa                                               <  1:3574765/93‑1 (MQ=14) aTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGAAGACGCAGACGTTCAGAGAGGTCa                                               >  2:1691455/1‑93 (MQ=9)  tAATGCGCATCACCACAATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                              <  1:4216543/93‑1 (MQ=12)  aaaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                              <  2:3071503/92‑1 (MQ=17)   aaTGCGCTGCGCCAGGATATCGGTGACTTCACCGCCGAGAGTGAGGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                             <  2:2478408/93‑1 (MQ=11)   aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                             <  2:2478407/93‑1 (MQ=18)   aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                             <  1:394658/93‑1 (MQ=18)    aTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATAc                                            >  1:4336029/1‑93 (MQ=18)     tGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACa                                           <  2:42243/93‑1 (MQ=18)       cgcATCCCCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAgg                                         <  1:1353497/93‑1 (MQ=255)       cgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAgg                                         <  1:1353499/93‑1 (MQ=21)       cgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAgg                                         <  2:3219050/93‑1 (MQ=21)        gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                        <  1:364783/93‑1 (MQ=21)        gcATCACCACGAGGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                        <  2:4584583/93‑1 (MQ=255)         cATCACCACGATCTCGGTGGCTTCACCACCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTAAGAGAGGTCATACAGGgt                                       <  1:2735651/93‑1 (MQ=255)          aTCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtg                                      >  2:915361/1‑93 (MQ=18)           tctccacGATGTGGGTGGGTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                     <  1:379979/90‑1 (MQ=11)           tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCGGACGTTCAGAGAGGTCATACAGGgggt                                     >  2:70962/1‑93 (MQ=255)           tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                     >  2:4494388/1‑93 (MQ=18)           tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                     >  2:1300494/1‑93 (MQ=18)            caccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTa                                    >  2:1576332/1‑93 (MQ=18)            caccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTa                                    >  1:4287240/1‑93 (MQ=18)            caccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTa                                    <  1:2256625/93‑1 (MQ=21)             accacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTAt                                   >  2:3089621/1‑93 (MQ=18)             accacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTAt                                   >  2:3619726/1‑93 (MQ=18)              ccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATc                                  <  2:2544050/93‑1 (MQ=21)               cacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATcc                                 <  2:840755/93‑1 (MQ=21)                acGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                <  2:748111/93‑1 (MQ=21)                acGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                <  1:1613979/93‑1 (MQ=18)                 cGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCg                               <  1:3129331/93‑1 (MQ=21)                  gATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCgg                              >  1:1330051/1‑93 (MQ=18)                  gATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCgg                              >  1:2594707/1‑93 (MQ=18)                   aTGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGa                             >  2:1964150/1‑93 (MQ=18)                   aTGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGa                             >  1:4372907/1‑93 (MQ=18)                    tGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGAt                            <  2:561328/93‑1 (MQ=21)                    tGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGAt                            >  2:3145709/1‑93 (MQ=18)                     gTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATc                           <  1:2306989/93‑1 (MQ=21)                     gTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATc                           <  1:931206/93‑1 (MQ=18)                      tCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCt                          <  2:400802/93‑1 (MQ=21)                       cGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCtt                         >  2:3710641/1‑93 (MQ=18)                        ggTGGCTGCGCCGCCGAGAGGGGAAGTCGAAACCTGCCAGCCGCGCAGACACAGCCGTTCCGAGGGGTCATACGCGGTGTATCCCTGATCTTc                        <  2:3173991/93‑1 (MQ=2)                         gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                       >  1:2815722/1‑93 (MQ=18)                            gCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCa                    <  2:4131689/93‑1 (MQ=21)                               tCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCt                 >  2:781303/1‑93 (MQ=18)                               tCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCt                 >  2:3593549/1‑93 (MQ=18)                               tCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCt                 <  2:49176/93‑1 (MQ=21)                                ccccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCtt                <  2:4506287/91‑1 (MQ=16)                                  ccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTc              <  1:2247700/93‑1 (MQ=21)                                  ccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCACGGATCTTCACCATCTTTc              >  2:3437787/1‑93 (MQ=16)                                   cgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCa             <  2:2296154/93‑1 (MQ=21)                                     ccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGt           <  1:3808422/93‑1 (MQ=21)                                      cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATATTCACCATCTTTCAGtt          >  2:4652127/1‑93 (MQ=16)                                       gagaGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGttt         <  1:4499675/93‑1 (MQ=21)                                       gagaGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGttt         >  2:4543342/1‑93 (MQ=18)                                       gagaGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGttt         >  1:3310559/1‑93 (MQ=18)                                       gagaGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGttt         <  2:2610735/93‑1 (MQ=21)                                       gagaGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGAATCTTCACCATCTTTCAGttt         <  1:212386/93‑1 (MQ=255)                                         gagTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGa       <  2:4670560/93‑1 (MQ=18)                                          agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa      <  1:3486298/93‑1 (MQ=21)                                          agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa      <  2:624167/93‑1 (MQ=21)                                          agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa      <  2:791488/93‑1 (MQ=21)                                          agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGCTTGaa      <  1:481080/93‑1 (MQ=16)                                           gTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAg     <  1:740609/93‑1 (MQ=18)                                           gTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAg     >  2:2819474/1‑93 (MQ=18)                                            tGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGc    <  1:1086236/93‑1 (MQ=21)                                             gAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa   <  2:4092799/93‑1 (MQ=18)                                             gAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa   <  2:2544084/93‑1 (MQ=18)                                              aaGGTCGGAACCTGCCAGCCGCGCAGACGCAGACTTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaa  <  2:1722394/93‑1 (MQ=16)                                                  |                                                                                         ATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAA  >  NC_000913/1568799‑1568936

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