Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,004,526	A→T	9.0%	Q184L (CAG→CTG)	PP_0866 →	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,004,526	0	A	T	9.0%	148.6 / 4.0	53	Q184L (CAG→CTG)	PP_0866	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base T (2/3);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004382‑1004674                                                                                                                                                    |                                                                                                                                                     aTGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGc                                                                                                                                                  <  6:117094/149‑1 (MQ=255)    gCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAgt                                                                                                                                               <  2:14317/149‑1 (MQ=255)      tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCTGGGCAgtgt                                                                                                                                             >  7:177893/1‑149 (MQ=255)         cATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCTGGGCAgtgtcga                                                                                                                                          <  8:177893/149‑4 (MQ=255)          aTCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                         <  5:29131/149‑1 (MQ=255)          aTCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                         <  5:184481/149‑1 (MQ=255)          aTCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                         >  6:284279/1‑149 (MQ=255)              ggAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTca                                                                                                                                     >  2:107084/1‑149 (MQ=255)               gAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTca                                                                                                                                     <  4:105831/148‑1 (MQ=255)                  cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGGCCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAg                                                                                                                                  <  1:213804/149‑1 (MQ=255)                  cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCTGGGCAgtg                                                                                                                                              <  5:296603/136‑1 (MQ=255)                  cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCTGGGCAgtg                                                                                                                                              >  6:296603/1‑136 (MQ=255)                  cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAgg                                                                                                                                 >  5:119597/1‑149 (MQ=255)                     aaTGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAc                                                                                                                              >  7:215/1‑149 (MQ=255)                            acacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       >  7:141655/1‑149 (MQ=255)                               cGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGAc                                                                                                                    <  1:155644/149‑1 (MQ=255)                                gCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGc                                                                                                                          <  7:160739/142‑1 (MQ=255)                                gCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGc                                                                                                                          >  8:160739/1‑142 (MQ=255)                                  aGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTc                                                                                                                 >  4:125603/1‑149 (MQ=255)                                       cAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGcc                                                                                                             <  7:143550/148‑1 (MQ=255)                                        aGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTg                                                                                                           >  6:71293/1‑149 (MQ=255)                                            cGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGAcc                                                                                                       >  6:178153/1‑149 (MQ=255)                                               gTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       <  4:43962/130‑1 (MQ=255)                                               gTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       >  3:43962/1‑130 (MQ=255)                                                   aGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagcca                                                                                                 <  5:71293/148‑1 (MQ=255)                                                          gCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCgt                                                                                         <  6:119597/149‑1 (MQ=255)                                                             aGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                            <  8:214103/95‑1 (MQ=255)                                                             aGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAG‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGtgtg                                                                                                                                            >  7:214103/1‑95 (MQ=255)                                                                 gCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCgg                                                                                  >  3:141299/1‑149 (MQ=255)                                                                    cGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTgg                                                                               <  7:78401/149‑1 (MQ=255)                                                                     gTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCa                                                                                              <  1:86961/133‑1 (MQ=255)                                                                     gTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCa                                                                                              >  2:86961/1‑133 (MQ=255)                                                                                   tGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGc                                                                >  7:161809/1‑149 (MQ=255)                                                                                             aTGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                       <  5:284279/148‑1 (MQ=255)                                                                                                  cAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAg                                                  <  1:107084/148‑1 (MQ=255)                                                                                                   aGCGTCGCCCCTCCGACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCaa                                                   <  1:105783/148‑1 (MQ=255)                                                                                                      gTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGa                                             >  8:200478/1‑149 (MQ=255)                                                                                                         gCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTc                                          <  5:178153/149‑1 (MQ=255)                                                                                                           ccAT‑‑GCCAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACACGCTTTGACGTCCGCCTGGACGCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGc                                        >  4:175397/1‑149 (MQ=255)                                                                                                                   cAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGgcgc                                  <  7:200478/149‑1 (MQ=255)                                                                                                                     gTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCaa                                <  3:125603/149‑1 (MQ=255)                                                                                                                          cGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCCCCCGTGTGGCGGTGGCACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGCCTCGCCGGCGCACCTGAc                           >  6:49451/1‑149 (MQ=255)                                                                                                                            g‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACgg                         >  3:96320/1‑149 (MQ=255)                                                                                                                            g‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACgg                         >  8:213562/1‑149 (MQ=255)                                                                                                                                gTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCg                      >  6:312395/1‑149 (MQ=255)                                                                                                                                gTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCg                      >  3:85792/1‑149 (MQ=255)                                                                                                                                 tttGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCgg                     <  4:141299/149‑1 (MQ=255)                                                                                                                                            cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAgggg              >  2:46646/1‑145 (MQ=255)                                                                                                                                            cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAgggg              <  1:46646/145‑1 (MQ=255)                                                                                                                                                 cccTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCatat      <  7:207906/148‑1 (MQ=255)                                                                                                                                                 cccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata       >  3:116853/1‑147 (MQ=255)                                                                                                                                                 cccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata       >  5:196071/1‑147 (MQ=255)                                                                                                                                                   cAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGa   <  6:196071/149‑1 (MQ=255)                                                                                                                                                    aGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  <  4:96320/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCAT‑‑GACAGTAACCGG‑CCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004382‑1004674

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