Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	593,434	A→C	6.7%	E391A (GAG→GCG) ‡	PP_0503 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	593,434	0	A	C	6.7%	224.5 / 6.6	74	E391A (GAG→GCG) ‡	PP_0503	MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base C (3/2);  total (41/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 = 8.47e-01

CTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAC  >  NC_002947/593288‑593579                                                                                                                                                    |                                                                                                                                                   cTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAg                                                                                                                                                   <  6:71579/148‑1 (MQ=255)  tATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgc                                                                                                                                                 >  3:17553/1‑149 (MQ=255)       gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                            >  8:5643/1‑149 (MQ=255)       gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACCGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                            >  2:146745/1‑149 (MQ=255)        cTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCt                                                                                                                                           >  4:193338/1‑149 (MQ=255)           ggCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGc                                                                                                                                         >  6:127845/1‑148 (MQ=255)                 ggTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTc                                                                                                                                  <  1:290822/149‑1 (MQ=255)                   tttGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACCGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                                                                <  1:146745/149‑1 (MQ=255)                      gggggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGt                                                                                                                              <  6:220461/148‑1 (MQ=255)                         ggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAg                                                                                                                           >  2:65185/1‑148 (MQ=255)                              tgtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTg                                                                                                                     >  4:142622/1‑149 (MQ=255)                                   gCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGcc                                                                                                                <  5:319859/149‑1 (MQ=255)                                   gCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGc                                                                                                                 >  5:250814/1‑148 (MQ=255)                                     gAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGcc                                                                                                                >  6:297641/1‑147 (MQ=255)                                          gggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgacc                                                                                                          >  1:296498/1‑148 (MQ=255)                                           ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccga                                                                                                        >  5:49993/1‑149 (MQ=255)                                           ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccga                                                                                                        >  4:234680/1‑149 (MQ=255)                                             tgtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccgac                                                                                                       >  1:183862/1‑148 (MQ=255)                                              gtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCa                                                                                                     >  2:239681/1‑149 (MQ=255)                                               tttGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCAt                                                                                                    >  5:133704/1‑149 (MQ=255)                                                 tGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATg                                                                                                   <  6:250814/148‑1 (MQ=255)                                                   cacaCTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGtat                                                                                                <  1:65185/149‑1 (MQ=255)                                                   cacaCTGGGGCTGGGGTGCGGTGGGTGGGCTTGGCAGGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGtat                                                                                                <  4:17553/149‑1 (MQ=255)                                                     cacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAc                                                                                              <  3:234680/149‑1 (MQ=255)                                                              tGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                     >  4:257097/1‑149 (MQ=255)                                                              tGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                     >  3:54548/1‑149 (MQ=255)                                                                 gTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgac                                                                                  <  6:22123/149‑1 (MQ=255)                                                                  ttGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacg                                                                                 >  1:210161/1‑149 (MQ=255)                                                                     cGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAAAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtc                                                                               >  8:157803/1‑148 (MQ=255)                                                                       gtgtgtGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                            >  6:310716/1‑149 (MQ=255)                                                                        tgtgtgGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacg                                                                           <  3:257097/149‑1 (MQ=255)                                                                          tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCAACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACAATGTATACCCACGGgtcgacgtcgacgt                                                                          >  1:213834/1‑149 (MQ=255)                                                                            tgGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcga                                                                       >  4:164876/1‑149 (MQ=255)                                                                              ggCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAcc                                                                     >  1:287972/1‑149 (MQ=255)                                                                               gCTTGGCATGGGCATCTGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCt                                                                    <  7:5643/149‑1 (MQ=255)                                                                                 ttGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTg                                                                  >  3:253877/1‑149 (MQ=255)                                                                                 ttGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTg                                                                  <  6:133704/149‑1 (MQ=255)                                                                                         ggCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGc                                                           <  2:213834/148‑1 (MQ=255)                                                                                         ggCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGc                                                           <  5:127845/148‑1 (MQ=255)                                                                                           cATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGt                                                        <  5:310716/149‑1 (MQ=255)                                                                                             tCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg                                                      <  2:210161/149‑1 (MQ=255)                                                                                             acAGTGCGGCGGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg                                                      >  1:35433/2‑149 (MQ=255)                                                                                                   ccgccgC‑ACTGTGCTACTGGCTATGGCTGCACCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGGCGACCTTGGTCTTGCTGTAGCC‑AAgg                                               <  1:220703/149‑1 (MQ=255)                                                                                                   ccgccgC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAGTGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAgg                                               >  3:43261/1‑149 (MQ=255)                                                                                                   ccgccgA‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCATAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAgg                                               <  3:193338/149‑1 (MQ=255)                                                                                                    cgccgc‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAggg                                              >  7:59265/1‑149 (MQ=255)                                                                                                     gccgc‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGt                                                                                                  <  1:112451/97‑1 (MQ=255)                                                                                                     gccgc‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGt                                                                                                  >  2:112451/1‑97 (MQ=255)                                                                                                        gc‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGt                                          <  3:164876/149‑1 (MQ=255)                                                                                                         c‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTc                                         >  8:119063/1‑149 (MQ=255)                                                                                                            cTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCgg                                       <  2:35433/149‑1 (MQ=255)                                                                                                             tgtgCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCggg                                      >  6:63201/1‑149 (MQ=255)                                                                                                                     tGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGAc                              >  5:137467/1‑149 (MQ=255)                                                                                                                             gCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCt                                                          >  3:259136/1‑114 (MQ=255)                                                                                                                             gCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCt                                                          <  4:259136/114‑1 (MQ=255)                                                                                                                              cTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCtt                     <  2:283775/149‑1 (MQ=255)                                                                                                                               tGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                    >  2:336903/1‑149 (MQ=255)                                                                                                                               tGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                    >  4:60578/1‑149 (MQ=255)                                                                                                                                 cAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgt                  <  4:54548/149‑1 (MQ=255)                                                                                                                                  aaCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTTGTCTTGCTGTa                                                       >  2:100532/1‑112 (MQ=255)                                                                                                                                   tcgcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgtgt                <  2:284536/148‑1 (MQ=255)                                                                                                                                   aCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgtgt                <  6:49993/149‑1 (MQ=255)                                                                                                                                    cgcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCCGGCCTTGGtgtc               <  7:119063/149‑1 (MQ=255)                                                                                                                                     gcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtc               >  3:69146/1‑148 (MQ=255)                                                                                                                                     attgcaGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtc               >  2:232416/6‑148 (MQ=255)                                                                                                                                        gcaGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcga             <  1:232416/146‑1 (MQ=255)                                                                                                                                        cgatttcggcCGAGGGCGCCCTGATATGTTCAGGTCAGAACTTGTAGGAAACACCGACCATGTTTTCCCCCGGGTCGAGGTGGTCGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatg           <  3:66374/143‑1 (MQ=255)                                                                                                                                          aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                    <  6:171487/138‑1 (MQ=255)                                                                                                                                          aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                    >  5:171487/1‑138 (MQ=255)                                                                                                                                          aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgtc         >  2:272096/1‑149 (MQ=255)                                                                                                                                          aGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgtc         <  2:287972/149‑1 (MQ=255)                                                                                                                                             gcggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGAc      <  2:183862/149‑1 (MQ=255)                                                                                                                                               ggcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTGTTGCTGTAGCC‑AAGGCGGGTCGGGCCGTTGCCCCTGGCCTTGGtgtcgatgtagatcg    >  4:66374/1‑146 (MQ=255)                                                                                                                                                   gAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAc  >  6:195205/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGC‑ACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCC‑AAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAC  >  NC_002947/593288‑593579

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