Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,661,522	T→A	13.8%	L51Q (CTG→CAG)	PP_1456 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,661,522	0	T	A	13.8%	77.4 / 5.8	29	L51Q (CTG→CAG)	PP_1456	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (12/13);  new base A (2/2);  total (14/15)
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.59e-01

ACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTG  >  NC_002947/1661390‑1661669                                                                                                                                     |                                                                                                                                                    aCATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggtgg                                                                                                                                    >  2:384324/1‑150 (MQ=255)      agcagcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTg                                                                                                                               <  1:80227/150‑1 (MQ=255)            gTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGTTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTg                                                                                                                         <  1:231000/150‑1 (MQ=255)                tGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCt                                                                                                                     <  1:326853/150‑1 (MQ=255)                tGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATga                                                                                                                                          >  1:454832/1‑128 (MQ=255)                tGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATga                                                                                                                                          <  2:454832/129‑2 (MQ=255)                     aTCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                                >  1:565860/1‑150 (MQ=255)                           tACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCg                                                                                                          >  2:501328/1‑150 (MQ=255)                                             ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTg                                                                                        >  1:383447/1‑150 (MQ=255)                                             ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTg                                                                                        >  1:302877/1‑150 (MQ=255)                                                 aTCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCa                                                                                    >  2:350905/1‑150 (MQ=255)                                                           catcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTACCTTCTGAGCCGGCCATCGCCCTGGCCAGCGCCACTGc                                                                          <  1:356147/150‑1 (MQ=255)                                                               tGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttc                                                                      >  1:548527/1‑150 (MQ=255)                                                                    tGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCAt                                                                 <  1:354171/150‑1 (MQ=255)                                                                            gcggcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGc                                                         <  1:350905/150‑1 (MQ=255)                                                                              ggcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTg                                                       <  2:302877/150‑1 (MQ=255)                                                                                 cTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATc                                                    <  2:565860/150‑1 (MQ=255)                                                                                       ttCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActgg                                              >  1:159379/1‑150 (MQ=255)                                                                                       ttCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATATCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActgg                                              >  1:106773/1‑150 (MQ=255)                                                                                        tCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggc                                             <  2:383447/150‑1 (MQ=255)                                                                                          aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATc                                                    <  2:501826/141‑1 (MQ=255)                                                                                          aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATc                                                    >  1:501826/1‑141 (MQ=255)                                                                                            ccTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                         <  1:384324/150‑1 (MQ=255)                                                                                                    aCATGGTGCAGACCCGTTTTGGCCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCa                                 >  1:490697/1‑150 (MQ=255)                                                                                                                        cagcATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCa                                    <  1:15751/124‑1 (MQ=255)                                                                                                                        cagcATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCa                                    >  2:15751/4‑127 (MQ=255)                                                                                                                          tCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGa           >  2:571396/1‑150 (MQ=255)                                                                                                                                   ccTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCtgtg  <  2:548527/150‑1 (MQ=255)                                                                                                                                   ccTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATATCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCtgtg  <  2:106773/150‑1 (MQ=255)                                                                                                                                     |                                                                                                                                                    ACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTG  >  NC_002947/1661390‑1661669

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