Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,630,918	G→A	13.7%	D114N (GAC→AAC)	lon‑II →	DNA‑binding ATP‑dependent protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,630,918	0	G	A	13.7%	65.7 / 2.8	29	D114N (GAC→AAC)	lon‑II	DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (13/12);  new base A (2/2);  total (15/14)
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.89e-01

CGGTACCGTTGCCACCGTGCTGCAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGGCGCCT  >  NC_002947/2630779‑2631060                                                                                                                                            |                                                                                                                                               cGGTACCGTTGCCACCGTGCTGCAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGa                                                                                                                                      <  2:605308/150‑1 (MQ=255)                       cAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGaa                                                                                                                >  1:477472/1‑150 (MQ=255)                        acctgctgAAGCTGCCTGCTGACACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTGCGAAGTTGAAGGGCACATCCGTGCCGATGTCTCCCTGATCGACCAAACCGATACTGCCGAGCGCGAGTCGGCAg                                                                                                               >  2:345680/3‑150 (MQ=255)                            gctgAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGtt                                                                                                           <  2:477472/150‑1 (MQ=255)                              tgAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCg                                                                                                         <  2:263649/150‑1 (MQ=255)                                    tGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGcgcg                                                                                                                        <  2:608304/129‑1 (MQ=255)                                    tGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGcgcg                                                                                                                        >  1:608304/1‑129 (MQ=255)                                                 cttgacggTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGaa                                                                                                                <  2:415474/118‑1 (MQ=255)                                                 cttgacggTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGaa                                                                                                                >  1:415474/7‑124 (MQ=255)                                                  gTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAg                                                                                     <  2:138793/150‑1 (MQ=255)                                                    cAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGtt                                                                                   <  2:352373/150‑1 (MQ=255)                                                            tGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCAACGAAACCGATGCTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTGGAACAGt                                                                           >  1:111104/1‑150 (MQ=255)                                                                      cGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCt                                                                 <  2:472073/150‑1 (MQ=255)                                                                        agcagcGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTgg                                                               <  1:193441/150‑1 (MQ=255)                                                                                             gcctgtcCGTAGTTGAAGGGCAGGTCCTTGCCGAAGTCTCCCTGATCAACGAAACCGATACTGCCGGGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAg                                          <  2:111104/145‑1 (MQ=255)                                                                                               ttCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTg                                        >  1:236842/1‑150 (MQ=255)                                                                                                                                 tctcCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGgc      <  1:56468/150‑1 (MQ=255)                                                                                                                                   tcCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGgcgc    >  1:314306/1‑150 (MQ=255)                                                                                                                                   tcCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGgcgc    >  2:482015/1‑150 (MQ=255)                                                                                                                                   tcCCTGATCAACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGaa                                                                               <  1:598093/75‑1 (MQ=255)                                                                                                                                   tcCCTGATCAACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGaa                                                                               >  2:598093/1‑75 (MQ=255)                                                                                                                                     ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTTCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCAAGCCGAAGTCCTGTCTTTGCTGAACAGCATCGACGAGCGAGGGCGCCt  >  1:299317/1‑150 (MQ=255)                                                                                                                                     ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTTAACAGCATCGACGAGCCAGGGCGCCt  >  2:494644/1‑150 (MQ=255)                                                                                                                                     ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGGCGCCt  >  1:610636/1‑150 (MQ=255)                                                                                                                                     ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCCTCGCTGAACAGCATCGACGAGCCAGGGCGCCt  >  1:412249/1‑150 (MQ=255)                                                                                                                                      cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCa                                                 >  2:206599/1‑102 (MQ=255)                                                                                                                                      cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCa                                                 <  1:206599/102‑1 (MQ=255)                                                                                                                                      cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCt                             <  2:243345/122‑1 (MQ=255)                                                                                                                                      cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCt                             >  1:243345/1‑122 (MQ=255)                                                                                                                                            |                                                                                                                                               CGGTACCGTTGCCACCGTGCTGCAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGGCGCCT  >  NC_002947/2630779‑2631060

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