Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,569,555	G→A	P312L (CCT→CTT)	gadC ←	L‑glutamate:4‑aminobutyrate antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,569,555	0	G	A	97.0%	108.2 / NA	33	P312L (CCT→CTT)	gadC	L‑glutamate:4‑aminobutyrate antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/23);  minor base C (1/0);  total (10/23)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAGGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGC  >  NC_000913/1569308‑1569795                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 aaGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATAAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAACGAc                                                                                                                                                                                                                                               >  2:376533/1‑251 (MQ=255)    aaCCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACcaa                                                                                                                                                                                                                                            <  2:338622/251‑1 (MQ=255)                  aGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCCCCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGc                                                                                                                                                                                                                               <  2:210947/250‑1 (MQ=255)                    ccAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGAt                                                                                                                                                                                                                            <  1:135131/251‑1 (MQ=255)                    ccAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGGGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCGGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGAt                                                                                                                                                                                                                            <  2:225416/251‑1 (MQ=255)                                aTGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTACGCCAGAac                                                                                                                                                                                                                >  2:237869/1‑251 (MQ=255)                                   aaaTAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACAcc                                                                                                                                                                                                             <  2:119942/251‑1 (MQ=255)                                                              gTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCAcgcg                                                                                                                                                                                  >  1:314966/1‑251 (MQ=255)                                                                cAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGaa                                                                                                                                                                                <  1:34743/251‑1 (MQ=255)                                                                       aGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTCTTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTcc                                                                                                                                                                          >  2:47998/1‑250 (MQ=255)                                                                         tGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCAcc                                                                                                                                                                       >  2:283118/1‑250 (MQ=255)                                                                         ggCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGTTTTTTCTGCGCTGTTACATACATCCCTCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTTCGGAGATCACGCGAACCGTCCACt                                                                                                                                                                       <  1:281644/250‑1 (MQ=255)                                                                                  gaaggaCATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTg                                                                                                                                                              <  1:135034/251‑1 (MQ=255)                                                                                    aggaCATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGt                                                                                                                                                            >  2:169153/1‑251 (MQ=255)                                                                                     ggaCATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTg                                                                                                                                                           <  1:7082/251‑1 (MQ=255)                                                                                     gaaCATGTTGTGACCGCCACCGGGATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGGTATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCTGGTTTTTCTGCGCTGTTACATACAGCCCGCGAGAAAGACCAACAAGCCAGGAGGCGATTTCCGCCAGAACACCCAGCAGCAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTg                                                                                                                                                           <  2:260655/249‑1 (MQ=255)                                                                                                     ccACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGAcc                                                                                                                                                                                                                                              <  1:227796/152‑1 (MQ=255)                                                                                                     ccACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGAcc                                                                                                                                                                                                                                              >  2:227796/1‑152 (MQ=255)                                                                                                      cACCGGTAGTGGCGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGa                                                                                                                                          <  2:346825/251‑1 (MQ=255)                                                                                                         cGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGAC‑AACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACgg                                                                                                                                      >  1:58414/1‑251 (MQ=255)                                                                                                                         gatCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCAtt                                                                                                                        <  2:314966/250‑1 (MQ=255)                                                                                                                           tCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTAcc                                                                                                                     >  2:313397/1‑251 (MQ=255)                                                                                                                                     aGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACAACCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGgaga                                                                                                           <  2:237851/251‑1 (MQ=255)                                                                                                                                        cGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGt                                                                                                        <  1:313397/251‑1 (MQ=255)                                                                                                                                                     tGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTAcc                                                                                           >  1:232206/1‑251 (MQ=255)                                                                                                                                                        gAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCgg                                                                                        <  2:145919/251‑1 (MQ=255)                                                                                                                                                          aaTGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGAC‑AACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAAt                                                                                     <  2:58414/251‑1 (MQ=255)                                                                                                                                                                  gCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAATGACCATc                                                                              <  2:114996/251‑1 (MQ=255)                                                                                                                                                                                 cATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAATGACCATCGCAATAGACAAacca                                                               <  1:122082/251‑1 (MQ=255)                                                                                                                                                                                                  gCGAATGCCGCAGGCAGCAGGGTTTTCTGCGCGGTTACATACATCCCGCGAGAAAGACCAAAAATCCAGGAGGCGGTTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGGCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGGTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACGCGGAATGACCATCGCAATAGACAAGCCACCAACAGAGCTTAAGCa                                              <  2:183721/251‑1 (MQ=255)                                                                                                                                                                                                      aTGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGAGCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACAGCAGAACGGTAAAGGTGTGCATTACCCCTGCGGAGAGGTTGATTTCATTACGCGGAATGACCATCGCAATAGACAAACCACCAACAGCGCTTAAGCAGAtt                                          <  2:381402/251‑1 (MQ=255)                                                                                                                                                                                                                              gcgcTGTTACATACATCCCGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACa                    <  1:246184/249‑1 (MQ=255)                                                                                                                                                                                                                                              cccGCGAGAAAGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTcgc  <  1:169153/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 AAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATTTTTGTTCATTTTAGCGAATGCCGCTGGCAGCAGGTTTTTCTGCGCTGTTACATACATCCCGCGAGAAGGACCAACAATCCAGGAGGCGATTTCCGCCAGAACACCCAGCAACAGCAGTGCGGAGATCACGCGAACCGTCCACTCAATTTCTGGTGCCACATGGGACATCAGAACGGTAAAGGTTTGCATTACCCCTGCGGAGAGGTTGATTTCATTACCCGGAATGACCATCGCAATAGACAAACCACCAACAGAGCTTAAGCAGATTGCCGCCACCATCAGCAGTAACATAGCCAACGGATAGTCGC  >  NC_000913/1569308‑1569795

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