Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,796,131	A→C	5.1%	L108R (CTG→CGG)	fcs ←	feruloyl‑CoA‑synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,796,131	0	A	C	5.1%	246.7 / 3.9	78	L108R (CTG→CGG)	fcs	feruloyl‑CoA‑synthetase
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base C (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GTGAGCACCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCT  >  NC_002947/3795984‑3796272                                                                                                                                                    |                                                                                                                                              gTGAGCACCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAg                                                                                                                                              <  5:171443/149‑1 (MQ=255)   gAGCACCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCg                                                                                                                                            >  5:390026/1‑149 (MQ=255)    aGCACCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCgg                                                                                                                                           <  3:319645/149‑1 (MQ=255)      cACCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCggg                                                                                                                                          <  5:101164/148‑1 (MQ=255)       aCCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGgcg                                                                                                                                        >  1:248420/1‑149 (MQ=255)            gCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCCGCTGGCGCtgct                                                                                                                                   <  7:38956/149‑4 (MQ=255)              aGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCCGCTGGCGCtgcttt                                                                                                                                 <  7:249157/149‑6 (MQ=255)                 cATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACt                                                                                                                               >  6:188574/1‑148 (MQ=255)                   tGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAg                                                                                                                            <  4:9384/149‑1 (MQ=255)                    ggCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAgg                                                                                                                           >  4:84429/1‑149 (MQ=255)                      cgcAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGcc                                                                                                                         >  6:261135/1‑149 (MQ=255)                       gcAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGcc                                                                                                                         >  3:161158/1‑148 (MQ=255)                        cAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCa                                                                                                                       <  5:188574/149‑1 (MQ=255)                            ttGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTcc                                                                                                                   <  7:285322/149‑1 (MQ=255)                             tGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTcc                                                                                                                   >  6:162742/1‑148 (MQ=255)                               gCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTa                                                                                                                 <  4:263670/148‑1 (MQ=255)                               gCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAg                                                                                                                <  2:74843/149‑1 (MQ=255)                                 gAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGc                                                                                                               >  8:243531/1‑148 (MQ=255)                                 gAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGc                                                                                                               >  7:258674/1‑148 (MQ=255)                                 gAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGc                                                                                                               >  1:230879/1‑148 (MQ=255)                                       ttGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTgg                                                                                                        >  6:158235/1‑149 (MQ=255)                                       ttGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTgg                                                                                                        >  5:129650/1‑149 (MQ=255)                                        gacGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTACCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTggc                                                                                                       <  6:394207/147‑1 (MQ=255)                                           gACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAGAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTggcggc                                                                                                    >  2:184645/1‑149 (MQ=255)                                              aacaGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGAt                                                                                                 >  6:157604/1‑149 (MQ=255)                                                caGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGAGTTGCAGTTGTTCGATGTCGTTGCCGGAAAGCAGAGCCAGCGGGAGCTCGGCACTGAGGCCCAGTCCAAGCAAGGTGGCGGCGATgg                                                                                               >  8:322367/1‑149 (MQ=255)                                                 aGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGt                                                                                              >  4:66892/1‑149 (MQ=255)                                                   cgcgTAGGGCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTTTTCGATCTCGTTGCCTGAAAGCAGCGCGAGCGGGCGCTCGGCACTGAGGCCCAGTCCTATCAAGTTGTCGGCTATGGt                                                                                              <  1:198914/147‑1 (MQ=255)                                                     cgTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTgcgc                                                                                          <  5:50809/149‑1 (MQ=255)                                                     cgTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTgcgc                                                                                          <  3:66892/149‑1 (MQ=255)                                                      gTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGca                                                                                         >  4:390784/1‑149 (MQ=255)                                                          gCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATc                                                                                     >  5:39779/1‑149 (MQ=255)                                                          gCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATc                                                                                     >  6:397686/1‑149 (MQ=255)                                                          gCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATc                                                                                     >  5:34210/1‑149 (MQ=255)                                                           ccGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCg                                                                                    <  1:68046/149‑1 (MQ=255)                                                             ggCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGc                                                                                  <  2:230879/149‑1 (MQ=255)                                                               cGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGa                                                                                <  3:84429/149‑1 (MQ=255)                                                               cGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGa                                                                                <  7:235206/149‑1 (MQ=255)                                                                 acacCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGc                                                                              <  5:162742/149‑1 (MQ=255)                                                                     cGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCt                                                                          <  2:129615/149‑1 (MQ=255)                                                                     cGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCt                                                                          <  6:41394/149‑1 (MQ=255)                                                                          aataggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACg                                                                     >  7:362528/1‑149 (MQ=255)                                                                           ataggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGt                                                                    <  5:261135/149‑1 (MQ=255)                                                                            taggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACg                                                                     >  7:327124/1‑147 (MQ=255)                                                                                  ataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGAt                                                             <  5:158235/149‑1 (MQ=255)                                                                                   taCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATc                                                            <  3:390784/149‑1 (MQ=255)                                                                                      cGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGaa                                                         >  2:201330/1‑149 (MQ=255)                                                                                        gCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTCCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACg                                                       >  3:328321/1‑149 (MQ=255)                                                                                         cATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAGAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGc                                                      <  1:184645/149‑1 (MQ=255)                                                                                              aTGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGc                                                 <  7:322367/149‑1 (MQ=255)                                                                                               tGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCa                                                >  2:358392/1‑149 (MQ=255)                                                                                                   gcCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGACACCGt                                            >  4:391719/1‑149 (MQ=255)                                                                                                    ccGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTc                                           >  6:11498/1‑149 (MQ=255)                                                                                                     cGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCt                                          >  8:152270/1‑149 (MQ=255)                                                                                                       aGGGCGATTTGCAGGTGTTCGATGTCGTTGCCTGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGc                                        <  7:227360/149‑1 (MQ=255)                                                                                                         ggCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCt                                      <  4:153972/149‑1 (MQ=255)                                                                                                         ggCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCt                                      >  2:179538/1‑149 (MQ=255)                                                                                                           cGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTgg                                    <  5:397686/149‑1 (MQ=255)                                                                                                            gATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGc                                   <  1:179538/149‑1 (MQ=255)                                                                                                            gATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGc                                   <  6:129650/149‑1 (MQ=255)                                                                                                                  cAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCcg                             >  3:228426/1‑149 (MQ=255)                                                                                                                   aGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCcgc                            <  3:112365/149‑1 (MQ=255)                                                                                                                   aGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCcgc                            >  1:351431/1‑149 (MQ=255)                                                                                                                   aGGTGGTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCCGGCACCGTCTGCCTGGCGGGCcgc                            >  4:55130/1‑149 (MQ=255)                                                                                                                     gtgtTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGa                          >  5:353935/1‑149 (MQ=255)                                                                                                                      tgtTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGAt                         <  7:243531/149‑1 (MQ=255)                                                                                                                       gtTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCATGCACCGTCTGCCTGGCGTGCCGCGAt                         >  2:77512/1‑148 (MQ=255)                                                                                                                           gATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAg                    >  8:16665/1‑149 (MQ=255)                                                                                                                           gATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGCCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAg                    >  1:278169/1‑149 (MQ=255)                                                                                                                            aTGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgt                   <  4:161158/149‑1 (MQ=255)                                                                                                                             tGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgtg                  <  2:68635/149‑1 (MQ=255)                                                                                                                               tCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgtggt                <  4:160452/149‑1 (MQ=255)                                                                                                                                cGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgtggtg               >  4:42763/1‑149 (MQ=255)                                                                                                                                    gCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCg           >  7:50094/1‑149 (MQ=255)                                                                                                                                           agcagcGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCAc    <  1:358392/149‑1 (MQ=255)                                                                                                                                            cgagcGCCCGCTGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCAcc   >  4:365145/3‑149 (MQ=255)                                                                                                                                            cgagcGCCCGCTGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCAcc   >  8:250602/3‑149 (MQ=255)                                                                                                                                             cagcGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCt  <  4:328321/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              GTGAGCACCTCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCT  >  NC_002947/3795984‑3796272

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