Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,006,474	A→C	7.3%	K1019Q (AAG→CAG)	PP_2627 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,006,474	0	A	C	7.3%	169.6 / 4.8	55	K1019Q (AAG→CAG)	PP_2627	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/22);  new base C (2/2);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

GTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAG  >  NC_002947/3006327‑3006620                                                                                                                                                    |                                                                                                                                                   gtgtACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGaa                                                                                                                                                   <  7:354920/149‑1 (MQ=255)      cAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGAc                                                                                                                                               >  6:282820/1‑148 (MQ=255)       aGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACaa                                                                                                                                             <  3:282033/149‑1 (MQ=255)       aGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACaa                                                                                                                                             >  2:222010/1‑149 (MQ=255)         ttttttCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAAcc                                                                                                                                           >  8:341194/1‑149 (MQ=255)              tCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                       <  3:334199/148‑1 (MQ=255)               ccAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAg                                                                                                                                      <  5:282820/148‑1 (MQ=255)                cAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGc                                                                                                                                    >  7:22593/1‑149 (MQ=255)                cAGCAACGACTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTGTCCTTCAGGAAAACCTTGAga                                                                                                                                     <  8:271216/148‑2 (MQ=255)                  gCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTACCTCAAGGTTGTCCTTCAGGAAAACCTTGAgatac                                                                                                                                  <  3:15066/149‑5 (MQ=39)                      cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                              <  6:75750/149‑1 (MQ=255)                               ggCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCa                                                                                                                     <  7:341194/149‑1 (MQ=255)                                   gCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                       <  2:137836/127‑1 (MQ=255)                                   gCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGa                                                                                                                                       >  1:137836/1‑127 (MQ=255)                                                 gtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGAt                                                                                                   >  3:132374/1‑149 (MQ=255)                                                         aCCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGcc                                                                                            >  8:306584/1‑148 (MQ=255)                                                          cccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGa                                                                                          >  8:72676/1‑149 (MQ=255)                                                                  ttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTc                                                                                                 <  2:82792/134‑1 (MQ=255)                                                                  ttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTc                                                                                                 >  1:82792/1‑134 (MQ=255)                                                                  ttcttTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATc                                                                                   >  1:332882/1‑148 (MQ=255)                                                                       tGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGag                                                                                <  1:222010/146‑1 (MQ=255)                                                                          aTGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTc                                                                          >  7:233610/1‑149 (MQ=255)                                                                             aCTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                                        <  7:72676/148‑1 (MQ=255)                                                                               ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                                        >  7:150602/1‑146 (MQ=255)                                                                               ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                                        <  8:150602/146‑1 (MQ=255)                                                                               ttCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCt                                                                     >  7:84221/1‑149 (MQ=255)                                                                                     gCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGa                                                               <  2:332882/149‑1 (MQ=255)                                                                                       ggtcttcGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                             <  7:306584/147‑1 (MQ=255)                                                                                       gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCg                                                                                  >  1:341243/1‑128 (MQ=255)                                                                                       gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCg                                                                                  <  2:341243/128‑1 (MQ=255)                                                                                       gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAgc                                                             >  6:217592/1‑149 (MQ=255)                                                                                                gCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCcgc                                                     >  2:374259/1‑148 (MQ=255)                                                                                                  cGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGa                                                   >  6:194445/1‑148 (MQ=255)                                                                                                  cGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGa                                                   >  1:217357/1‑148 (MQ=255)                                                                                                   gAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGaaa                                                 >  5:79160/1‑149 (MQ=255)                                                                                                        ggCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttc                                            <  2:217357/149‑1 (MQ=255)                                                                                                          cgcCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                          <  1:374259/149‑1 (MQ=255)                                                                                                                     gTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                               <  4:132374/149‑1 (MQ=255)                                                                                                                     gTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTg                               <  8:84221/149‑1 (MQ=255)                                                                                                                          aTGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAg                          >  8:345234/1‑149 (MQ=255)                                                                                                                              tCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTaa                      <  6:79160/149‑1 (MQ=255)                                                                                                                                      tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg              >  4:230937/1‑149 (MQ=255)                                                                                                                                      tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCGTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg              >  6:269843/1‑149 (MQ=255)                                                                                                                                      tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg              >  5:307630/1‑149 (MQ=255)                                                                                                                                      tcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg              >  8:14706/1‑149 (MQ=255)                                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCa                                                          >  8:146558/1‑104 (MQ=255)                                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCa                                                          <  7:146558/104‑1 (MQ=255)                                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                          >  4:23143/1‑120 (MQ=255)                                                                                                                                       cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACcttctt                                          <  3:23143/120‑1 (MQ=255)                                                                                                                                        aaGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCa               <  8:233610/146‑1 (MQ=255)                                                                                                                                        aaGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt             >  6:272957/1‑148 (MQ=255)                                                                                                                                          ggTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCa          <  3:230937/149‑1 (MQ=255)                                                                                                                                           gTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAg         >  5:171223/1‑149 (MQ=255)                                                                                                                                                  tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAg  >  4:167194/1‑149 (MQ=255)                                                                                                                                                  tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAg  >  2:90871/1‑149 (MQ=255)                                                                                                                                                  tGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAg  >  7:160694/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGGGCAATCCCTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAGCGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTTCAGCATCGAG  >  NC_002947/3006327‑3006620

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