Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,085,264	A→T	10.4%	E125V (GAG→GTG)	aspA →	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,085,264	0	A	T	10.4%	126.6 / 5.8	48	E125V (GAG→GTG)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base T (2/3);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

ATGCCAAGCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAGC  >  NC_002947/6085117‑6085412                                                                                                                                                     |                                                                                                                                                     aTGCCAAGCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAg                                                                                                                                                     <  8:160135/149‑1 (MQ=255)      aaGCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGc                                                                                                                                                   >  1:193015/1‑146 (MQ=255)      aaGCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGc                                                                                                                                                   <  2:193015/146‑1 (MQ=255)        gCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATgg                                                                                                                                              <  1:231906/149‑1 (MQ=255)        gCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCAGTGCCAACGTTGCGCTGGAGGCCATgg                                                                                                                                              >  6:120994/1‑149 (MQ=255)          aCGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGc                                                                                                                                            >  3:212442/1‑149 (MQ=255)            gCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGcca                                                                                                                                          >  4:103722/1‑149 (MQ=255)             cTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGccac                                                                                                                                         <  3:70216/149‑1 (MQ=255)              tGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATgg                                                                                                                                              >  8:161846/1‑143 (MQ=255)              tGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATgg                                                                                                                                              <  7:161846/143‑1 (MQ=255)              tGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGccac                                                                                                                                         >  2:97751/1‑148 (MQ=255)                ccATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGTGGCCCATGGCCTCCAg                                                                                                                                      >  2:292770/1‑149 (MQ=255)                 cATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGTGGCCCATGGCCTCCAgc                                                                                                                                     <  7:246891/149‑2 (MQ=255)                 cATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGTGGCCCATGGCCTCCAgc                                                                                                                                     <  1:292770/149‑2 (MQ=255)                  aTTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAg                                                                                                                                      >  5:189845/1‑147 (MQ=255)                   tttCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGATCGTACTTCTACAAACATTAACGCCAAAGAACTCATTGCCAACGTTGCGCTGGAGGCCATTGGccaccac                                                                                                                                      >  7:149871/1‑145 (MQ=255)                       cGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGt                                                                                                                                <  1:97751/148‑1 (MQ=255)                          aGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGagt                                                                                                                            <  5:319570/149‑1 (MQ=255)                                           tCAAAGGCGATTTCCACGACCAGTTCGTGGTAGAGCATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGagtacc                                                                                                                         >  1:205830/1‑136 (MQ=255)                                           tCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGagtacc                                                                                                                         <  2:205830/135‑1 (MQ=255)                                           tCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCAccc                                                                                                            >  1:24717/1‑148 (MQ=255)                                            cAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCt                                                                                                           >  7:233053/1‑148 (MQ=255)                                                 gCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTaacaac                                                                                                     >  3:128328/1‑149 (MQ=255)                                                     ctaccACGAACTGGTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCGCCAGAAGGGTGAGTACCAGTACCTGCAc                                                                                                              >  5:277291/4‑136 (MQ=255)                                                     ctaccACGAACTGGTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCAc                                                                                                              <  6:277291/133‑1 (MQ=255)                                                         cACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGTGGCCCATGGCCTCCAgc                                                                                                                                     >  6:89197/1‑108 (MQ=255)                                                         cACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGTGGCCCATGGCCTCCAgc                                                                                                                                     <  5:89197/109‑2 (MQ=255)                                                         cACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGa                                                                                              <  4:31152/148‑1 (MQ=255)                                                             aCCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTAATTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATg                                                                                         <  8:149871/149‑1 (MQ=255)                                                                aGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTTCGCTGGAGCCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGgcg                                                                                      <  1:302634/149‑1 (MQ=255)                                                                     gTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTc                                                                                 >  6:246512/1‑149 (MQ=255)                                                                     gTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTc                                                                                 <  5:246512/149‑1 (MQ=255)                                                                       ggTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGCCCATGGCCCACCAGAAGGGTGACTACCAGTACCTGCACCCAAACAACGACGAGAACATGGCGCAGTCGa                                                                               >  5:68273/1‑149 (MQ=255)                                                                        gTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGa                                                                               <  3:103722/148‑1 (MQ=255)                                                                           gA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCaa                                                                           <  6:189845/149‑1 (MQ=255)                                                                                 gATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAacgacg                                                                      >  4:33207/1‑149 (MQ=255)                                                                                          cggcgCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGTGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGa                                                             <  4:212442/149‑1 (MQ=255)                                                                                                  gTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCaa                                                       <  1:186229/147‑1 (MQ=255)                                                                                                    aCTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACACGACCGCAATACg                                                   >  7:241289/1‑149 (MQ=255)                                                                                                         tACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTgg                                              <  3:33207/149‑1 (MQ=255)                                                                                                           ccAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGc                                            <  2:24717/149‑1 (MQ=255)                                                                                                              aCATGAACGCCAAAGAAGTCATTGCCAACGTTGCGCTGGAGGGCATGGGCCATCAGAAGGGTGAGTATCATTACCTGCACCCTATCAACGAGTTTTACATGGGGCAGTCGACCATCGACGCCTACCCGACCGCAATCCGCCTGGGCctg                                         <  6:68273/149‑1 (MQ=255)                                                                                                                        cAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGcc                               >  3:128953/1‑149 (MQ=255)                                                                                                                                  aTTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACgcgc                      <  3:96786/148‑1 (MQ=255)                                                                                                                                       cAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTgg                >  6:292203/1‑149 (MQ=255)                                                                                                                                            ttGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGc           >  3:233022/1‑149 (MQ=255)                                                                                                                                               cgcTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCt         <  5:292203/148‑1 (MQ=255)                                                                                                                                                    gAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAg   >  4:292435/1‑149 (MQ=255)                                                                                                                                                     aGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAGc  >  4:115101/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATGCCAAGCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAGC  >  NC_002947/6085117‑6085412

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