Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,714,136	A→G	6.8%	M529V (ATG→GTG)	ubiB →	2‑octaprenylphenol hydroxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,714,136	0	A	G	6.8%	170.8 / 2.8	59	M529V (ATG→GTG)	ubiB	2‑octaprenylphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (24/31);  new base G (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

CCACAATTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACA  >  NC_002947/5713989‑5714280                                                                                                                                                    |                                                                                                                                                  ccACAATTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGAt                                                                                                                                                  >  8:309675/1‑149 (MQ=255) acacaATTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGAt                                                                                                                                                  >  3:312467/2‑149 (MQ=255)   acaATTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATgc                                                                                                                                                >  3:120457/1‑149 (MQ=255)     aaTTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATgctg                                                                                                                                              <  4:162772/149‑1 (MQ=255)      aTTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATgctgg                                                                                                                                             <  4:120457/149‑1 (MQ=255)           cGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCg                                                                                                                                        <  7:309675/149‑1 (MQ=255)                    gccgccAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTAc                                                                                                                               <  5:319009/149‑1 (MQ=255)                           gggCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCg                                                                                                                        >  8:48624/1‑149 (MQ=255)                              cGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTgc                                                                                                                     <  7:48624/149‑1 (MQ=255)                                aCCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTgcgc                                                                                                                   >  8:143180/1‑149 (MQ=255)                                aCCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTgcgc                                                                                                                   <  4:130758/149‑1 (MQ=255)                                  cGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTa                                                                                                                                >  5:191676/1‑134 (MQ=255)                                  cGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTa                                                                                                                                <  6:191676/134‑1 (MQ=255)                                   gCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTAGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGAt                                                                                                                          <  6:26995/139‑1 (MQ=255)                                   gCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTAGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGAt                                                                                                                          >  5:26995/1‑139 (MQ=255)                                         cgcTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATagcca                                                                                                          <  6:146614/149‑1 (MQ=255)                                             gcgcCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCcg                                                                                                      >  3:99930/1‑149 (MQ=255)                                              cgcCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATgc                                                                                                                                                >  3:140111/1‑106 (MQ=255)                                              cgcCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATgc                                                                                                                                                <  4:140111/106‑1 (MQ=255)                                              cgcCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATagccagcc                                                                                                       <  8:147716/147‑1 (MQ=255)                                                ccTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCgcgc                                                                                                   <  1:108001/149‑1 (MQ=255)                                                  tgctgcGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGAGGGTGATGCTGGCTGCGGGCCTTAACCTGAACGTGCGCCAATACCCAGCCGCGCCt                                                                                                 >  6:3426/1‑149 (MQ=255)                                                            cGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGca                                                                                       >  2:109355/1‑149 (MQ=255)                                                                cTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    >  3:45305/1‑148 (MQ=255)                                                                     cggcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGGTGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    >  2:50947/1‑143 (MQ=255)                                                                     cggcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGGTGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    <  1:50947/143‑1 (MQ=255)                                                                     cggcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGGTGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    <  1:45585/143‑1 (MQ=255)                                                                     cggcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGGTGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    >  2:45585/1‑143 (MQ=255)                                                                      ggcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    <  8:263961/142‑1 (MQ=255)                                                                      ggcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCacac                                                                                    >  7:263961/1‑142 (MQ=255)                                                                       gcggTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCaa                                                                            <  6:281378/149‑1 (MQ=255)                                                                          gTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAgg                                                                         >  7:72033/1‑149 (MQ=255)                                                                               gcAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCg                                                                    >  1:67220/1‑149 (MQ=255)                                                                                   tCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCAGGCCGCCTGGCCCGCGGGGCTGATGCTGGCTGCGGAGCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGca                                                                 >  5:257960/1‑148 (MQ=255)                                                                                     ttGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑g                                                             >  1:68770/1‑149 (MQ=255)                                                                                        gCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGa                                                           <  2:42820/148‑1 (MQ=255)                                                                                          cGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCg                                                        >  6:118059/1‑149 (MQ=255)                                                                                            cAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCggg                                                      >  7:252829/1‑149 (MQ=255)                                                                                                     ccGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTg                                             <  8:252829/149‑1 (MQ=255)                                                                                                     ccGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTg                                             >  8:308408/1‑149 (MQ=255)                                                                                                     ccGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTg                                             <  2:67220/149‑1 (MQ=255)                                                                                                     ccGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTg                                             <  5:308053/149‑1 (MQ=255)                                                                                                       gAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGa                                           <  4:99930/149‑1 (MQ=255)                                                                                                       gAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGa                                           <  2:68770/149‑1 (MQ=255)                                                                                                        aaaCCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAg                                          <  7:308408/149‑1 (MQ=255)                                                                                                          aCCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTc                                        <  4:45305/149‑1 (MQ=255)                                                                                                             gCCAGCCTGGCCGCCCCTGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTcgac                                     <  6:280841/149‑1 (MQ=255)                                                                                                             gCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTcgac                                     >  5:35162/1‑149 (MQ=255)                                                                                                                       ccgccCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGAc                           <  5:118059/149‑1 (MQ=255)                                                                                                                             gggCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAATGGGCCGGGCACA‑GGCGCGGGCGCGCTTTGGAGTCGACGATGAAAGActggctg                     <  4:312467/148‑1 (MQ=255)                                                                                                                                 cgccTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGAc                  >  2:272081/1‑148 (MQ=255)                                                                                                                                 cgccTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACg                 <  1:109355/149‑1 (MQ=255)                                                                                                                                       gCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACATGGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATc            <  1:272081/149‑1 (MQ=255)                                                                                                                                          cgcgTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGt        <  8:21183/149‑1 (MQ=255)                                                                                                                                                gatGATGCTGGATGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGaa    >  8:207430/3‑147 (MQ=255)                                                                                                                                                gCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACa  <  6:35162/149‑1 (MQ=255)                                                                                                                                                gCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACa  <  7:207430/149‑1 (MQ=255)                                                                                                                                                 cTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAg         >  2:328776/1‑141 (MQ=255)                                                                                                                                                 cTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAg         <  1:328776/141‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCACAATTGCCGGAACGGCGCCGCCAGGGCGACCGCTGGGCGCTGCGCCTGCTGGGCGCCGGCCTGCTCGGCGGTGGCGCAGTCTTGGCCGCAGGCGCTGCCGAAACCGCCAGCCTGGCCGCCCCGGCCGCCTGGCCCGCGTGGCTGATGCTGGCTGCGGGCCTTTACCTGATCGTGCGCCAATAGCCAGCCGCGCCTGTGGCTGGCACACTAGCGCAAAGGGCCCGGCACA‑GGAGCGGGCCCGCTTTGGAGTCGACGATGAAAGACTGGCTGGACGAGATCAAGTGGAACA  >  NC_002947/5713989‑5714280

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