Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,220,452	T→A	6.7%	D187V (GAC→GTC)	PP_5539 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,220,452	0	T	A	6.7%	176.4 / 2.5	60	D187V (GAC→GTC)	PP_5539	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base A (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TAGGCGAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCT  >  NC_002947/3220306‑3220587                                                                                                                                                   |                                                                                                                                          tAGGCGAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTc                                                                                                                                          <  3:261221/148‑1 (MQ=255) tAGGCGAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCa                                                                                                                                         >  2:176866/1‑149 (MQ=255)    gCGAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATc                                                                                                                                       >  1:30561/1‑148 (MQ=255)    gCGAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATc                                                                                                                                       >  1:56086/1‑148 (MQ=255)      gAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCAt                                                                                                                                    <  1:176866/149‑1 (MQ=255)          ccccTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCg                                                                                                                                <  1:338918/149‑1 (MQ=255)           cccTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCg                                                                                                                                >  4:114237/1‑148 (MQ=255)           cccTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCg                                                                                                                                >  2:241550/1‑148 (MQ=255)            ccTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGaa                                                                                                                              >  8:147698/1‑149 (MQ=255)             cTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATg                                                                                                                                   >  3:79246/1‑143 (MQ=255)             cTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATg                                                                                                                                   <  4:79246/143‑1 (MQ=255)                 ggCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGt                                                                                                                         <  4:44942/149‑1 (MQ=255)                           gAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGACATCCATGGCGAACAAGTTGCGGTTggc                                                                                                               >  1:14684/1‑149 (MQ=255)                                gcTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCt                                                                                                          <  3:114237/149‑1 (MQ=255)                                 cTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTg                                                                                                         <  3:55695/149‑1 (MQ=255)                                   gAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGgc                                                                                                       >  4:46658/1‑149 (MQ=255)                                               gATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAATATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCt                                                                                           >  5:114114/1‑149 (MQ=255)                                               gATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGACATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCt                                                                                           >  7:42927/1‑149 (MQ=255)                                                  gTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCAGgccgccaa                                                                                                 <  4:36437/140‑3 (MQ=255)                                                  gTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCAGgccgccaa                                                                                                 >  3:36437/1‑138 (MQ=255)                                                       cTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTc                                                                                   <  1:380361/149‑1 (MQ=255)                                                        ttCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCa                                                                                  >  6:125617/1‑149 (MQ=255)                                                               tACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGcc                                                                           >  8:174512/1‑149 (MQ=255)                                                                   cGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATg                                                                       >  8:229395/1‑149 (MQ=255)                                                                      gaCGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCg                                                                    <  7:174512/149‑1 (MQ=255)                                                                         gCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGccgatgccgatg                                                                 >  2:205685/1‑149 (MQ=255)                                                                             cGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGACATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCa                                                            <  2:14684/149‑1 (MQ=255)                                                                               tGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACg                                                          <  2:51690/149‑1 (MQ=255)                                                                                  cGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgtg                                                       >  6:37383/1‑149 (MQ=255)                                                                                  cGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgtg                                                       <  3:46658/149‑1 (MQ=255)                                                                                   gAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTGTTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgtg‑t                                                     <  1:316735/149‑1 (MQ=255)                                                                                     aGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAACTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCATCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgtg                                                       >  4:307991/1‑146 (MQ=255)                                                                                      gTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCATGGAGCCGATGCCGGTG‑GTCACGGTG‑Taga                                                  <  4:311857/149‑1 (MQ=255)                                                                                       tCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑Tagag                                                 <  2:109696/149‑1 (MQ=255)                                                                                         aTAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAggg                                               >  6:281040/1‑149 (MQ=255)                                                                                          tAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGGCGGTG‑GTCACGgtg                                                       >  3:83686/1‑141 (MQ=255)                                                                                          tAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgtg                                                       <  4:83686/141‑1 (MQ=255)                                                                                          tAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgt                                                        <  1:126184/140‑1 (MQ=255)                                                                                          tAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGgt                                                        >  2:126184/1‑140 (MQ=255)                                                                                             gTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATc                                           >  2:184432/1‑149 (MQ=255)                                                                                                tCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGAG‑GTCACGGTG‑TAGAGGGCATCgg                                         >  5:141188/1‑148 (MQ=255)                                                                                                   aGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCAcc                                     >  8:167548/1‑149 (MQ=255)                                                                                                   aGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCAcc                                     <  2:56086/149‑1 (MQ=255)                                                                                                     gCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCtt                                   >  1:304718/1‑149 (MQ=255)                                                                                                                ccgcGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCt                         >  2:71764/1‑148 (MQ=255)                                                                                                                ccgcGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTc                        <  5:281040/149‑1 (MQ=255)                                                                                                                    gAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCt                    <  7:167548/149‑1 (MQ=255)                                                                                                                    gAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCAAGGTGCTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGc                     <  1:241550/149‑1 (MQ=255)                                                                                                                       gATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGa                 <  2:30561/149‑1 (MQ=255)                                                                                                                         tcttctTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGAc                <  7:229395/148‑1 (MQ=255)                                                                                                                          cttcttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCt              >  5:81468/1‑149 (MQ=255)                                                                                                                           ttcttGTAGTCGATGTCGTACAAGTCATCCATGGCGCACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTc                                  <  3:41715/128‑1 (MQ=255)                                                                                                                           ttcttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTc                                  >  4:41715/1‑128 (MQ=255)                                                                                                                            tcttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCggtg                                                                 <  6:285773/98‑1 (MQ=255)                                                                                                                            tcttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCggtg                                                                 >  5:285773/1‑98 (MQ=255)                                                                                                                             cttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAgg                                                <  8:303114/112‑1 (MQ=255)                                                                                                                             cttGTAGTCGATGTCGTACAAGACATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGACACCTTCGTAGTTGCTCAGCTGGACCTg             <  8:42927/147‑1 (MQ=255)                                                                                                                                   gTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGa     >  3:365643/1‑149 (MQ=255)                                                                                                                                      gATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGAGTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGAc    <  1:108222/148‑1 (MQ=255)                                                                                                                                      gATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCt  >  3:13410/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                          TAGGCGAGCCCCCTGGGGCGACGACCGAGGCGCTGAGGTTCGTGTCGATGTCCGCTTCGGCGTACCCGAGACGCAGCGTGCCGAAGTCATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTG‑GTCACGGTG‑TAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCT  >  NC_002947/3220306‑3220587

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