Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,678,243	C→G	8.7%	A111P (GCC→CCC)	ldcC ←	lysine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,678,243	0	C	G	8.7%	140.2 / 5.4	46	A111P (GCC→CCC)	ldcC	lysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base G (2/2);  total (21/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 = 9.84e-01

TGCTGTACCAGCGCCTTGAAGAACGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCAGCGCCTTCGGCCGCGATGAGCA  >  NC_002947/4678101‑4678389                                                                                                                                               |                                                                                                                                                   tgatGTACCAGCGCCTTGAAGAACGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACTCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTcgg                                                                                                                                               <  7:153787/145‑1 (MQ=255)        ccAGCGCCTTGAAGAACGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCa                                                                                                                                       >  5:115610/1‑149 (MQ=255)        ccAGCGCCTTGAAGAACGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCa                                                                                                                                       >  3:116747/1‑149 (MQ=255)                  gaagaaCGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAgg                                                                                                                              <  4:217280/148‑1 (MQ=255)                  gaagaaCGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAgg                                                                                                                              <  4:84406/148‑1 (MQ=255)                    agaaCGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTc                                                                                                                           <  5:98183/149‑1 (MQ=255)                      aaCGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCAc                                                                                                                         <  7:218277/149‑1 (MQ=255)                               agcagGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGcc                                                                                                                <  5:126461/149‑1 (MQ=255)                                   ggCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAgg                                                                                                            >  3:99797/1‑149 (MQ=255)                                       gTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCg                                                                                                         >  7:296311/1‑148 (MQ=255)                                        tCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGGTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCg                                                                                                         >  3:24576/1‑147 (MQ=255)                                        tCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGGTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCg                                                                                                         >  3:20434/1‑147 (MQ=255)                                        tCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGaa                                                                                                       <  4:29558/149‑1 (MQ=255)                                             gtaggtGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAAGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCg                                                                                                  <  5:164992/149‑1 (MQ=255)                                                 gtgtgtGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCAcc                                                                                                                        <  4:146175/123‑1 (MQ=255)                                                 gtgtgtGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCAcc                                                                                                                        >  3:146175/1‑123 (MQ=255)                                                     gtgCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTg                                                                                          >  3:154713/1‑149 (MQ=255)                                                        cagcaCGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGGTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGt                                                                                       <  4:20434/149‑1 (MQ=255)                                                        cagcaCGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGGTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGt                                                                                       <  4:24576/149‑1 (MQ=255)                                                        cagcaCGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGt                                                                                       <  2:170551/149‑1 (MQ=255)                                                                    ccTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACgg                                                                           <  4:116747/149‑1 (MQ=255)                                                                    ccTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACgg                                                                           >  8:120427/1‑149 (MQ=255)                                                                              gagaAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGa                                                                 >  8:189146/1‑149 (MQ=255)                                                                                  aaCGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAg                                                             <  6:115610/149‑1 (MQ=255)                                                                                     ggCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTc                                                          >  8:81863/1‑149 (MQ=255)                                                                                     ggCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTc                                                          <  8:296311/149‑1 (MQ=255)                                                                                      gCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCg                                                         <  4:154713/149‑1 (MQ=255)                                                                                                ttCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCt                                               >  4:274188/1‑149 (MQ=255)                                                                                                ttCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCt                                               >  6:229826/1‑149 (MQ=255)                                                                                                ttCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGAACGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCt                                               >  1:259330/1‑149 (MQ=255)                                                                                                   gAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCa                                            <  4:99797/149‑1 (MQ=255)                                                                                                       aGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAg                                                             >  8:128573/1‑128 (MQ=255)                                                                                                       aGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAg                                                             <  7:128573/128‑1 (MQ=255)                                                                                                         gTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAAt                                      <  7:81863/149‑1 (MQ=255)                                                                                                         gTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGGTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTaaag                                      >  4:16431/1‑148 (MQ=255)                                                                                                            caGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGCGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATgg                                    <  4:179365/147‑1 (MQ=255)                                                                                                            aaGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATgg                                    <  7:120427/148‑1 (MQ=255)                                                                                                                    cACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACgg                                                                           >  4:115955/1‑101 (MQ=255)                                                                                                                    cACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACgg                                                                           <  3:115955/101‑1 (MQ=255)                                                                                                                         aGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCAgc                      <  7:189146/149‑1 (MQ=255)                                                                                                                                       gctcATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCAGCGCCTTCGGCCGCGa        >  1:204030/1‑149 (MQ=255)                                                                                                                                           aTGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCAGCGCCTTCGGCCGCGAt       >  1:286747/1‑146 (MQ=255)                                                                                                                                            tGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCAGCGCCTTCGGCCGCGATGAGc   <  2:204030/149‑1 (MQ=255)                                                                                                                                             ggCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGTCGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCCGCGCCTTCGGCCGCGATGAGCa  >  6:300897/1‑149 (MQ=255)                                                                                                                                             ggCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTg                                              >  8:115793/1‑105 (MQ=255)                                                                                                                                             ggCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTg                                              <  7:115793/105‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   TGCTGTACCAGCGCCTTGAAGAACGGCGGCAGCAGGCCGTCGAGGTAGGTGTGTGCAGCACGCGCCACCTGGCGGGCGAGAAACGGCACGGTGTCTTCGAACAGGTAAAGGATGCCACGCAGCTGGTTGAGCTCGCTCATGGCTTCGGCGGGGGCATTTTCCAGGGTCACCTGCTCGCCCAGGGCGAAGATCGGCAGGTTGGGTGCACGCAGACGGGCCAGGCGGATCAGCTCGGCCATGTTCTGCAGTAAATGGGTATTCTCGCCAGCGCCTTCGGCCGCGATGAGCA  >  NC_002947/4678101‑4678389

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