Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,917,129	T→C	6.4%	L304P (CTG→CCG)	mexB →	multidrug resistance protein MexB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,917,129	0	T	C	6.4%	180.1 / 2.6	62	L304P (CTG→CCG)	mexB	multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base C (2/2);  total (30/32)
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.65e-01

GGCTCCACGGTCACCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCA  >  NC_002947/3916984‑3917274                                                                                                                                                   |                                                                                                                                                   ggCTCCACGGTCACCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTgg                                                                                                                                                  >  5:283985/1‑148 (MQ=255) ggCTCCACGGTCACCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTgg                                                                                                                                                  >  7:304781/1‑148 (MQ=255)   cTCCACGGTCACCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTAc                                                                                                                                               >  2:172514/1‑149 (MQ=255)    tCCACGGTCACCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACg                                                                                                                                              >  8:104329/1‑149 (MQ=255)             aCCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGa                                                                                                                                      <  7:104329/148‑1 (MQ=255)             aCCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGAt                                                                                                                                     >  6:29504/1‑149 (MQ=255)                aTCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCa                                                                                                                                  <  4:99383/149‑1 (MQ=255)                   ggggATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGCACGGGCGAAGATGCAgg                                                                                                                                <  5:336393/148‑1 (MQ=255)                    gggATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGAc                                                                                                                              <  2:119353/149‑1 (MQ=255)                     ggATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGAcc                                                                                                                             <  5:104476/149‑1 (MQ=255)                          tCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGgcg                                                                                                                        >  5:216010/1‑149 (MQ=255)                            gCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGgcgcg                                                                                                                      >  3:134520/1‑149 (MQ=255)                                     gAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTccc                                                                                                             >  5:196442/1‑149 (MQ=255)                                     gAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTcc                                                                                                              <  8:304781/148‑1 (MQ=255)                                      aaaTCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCg                                                                                                            >  3:129495/1‑149 (MQ=255)                                      aaaTCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCg                                                                                                            >  4:188345/1‑149 (MQ=255)                                         tCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGaa                                                                                                         <  5:29504/149‑1 (MQ=255)                                           ggTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAg                                                                                                        <  6:283985/148‑1 (MQ=255)                                               cccAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGaa                                                                                                                                        <  7:271026/112‑1 (MQ=255)                                               cccAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGaa                                                                                                                                        >  8:271026/1‑112 (MQ=255)                                                ccAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGAGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTc                                                                                                   >  1:278285/1‑148 (MQ=255)                                                 cAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTc                                                                                                   >  3:124017/1‑147 (MQ=255)                                                 cAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTc                                                                                                   >  4:25923/1‑147 (MQ=255)                                                       tACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGa                                                                                          <  3:25923/149‑1 (MQ=255)                                                          cAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAg                                                                                                        >  2:174923/1‑133 (MQ=255)                                                          cAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAg                                                                                                        <  1:174923/133‑1 (MQ=255)                                                           aGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGt                                                                                                                                                 >  2:32400/1‑91 (MQ=255)                                                            gTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATcc                                                                                     <  6:364998/149‑1 (MQ=255)                                                                 gCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTAc                                                                                <  6:196442/149‑1 (MQ=255)                                                                   aCGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGa                                                                              <  1:172514/149‑1 (MQ=255)                                                                     gcgcCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGaa                                                                                                                                        >  1:140445/1‑90 (MQ=255)                                                                     gcgcCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGaa                                                                                                                                        <  2:140445/90‑1 (MQ=255)                                                                           gAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCCGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGt                                                                                              <  8:13307/125‑1 (MQ=255)                                                                           gAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCCGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGt                                                                                              >  7:13307/1‑125 (MQ=255)                                                                            aaCGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGcc                                                                     >  1:13836/1‑149 (MQ=255)                                                                            aaCGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGcc                                                                     <  2:60763/149‑1 (MQ=255)                                                                                gCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGa                                                                                                                               <  4:92356/88‑1 (MQ=255)                                                                                gCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGa                                                                                                                               >  3:92356/1‑88 (MQ=255)                                                                                gCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTc                                                                 <  6:216010/149‑1 (MQ=255)                                                                                      cGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAg                                                           <  1:96805/149‑1 (MQ=255)                                                                                       ggCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAgg                                                          <  4:84696/149‑1 (MQ=255)                                                                                         ccACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCACAGTACGACATCCCCTACGACACCTCGCCGTTc                                                                 <  1:43683/141‑1 (MQ=255)                                                                                         ccACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTc                                                                 >  2:43683/1‑140 (MQ=255)                                                                                              gCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTTGACGAACGGCGAGGTGtcgta                                                   <  4:134520/149‑3 (MQ=255)                                                                                               ccTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTTGACGAACGGCGAGGTGtcgtag                                                  <  2:278285/149‑4 (MQ=255)                                                                                                cTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCGATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGtcgatc                                                  >  1:48486/1‑149 (MQ=255)                                                                                                  tCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCg                                                                <  7:92153/132‑1 (MQ=255)                                                                                                  tCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCg                                                                >  8:92153/1‑132 (MQ=255)                                                                                                       gTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGc                                                                                                      >  6:331179/1‑90 (MQ=255)                                                                                                          cAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATTGAGCAGGtcat                                       <  4:129495/149‑1 (MQ=255)                                                                                                               gTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAaca                                  >  6:177903/1‑149 (MQ=255)                                                                                                                 cACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACAcc                                >  3:99702/1‑149 (MQ=255)                                                                                                                 cACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACAcc                                >  5:209774/1‑149 (MQ=255)                                                                                                                          ccAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAg                       <  6:209774/149‑1 (MQ=255)                                                                                                                          ccAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAg                       <  4:99702/149‑1 (MQ=255)                                                                                                                           cAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAgg                      <  8:193586/149‑1 (MQ=255)                                                                                                                              cGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcatc                                      <  6:106783/130‑1 (MQ=255)                                                                                                                              cGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcatc                                      >  5:106783/1‑130 (MQ=255)                                                                                                                                   aTGGAAACCGCTACCCCGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGc                                              >  6:59603/1‑118 (MQ=255)                                                                                                                                   aTGGAAACCGCTACCCCGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGc                                              <  5:59603/118‑1 (MQ=255)                                                                                                                                      gAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATgctgct            >  8:121282/1‑149 (MQ=255)                                                                                                                                                cccTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCa  <  3:188345/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GGCTCCACGGTCACCATCGGGGATGTCGCCCGGGTTGAAATCGGTGCCCAGGAATACCAGTACGGCACGCGCCTGAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCC‑ATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCCCGGAAGGCGTCA‑AGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCA  >  NC_002947/3916984‑3917274

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