Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,085,270	T→A	10.6%	M127K (ATG→AAG) ‡	aspA →	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,085,270	0	T	A	10.6%	117.0 / 6.0	47	M127K (ATG→AAG) ‡	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base A (2/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

GCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAGCCTGA  >  NC_002947/6085124‑6085416                                                                                                                                                    |                                                                                                                                                   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                                                                                                                                         <  5:89197/109‑2 (MQ=255)                                                  cACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGTGGCCCATGGCCTCCAgc                                                                                                                                         >  6:89197/1‑108 (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                                                                                     <  5:246512/149‑1 (MQ=255)                                                              gTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTc                                                                                     >  6:246512/1‑149 (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)                                                                                                                                                  cATGGGCCACCAGAAGGGTGAGTACCAGTCCCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGAAGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAGCCTGa  >  7:80002/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCACGCTGCCATTTCCGCAGCCTGCGCCCGCCTGATCAAAGGCGATTTCCACGACCAGTTCGTGGTAGA‑CATGATCCAGGGCGGCGCTGGTACTTCTACCAACATGAACGCCAACGAAGTCATTGCCAACGTTGCGCTGGAGGCCATGGGCCACCAGAAGGGTGAGTACCAGTACCTGCACCCTAACAACGACGTGAACATGGCGCAGTCGACCAACGACGCCTACCCGACCGCAATCCGCCTGGGCCTGCTGCTGGGCCACGACGCGCTGCTGGCCAGCCTTGATAGCCTGA  >  NC_002947/6085124‑6085416

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