Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,004,526	A→T	7.5%	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	7.5%	189.9 / 5.2	67	Q184L (CAG→CTG)	PP_0866	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base T (3/2);  total (36/31)
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.73e-01

GGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATA  >  NC_002947/1004380‑1004671                                                                                                                                                    |                                                                                                                                                  ggATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAgg                                                                                                                                                 >  3:360971/1‑149 (MQ=255)       cTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgtg                                                                                                                                           <  7:159916/149‑1 (MQ=255)        tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgtgt                                                                                                                                          >  1:106243/1‑149 (MQ=255)        tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgtgt                                                                                                                                          >  3:344621/1‑149 (MQ=255)          gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgt                                                                                                                                          <  1:308950/147‑1 (MQ=255)            aTCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                      >  4:271991/1‑149 (MQ=255)                ggAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGc                                                                                                                                               <  8:254031/136‑1 (MQ=255)                ggAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGc                                                                                                                                               >  7:254031/1‑136 (MQ=255)                 gAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTca                                                                                                                                  >  3:310143/1‑148 (MQ=255)                    cTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAgg                                                                                                                              >  4:364377/1‑149 (MQ=255)                     tCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgt                                                                                                                                            >  2:144495/1‑134 (MQ=255)                     tCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTcac                                                                                                                                 >  5:317679/1‑145 (MQ=255)                     tCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTcac                                                                                                                                 <  6:317679/145‑1 (MQ=255)                       aaTGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAc                                                                                                                           >  6:23272/1‑149 (MQ=255)                            cAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCt                                                                                                                      >  3:37040/1‑149 (MQ=255)                              acacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                      >  2:154928/1‑131 (MQ=255)                              acacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATc                                                                                                                                      <  1:154928/131‑1 (MQ=255)                               cacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                    >  5:386302/1‑148 (MQ=255)                               cacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTg                                                                                                                   <  5:165583/149‑1 (MQ=255)                                acGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGa                                                                                                                  <  4:310143/149‑1 (MQ=255)                                     gTTGCAGGTCGAGTTCAGCGCCGTCCAGCGGCAC‑GTGCGCCTGAACGCTGGGGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTcc                                                                                                             <  4:342387/149‑1 (MQ=255)                                     gTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTcc                                                                                                             <  4:336890/149‑1 (MQ=255)                                      ttGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGCTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCg                                                                                                            >  6:69294/1‑149 (MQ=255)                                            gTCGAGTTCAGCGCCGGCCAGCGGCACGGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTg                                                                                                        >  1:209620/1‑148 (MQ=255)                                            gTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGa                                                                                                      <  7:306608/149‑1 (MQ=255)                                            gTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGa                                                                                                      <  5:23272/149‑1 (MQ=255)                                                   tCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagc                                                                                                >  7:130009/1‑148 (MQ=255)                                                     aGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCa                                                                                                                                                   >  2:176871/1‑95 (MQ=255)                                                     aGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCa                                                                                                                                                   <  1:176871/95‑1 (MQ=255)                                                     aGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagccag                                                                                             >  6:368855/1‑149 (MQ=255)                                                     aGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagcca                                                                                              <  4:152959/148‑1 (MQ=255)                                                                    cAC‑GTGCGCCTGAACGCTTGTGAAGTGATGGTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGGCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                              <  7:347831/149‑1 (MQ=255)                                                                    cAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                              <  7:318276/149‑1 (MQ=255)                                                                    cAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                              <  8:359398/149‑1 (MQ=255)                                                                    cAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                              <  3:271991/149‑1 (MQ=255)                                                                    cAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                              >  2:164024/1‑149 (MQ=255)                                                                     aC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACAGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTg                                                                             >  2:229320/1‑149 (MQ=255)                                                                            gcCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGCTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAg                                                                       <  5:69294/149‑1 (MQ=255)                                                                             ccTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTcac                                                                                                                                 <  4:360432/90‑1 (MQ=255)                                                                             ccTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTcac                                                                                                                                 >  3:360432/1‑90 (MQ=255)                                                                                             gTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTg                                                        >  3:370815/1‑147 (MQ=255)                                                                                                aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAATAGGGGTCGGTGCGGGTGCAgg                                                    <  6:386302/148‑1 (MQ=255)                                                                                                aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAggg                                                   >  8:24135/1‑149 (MQ=255)                                                                                                aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCAGCCTGGACCCAGCCAGCACCCGAGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAggg                                                   >  1:151077/1‑149 (MQ=255)                                                                                                        cGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTg                                                                             <  1:349573/115‑1 (MQ=255)                                                                                                        cGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATAGTCACAGGTAACCGCGTTGACGTACGCCTGGACCCAGCCCGCACCCGTGTGGCGGTg                                                                             >  2:349573/1‑115 (MQ=255)                                                                                                            gCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTc                                       <  5:368855/149‑1 (MQ=255)                                                                                                            gCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTc                                       <  1:164024/149‑1 (MQ=255)                                                                                                            gCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTc                                       <  2:209620/149‑1 (MQ=255)                                                                                                              ccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTcaca                                                                                                                                >  4:255031/1‑58 (MQ=255)                                                                                                              ccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTcaca                                                                                                                                <  3:255031/58‑1 (MQ=255)                                                                                                                aTGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCg                                   >  4:101676/1‑149 (MQ=255)                                                                                                                         aCCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACt                           >  5:20850/1‑148 (MQ=255)                                                                                                                            ggCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                           <  5:355991/113‑1 (MQ=255)                                                                                                                            ggCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                           >  6:355991/1‑113 (MQ=255)                                                                                                                            ggCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                       >  2:72843/1‑149 (MQ=255)                                                                                                                                   tgtgGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCggg                 >  4:240226/1‑148 (MQ=255)                                                                                                                                     tgGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCa               <  8:130009/148‑1 (MQ=255)                                                                                                                                     acgatcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGgcg                                >  5:191271/7‑131 (MQ=255)                                                                                                                                     acgatcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGgcg                                <  6:191271/125‑1 (MQ=255)                                                                                                                                         tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGt          >  4:137374/3‑149 (MQ=255)                                                                                                                                          gACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGACAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTc         <  1:263683/149‑1 (MQ=255)                                                                                                                                           acacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCg        >  3:230582/1‑149 (MQ=255)                                                                                                                                            cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGc       <  3:101676/149‑1 (MQ=255)                                                                                                                                            cacTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGc       <  4:370815/149‑1 (MQ=255)                                                                                                                                               tGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata    >  1:200841/1‑149 (MQ=255)                                                                                                                                                 cccAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCAtata  <  4:344621/149‑1 (MQ=255)                                                                                                                                                  ccAGGGCAGTGTGATAGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata    >  8:23484/1‑146 (MQ=255)                                                                                                                                                   cAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata    >  8:50704/1‑145 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GGATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCAC‑GTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATA  >  NC_002947/1004380‑1004671

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