Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	419,313	C→T	8.4%	A122V (GCC→GTC)	PP_0346 →	heptose kinase WapQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	419,313	0	C	T	8.4%	205.9 / 6.0	71	A122V (GCC→GTC)	PP_0346	heptose kinase WapQ
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (3/3);  total (34/37)
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.00e+00

CACCTTCGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTACCCCAAGCAC  >  NC_002947/419173‑419449                                                                                                                                              |                                                                                                                                         cACCTTCGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTc                                                                                                                                   >  8:254467/1‑148 (MQ=255) cACCTTCGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTc                                                                                                                                   <  7:254467/148‑1 (MQ=255)  aCCTTCGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGa                                                                                                                                 >  5:330719/1‑149 (MQ=255)    cTTCGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACa                                                                                                                               <  5:69592/149‑1 (MQ=255)       cGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGGGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGc                                                                                                                             <  4:161715/148‑1 (MQ=255)          ccGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCctgc                                                                                                                         >  2:43327/1‑149 (MQ=255)            gcgAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTCGCGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCctgct                                                                                                                        >  1:69959/1‑149 (MQ=255)                    cGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGTCGACCTCGACAGCCTGCTGGCGCAGTg                                                                                                               >  2:242247/1‑149 (MQ=255)                    cGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGAc                                                                                                                                      <  1:228660/126‑1 (MQ=255)                    cGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGAc                                                                                                                                      >  2:228660/1‑126 (MQ=255)                    cGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTg                                                                                                               >  7:256156/1‑149 (MQ=255)                        aCATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCcg                                                                                                           >  1:188002/1‑149 (MQ=255)                            cAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGt                                                                                                       >  2:27360/1‑149 (MQ=255)                              gTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGc                                                                                                                                          <  1:42088/112‑1 (MQ=255)                              gTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGc                                                                                                                                          >  2:42088/1‑112 (MQ=255)                               tCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTgg                                                                                                    <  6:330719/149‑1 (MQ=255)                               tCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTgg                                                                                                    >  7:87358/1‑149 (MQ=255)                               tCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTgg                                                                                                    <  8:87358/149‑1 (MQ=255)                                 gCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGTCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGt                                                                                                  <  1:242247/149‑1 (MQ=255)                                 gCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGt                                                                                                  >  5:143313/1‑149 (MQ=255)                                     ccAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGa                                                                                                <  6:351157/147‑1 (MQ=255)                                      cAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCCTCCTGATGACCCGCGCTCTGGACGAGTGGGTCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGc                                                                                             >  4:305172/1‑149 (MQ=255)                                      cAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGTCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGc                                                                                             >  4:310164/1‑149 (MQ=255)                                      cAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGAc                                                                                                                                <  3:332281/114‑1 (MQ=255)                                      cAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGAc                                                                                                                                >  4:332281/1‑114 (MQ=255)                                        gAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGcc                                                                                            >  2:24373/1‑148 (MQ=255)                                          aGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGa                                                                                          <  8:256156/148‑1 (MQ=255)                                                 cATTC‑CGGCATTTCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGGGAAGCCGAGCGCACTg                                                                                  <  6:164486/149‑1 (MQ=255)                                                 cATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTg                                                                                  <  1:27360/149‑1 (MQ=255)                                                 cATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTg                                                                                  <  2:188002/149‑1 (MQ=255)                                                       cGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCctgctg                                                                                                                       >  2:299350/1‑107 (MQ=255)                                                       cGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCctgctg                                                                                                                       <  1:299350/107‑1 (MQ=255)                                                       cGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATc                                                                             <  3:167969/149‑1 (MQ=255)                                                        ggCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATcc                                                                            >  2:316554/1‑149 (MQ=255)                                                         gCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCt                                                                           >  4:349973/1‑149 (MQ=255)                                                            ttGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCTGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCa                                                                        <  1:43327/149‑1 (MQ=255)                                                              gCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAgg                                                                      >  1:43545/1‑149 (MQ=255)                                                               cAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGaa                                                                                               <  7:251575/123‑1 (MQ=255)                                                               cAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGaa                                                                                               >  8:251575/1‑123 (MQ=255)                                                               cAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGc                                                                     >  3:213876/1‑149 (MQ=255)                                                               cAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGc                                                                     <  6:143313/149‑1 (MQ=255)                                                                   cGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGc                                                                 <  2:69959/149‑1 (MQ=255)                                                                            aCGGCGAGCGCAAGGATGGCGGCCATCACCGGGCCATCCGTCTTCCCCGCGCTGTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCt                                                         <  5:46879/148‑1 (MQ=255)                                                                             cGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGt                                                                                                                <  2:316303/92‑1 (MQ=255)                                                                             cGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGt                                                                                                                >  1:316303/1‑92 (MQ=255)                                                                                  aGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGt                                                  <  1:316554/149‑1 (MQ=255)                                                                                    cgcAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTAc                                                <  1:24373/149‑1 (MQ=255)                                                                                      cAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACGGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGgc                                                      >  1:262386/1‑141 (MQ=255)                                                                                      cAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACCGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGgc                                                      <  2:262386/141‑1 (MQ=255)                                                                                       aaGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGTCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCt                                             <  3:305172/149‑1 (MQ=255)                                                                                       aaGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGTCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCt                                             <  3:310164/149‑1 (MQ=255)                                                                                            gggCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGa                                                                                          >  6:333264/1‑99 (MQ=255)                                                                                            gggCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGa                                                                                          <  5:333264/99‑1 (MQ=255)                                                                                            gggCGGCAAGCACCAGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGca                                          >  5:264935/1‑147 (MQ=255)                                                                                                     gCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACa                                                                                                                               >  7:143996/1‑53 (MQ=255)                                                                                                     gCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACa                                                                                                                               <  8:143996/53‑1 (MQ=255)                                                                                                             ccATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCAt                       <  2:43545/149‑1 (MQ=255)                                                                                                              cATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATg                      >  7:205393/1‑149 (MQ=255)                                                                                                              cATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATg                      <  4:213876/149‑1 (MQ=255)                                                                                                                tCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGgc                    >  4:77755/1‑149 (MQ=255)                                                                                                                   tgatgaCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTg                                                                                  <  6:79453/84‑1 (MQ=255)                                                                                                                   tgatgaCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTg                                                                                  >  5:79453/1‑84 (MQ=255)                                                                                                                     atgaCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGgctgat                <  3:349973/148‑3 (MQ=255)                                                                                                                      tgaCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTc              <  8:205393/149‑1 (MQ=255)                                                                                                                       gaCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGcc                                   <  4:247985/127‑1 (MQ=255)                                                                                                                       gaCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGcc                                   >  3:247985/1‑127 (MQ=255)                                                                                                                        aCCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTa            >  6:336653/1‑149 (MQ=255)                                                                                                                          ccGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTAcc          >  6:85178/1‑149 (MQ=255)                                                                                                                           cgcgCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGgctgc                 >  7:148910/1‑141 (MQ=255)                                                                                                                           cgcgCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGgctgc                 <  8:148910/141‑1 (MQ=255)                                                                                                                                    gACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTACCCCAAGcac  <  5:85178/147‑1 (MQ=255)                                                                                                                                         gTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGTCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTAcc          >  1:268556/1‑134 (MQ=255)                                                                                                                                         gTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTAcc          <  2:268556/134‑1 (MQ=255)                                                                                                                                              |                                                                                                                                         CACCTTCGCCCGCGAATTCCGCAACATCAGTCGCTACCAGAAGCTGCGCATTC‑CGGCATTGCAGGCGGTGTTCTACGGCGAGCGCAAGCAGGGCGGCAAGCACCGGGCCATCCTGATGACCCGCGCTCTGGACGAGTGGGCCGACCTCGACAGCCTGCTGGCGCAGTGGCCGCGGTTGGGTGAAGCCGAGCGCACTGGCATCCTGCAGGCCTGCGGTCAGCTGGCGCGTACCCTGCACAGCGCCGGCCAGGTGCATGGCTGCTTCTACCCCAAGCAC  >  NC_002947/419173‑419449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑