Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,677,258	A→T	7.4%	L439Q (CTG→CAG)	ldcC ←	lysine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,677,258	0	A	T	7.4%	147.8 / 2.8	54	L439Q (CTG→CAG)	ldcC	lysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base T (2/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CTGCCAGATGCTGAACCACCAGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGG  >  NC_002947/4677113‑4677405                                                                                                                                                  |                                                                                                                                                     cTGCCAGATGCTGAACCACCAGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAgcg                                                                                                                                                   <  4:109561/149‑1 (MQ=255)  tGCCAGATGCTGAACCACCAGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCCgc                                                                                                                                                    >  6:279835/1‑147 (MQ=255)       gATGCTGAACCACCACTCCGCTGCCGCAATCTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACCCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAgcgaggcga                                                                                                                                             <  6:181036/149‑1 (MQ=255)          gCTGAACCACCAGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGACGCCACATACAGCGAGGCGAC‑aa                                                                                                                                         >  6:60549/1‑149 (MQ=255)          gCTGAACCACCAGTCCGCTGCGGCAATGTGCTACCGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑aa                                                                                                                                         >  7:54838/1‑149 (MQ=255)                ccaccaGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTg                                                                                                                                    >  8:90163/1‑148 (MQ=255)                ccaccaGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTg                                                                                                                                    <  7:90163/148‑1 (MQ=255)                     aGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTgc                                                                                                                              <  3:120949/149‑1 (MQ=255)                         cGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTgcggc                                                                                                                           <  3:322087/148‑1 (MQ=255)                          gCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTgcggc                                                                                                                           >  6:370763/1‑147 (MQ=255)                              cGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCTAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGt                                                                                                                      <  2:371480/148‑1 (MQ=255)                               ggCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAgg                                                                                                                       >  3:41619/1‑146 (MQ=255)                                gCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTa                                                                                                                     >  1:148653/1‑147 (MQ=255)                                        cTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGca                                                                                                           <  5:60549/149‑1 (MQ=255)                                           aCGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCCTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGcatca                                                                                                        >  1:57291/1‑149 (MQ=255)                                                      ccAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCg                                                                                             <  1:191171/149‑1 (MQ=255)                                                                 cgaAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATCGTCGAAGCCACATCCAGCGAGGGGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAAcaggt                                                                                  >  5:198583/1‑145 (MQ=255)                                                                   aaaaCTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAAcggtccc                                                                                <  8:54838/149‑1 (MQ=255)                                                                       cTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAAcggtcccggtc                                                                            >  7:208918/1‑149 (MQ=255)                                                                        tCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCTGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAAcggtcccggtc                                                                            >  2:340570/1‑148 (MQ=255)                                                                        tCAACGCCTCGTCGAACATCTCCAGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAAcggtcccggtc                                                                            >  1:89716/1‑148 (MQ=255)                                                                            cGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCt                                                                       >  5:269414/1‑149 (MQ=255)                                                                                tcgtcgAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTgccg                                                                   >  2:111649/1‑149 (MQ=255)                                                                                     gAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTc                                                                                              <  6:74034/117‑1 (MQ=255)                                                                                     gAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTc                                                                                              >  5:74034/1‑117 (MQ=255)                                                                                          tctcCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCACCTGCCGCCTGGCCCCg                                                         >  1:85318/1‑149 (MQ=255)                                                                                                   gcagcGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGcaca                                                <  5:370271/149‑1 (MQ=255)                                                                                                   gcagcGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGcaca                                                <  5:279835/149‑1 (MQ=255)                                                                                                    cagcGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCTGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACAt                                               <  1:340570/149‑1 (MQ=255)                                                                                                            gcCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAGAAATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGACCGCTTCGTTGAAACGGTCCCTGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAatcat                                        >  1:273356/1‑149 (MQ=255)                                                                                                              ccggccggCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAatcatc                                       >  1:203502/1‑147 (MQ=255)                                                                                                              ccggccggCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAatcatc                                       >  2:238634/1‑147 (MQ=255)                                                                                                                 gccggccCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACAGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGCCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACAt                                               <  7:151442/136‑1 (MQ=255)                                                                                                                 gccggccCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACAGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGCCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACAt                                               >  8:151442/1‑136 (MQ=255)                                                                                                                  ccggccCTTCCATCATGGTCGAAGCCACATCAAGCGAGGCGAG‑AATTCTCTACTGAGCCTAGGTAGAGATGTGCATCATGAACGATTCTTTCAAAAGGTCCCCTTAAATCTGCAGCCAGGCCCCGTCCTGCACATGAATCATCGAGGcc                                 >  4:240247/1‑149 (MQ=255)                                                                                                                   cggccCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCt                                >  4:323105/1‑149 (MQ=255)                                                                                                                   cggccCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCt                                >  3:149795/1‑149 (MQ=255)                                                                                                                     gccCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTgg                              <  4:41619/149‑1 (MQ=255)                                                                                                                       ccTTCCATCATGGTCGAAGCCACATCCTGCGAGGCGAG‑AATGCTGTACTGCGGCGACGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCAGGCCCCGTCCTGCACATCAATCAACGAGGCCtcg                              >  8:201755/1‑145 (MQ=255)                                                                                                                           ccATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAAc                        >  8:233483/1‑149 (MQ=255)                                                                                                                                atGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAAcgccgc                   >  7:306129/1‑149 (MQ=255)                                                                                                                                   gTCTAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGag                <  3:240247/149‑1 (MQ=255)                                                                                                                                    tCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGag                >  4:287301/1‑148 (MQ=255)                                                                                                                                      gAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGagcag             <  6:269414/149‑1 (MQ=255)                                                                                                                                         gCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTg                                                    <  1:251307/107‑1 (MQ=255)                                                                                                                                         gCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTg                                                    >  2:251307/1‑107 (MQ=255)                                                                                                                                           cacaTCCAGCGATGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGaa                         <  5:43404/132‑1 (MQ=255)                                                                                                                                           cacaTCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGaa                         >  6:43404/1‑132 (MQ=255)                                                                                                                                           cacaTCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTtgt        <  8:306129/149‑1 (MQ=255)                                                                                                                                            acaTCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTtgtg       >  3:369753/1‑149 (MQ=255)                                                                                                                                             caTCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTgtgg      <  1:238634/149‑1 (MQ=255)                                                                                                                                              aTCCTGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTgtggg     <  7:201755/149‑1 (MQ=255)                                                                                                                                              aTCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTgtggg     >  1:233571/1‑149 (MQ=255)                                                                                                                                               tCCAGCGAGGCGAG‑AATGCTGTACTGCGGAGATGTAGAGATGTGCATCATGAACTCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTgtgggt    <  8:208918/149‑1 (MQ=255)                                                                                                                                                  aGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTgtgggtgg  >  4:202993/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CTGCCAGATGCTGAACCACCAGTCCGCTGCGGCAATGTGCTCACGCAGGTTGGCCAGGGCACGACGAAAACTCAACGCCTCGTCGAACATCTCCTGCAGCAGCGAGCGCCCGGCCGGCCCTTCCATCATGGTCGAAGCCACATCCAGCGAGGCGAG‑AATGCTGTACTGCGGCGAGGTAGAGATGTGCATCATGAACGCTTCGTTGAAACGGTCCCGGTCCAGCTGCCGCCTGGCCCCGTCCTGCACATGAATCATCGAGGCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGG  >  NC_002947/4677113‑4677405

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: ‑