Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,734,042	A→C	18.7%	intergenic (‑68/+89)	proY ← / ← hutH	proline (histidine) APC transporter/histidine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,734,042	0	A	C	18.7%	90.0 / 2.6	45	intergenic (‑68/+89)	proY/hutH	proline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (14/22);  new base C (7/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 6.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACA  >  NC_002947/5733895‑5734189                                                                                                                                                     |                                                                                                                                                     gTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCAc                                                                                                                                                     <  1:252144/149‑1 (MQ=255)    ccGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGc                                                                                                                                                   <  4:230004/148‑1 (MQ=255)         ggCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCCCGCCCGCTc                                                                                                                                             >  3:99582/1‑149 (MQ=255)           cGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                            >  5:224798/1‑148 (MQ=255)                 cGAGGTCAATTAAACGTATGTGTCGGTCACTTAGCCCGCTCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTTTTTTTGTGGTCTTGAGATTCAGTGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAtgg                                                                                                                                     <  8:242/149‑1 (MQ=255)                     ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATg                                                                                                                                 <  6:224798/149‑1 (MQ=255)                        cATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCCCGCCCGCTCCCACAGGGGATGgtg                                                                                                                              >  4:228858/1‑149 (MQ=255)                            aaaCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtttt                                                                                                                          >  8:95918/1‑149 (MQ=255)                            aaaCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtttt                                                                                                                          <  3:207611/149‑1 (MQ=255)                            aaaCGGATGTGGAGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtttt                                                                                                                          <  6:9245/149‑1 (MQ=255)                                ggATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCAc                                                                                                                      <  7:95918/149‑1 (MQ=255)                                  aTGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCt                                                                                                                    <  5:37426/149‑1 (MQ=255)                                     tgGCCGGCACTTATCCCGCGCGGGCGCCCTTCAGCTTGTTGCATGTCACCTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGTGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTAc                                                                                                                  <  4:194613/148‑1 (MQ=255)                                       gCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTc                                                                                                               >  4:29769/1‑149 (MQ=255)                                          ggCACTTAGCCCGCGGTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGTGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGTGATGGTGTTTTGCACCTACTTCtgt                                                                                                            <  4:254366/149‑1 (MQ=255)                                                    ccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                                                  >  6:286159/1‑149 (MQ=255)                                                    ccGCGCTGGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                                                  <  1:209649/149‑1 (MQ=255)                                                     cgcgCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTg                                                                                                                         <  3:224995/125‑1 (MQ=255)                                                            gagaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCCCCCCCCCACCCACACGGGCTGGGGTTTTTACCCCAAATGTGTGGGCGCCGGCGCCCCcgc                                                                                          >  1:135914/1‑149 (MQ=255)                                                                 cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCCCCCCCCCGCCCCCACGGGATGGGGTTTTTCACCGACTTGAGGTGGAGCCGGCGCGCCCGCGaacg                                                                                     >  4:72528/1‑147 (MQ=255)                                                                  ttGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTATGTGGGAGCGGGCGCGCCCGCGAAggg                                                                                    >  7:167313/1‑149 (MQ=255)                                                                       tttgttgCTTTTCACGTCTTTTATTTTTTTTTGTGGGCGTTATATTCTGGGGGCTGCTGCGCA‑GCCCCTTCGGGGGCACGCCGGCTCCCACTGGGTTTGGTGTTTTGCACCTACTTGTGTGGGTGGGGGCGCGCCCGCGAAGGgccgc                                                                                <  4:99582/147‑1 (MQ=255)                                                                                 gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGccc                                                                     >  1:234438/1‑149 (MQ=255)                                                                                 gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACCCCCGCTCCCACAGGGGATGGTGTTTTGCACCGACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑CGCGGcccc                                                                    >  2:28922/1‑149 (MQ=255)                                                                                  tcccGTCCTTGAATTGTTTTTGTGGCCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGgccgc                                                                                <  3:219828/134‑1 (MQ=255)                                                                                  tCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCCCGACCACAGGGGATGGGGTTTTTGACCGACGTGTGTGGGGGCGGGCGCGCCCGCGCAGGGCCGCCCGAGCggga                                                                      >  4:219828/1‑145 (MQ=255)                                                                                   cACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCCCCCCCCCAACCACCGCGGGAGGGGTTTTGAACCAACTAATGTGGGACCGGGGCCGCCCCCGCACGGCACCCC‑CCCCGCCCCac                                                                  >  4:169907/1‑148 (MQ=255)                                                                                       tcgtttaaTGGTTTTTGGGGGCGTGACATCCAGGGGCCTGCCGCGCA‑GCCCCGTCGCGGGCCCGCCCGCTCCCACCGGGGAGGGTGTTTTGCACCTACTGCTGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGc                                                              <  4:267335/143‑1 (MQ=255)                                                                                       tCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGc                                                              <  6:226442/149‑1 (MQ=255)                                                                                         cTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGt                                                            >  4:198155/1‑149 (MQ=255)                                                                                             aaTTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACCCCCGCTCCCACAGGGGATGGTGGTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCCCGCAGGGCCGCCC‑AGCGGCCCCCGGGGGGTTa                                                          >  4:195554/1‑147 (MQ=255)                                                                                             aaTTGTTTTTGTGGTCGTGAGATGCCGGGGGCTGCTGCGCA‑GCCCCTGCGCGGGCACGCCCGCTCCCACAGGGGATGGGGTTTGGCACCTACTTGTGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACa                                                        <  3:72528/149‑1 (MQ=255)                                                                                                  tttttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCCCCCCCCCAACCACAGGGGATGGGGTTTTGGACCTACGTATGTGGGAGCGGGCCCGCCCCCCGAAGGCCCCCC‑CGCCGCCCCCCGGGCGTTAAAGGCtt                                                   >  4:110512/1‑149 (MQ=255)                                                                                                        tgGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGGGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGt                                             <  3:228858/149‑1 (MQ=255)                                                                                                           tagTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTaca                                          >  7:287852/3‑149 (MQ=255)                                                                                                           tCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCGCTTCGCGGGCCCGCCCGCTCCCACAGGGGATGGTGTTTTGGACCGACTTCTGTGGGAGCGGGCGCGCCCGCGCAGGGCCGCCC‑CGCGGCCCCCGGGGCGGTTCAGGCTTGGCAGAaca                                          >  1:26940/1‑149 (MQ=255)                                                                                                            cGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTacac                                         <  3:29769/149‑1 (MQ=255)                                                                                                              tGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATTGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                               >  5:73925/1‑141 (MQ=255)                                                                                                              tGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                               <  6:73925/141‑1 (MQ=255)                                                                                                                           gctgctGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCAACAGGGGATGGTGTTTGGCACCTACTTCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGc                          <  2:234438/149‑1 (MQ=255)                                                                                                                             tgctgcGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCCGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCCCCCGCACCAAGccc                        >  2:13135/1‑148 (MQ=255)                                                                                                                               gtgcGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGt                      <  8:167313/148‑1 (MQ=255)                                                                                                                                         cgcTTCGCGGGCCGGCCCGCTCCCCCAGGGTTGGGTGTTTTGCCCCTACTGCTGTGGGTGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTAcc             <  3:169907/147‑1 (MQ=255)                                                                                                                                                    cACGCCCGCTCCCACAGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACa  <  3:110512/149‑1 (MQ=255)                                                                                                                                                    cACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCAGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCaa    >  7:271796/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCA‑GCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCC‑AGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACA  >  NC_002947/5733895‑5734189

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