Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	593,434	A→C	6.1%	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.1%	242.4 / 6.9	82	E391A (GAG→GCG) ‡	PP_0503	MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (41/36);  new base C (2/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

ATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCA  >  NC_002947/593290‑593581                                                                                                                                                  |                                                                                                                                                     aTTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgcc                                                                                                                                                  >  1:178938/1‑149 (MQ=255) aTTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgcc                                                                                                                                                  >  7:285038/1‑149 (MQ=255)  ttCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgccg                                                                                                                                                 >  1:300617/1‑149 (MQ=255)     gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                              >  5:248662/1‑149 (MQ=255)     gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                              >  6:261335/1‑149 (MQ=255)     gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                              <  3:333534/149‑1 (MQ=255)     gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                              >  3:188008/1‑149 (MQ=255)     gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                                                                              >  2:389406/1‑149 (MQ=255)       tGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTg                                                                                                                                            <  1:169890/149‑1 (MQ=255)         ggCAGCGGTTTGTGGGCGGTGGTGGCGAGCTGTGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCtcgc                                                                                                                                          <  5:377536/149‑4 (MQ=255)           cAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAAt                                                                                                                                        <  3:116255/149‑1 (MQ=255)            aGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATg                                                                                                                                       <  2:178938/149‑1 (MQ=255)               ggTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGtt                                                                                                                                     >  4:331220/1‑148 (MQ=255)               ggTGTGGGGGCGGTGGTGGCGAGCGGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTc                                                                                                                                    <  7:282081/149‑1 (MQ=255)                 tttGGGGGCGGTGTTAGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                                                                  >  4:317485/1‑149 (MQ=255)                   tGGGGGGGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCGGGTGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCTAGGGCGGCCGAGGCCGCCCTGCAATGTTCAgg                                                                                                                                 <  1:328778/148‑1 (MQ=255)                    gggggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGt                                                                                                                                <  6:248662/148‑1 (MQ=255)                     ggggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGTGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑c                                                                                                                              <  5:369225/148‑1 (MQ=255)                     ggggCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑c                                                                                                                              <  2:300617/148‑1 (MQ=255)                        gCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGaa                                                                                                                          <  8:243362/149‑1 (MQ=255)                            tgtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTg                                                                                                                      <  1:298606/149‑1 (MQ=255)                            tgtTAGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTg                                                                                                                      <  3:317485/149‑1 (MQ=255)                              ttGGCGAGCTGGGTGTTTGCACACTGGGTCGTGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTa                                                                                                                    <  5:6618/149‑1 (MQ=255)                              ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGt                                                                                                                                <  5:359360/138‑1 (MQ=255)                              ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGt                                                                                                                                >  6:359360/1‑138 (MQ=255)                                 gCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGc                                                                                                                  >  3:283304/1‑148 (MQ=255)                                   gAGCTGGGTGTTTGCAGACTGGGGCTGGGTTGCGGTGTGTTGGATTGGCATGTGCATCAGTGCCGCCGCACTG‑TGCTACTGTCTATGCCTGAAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCaa                                                                                                               <  5:261335/149‑1 (MQ=255)                                    aGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAg                                                                                                              >  7:330773/1‑149 (MQ=255)                                         ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAgaccga                                                                                                         <  1:207263/149‑1 (MQ=255)                                             tttGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCAt                                                                                                     <  6:152611/149‑1 (MQ=255)                                                  acacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGCTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCt                                                                                                                                             <  1:297312/106‑1 (MQ=255)                                                  acacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCt                                                                                                                                             >  2:297312/1‑105 (MQ=255)                                                   cacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATAc                                                                                               >  4:393868/1‑149 (MQ=255)                                                    acTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATAcc                                                                                              >  7:87619/1‑149 (MQ=255)                                                                 tGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTAGTAGCCAAGACCGACCATGTATACCCACGGgtcgacg                                                                                  >  8:158775/1‑148 (MQ=255)                                                                    ggTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcga                                                                              >  6:201132/1‑149 (MQ=255)                                                                    ggTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcga                                                                              >  7:318537/1‑149 (MQ=255)                                                                        tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACAGACCATGTATACCCACGGgtcgacgtcgacgtc                                                                          >  8:313686/1‑149 (MQ=255)                                                                          tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcga                                                                        <  4:188008/149‑1 (MQ=255)                                                                          tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcga                                                                        >  4:90535/1‑149 (MQ=255)                                                                           gggCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcgac                                                                       <  2:405002/149‑1 (MQ=255)                                                                           gggCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcgac                                                                       >  4:28/1‑149 (MQ=255)                                                                            ggCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAcc                                                                      <  5:201132/149‑1 (MQ=255)                                                                             gCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAccc                                                                     <  3:393868/149‑2 (MQ=255)                                                                                   cATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCt                                                               >  3:145816/1‑149 (MQ=255)                                                                                    aTGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAcccgtgg                                                                 <  2:351604/146‑6 (MQ=255)                                                                                       ggCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                             >  1:98608/1‑131 (MQ=255)                                                                                       ggCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                             <  2:98608/131‑1 (MQ=255)                                                                                          aTCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTATTGCTg                                                          >  6:64549/1‑147 (MQ=255)                                                                                          aTCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGCCGGCCGAGGCCGCCCTGCAATGTTAAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCAACGTCGACGTCGACCTTGGTCTTGCTGTa                                                        <  7:169640/149‑1 (MQ=255)                                                                                           tCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg                                                       <  6:242725/149‑1 (MQ=255)                                                                                                  cgccgcACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAggg                                                <  8:330773/149‑1 (MQ=255)                                                                                                  cgccgcACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTCCAGGT‑CAGAACTCGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAggg                                                >  7:212007/1‑149 (MQ=255)                                                                                                        aCTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCg                                          >  8:212574/1‑149 (MQ=255)                                                                                                                tACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      >  4:236308/1‑98 (MQ=255)                                                                                                                tACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      <  3:236308/98‑1 (MQ=255)                                                                                                                   tGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGAc                                >  8:135529/1‑149 (MQ=255)                                                                                                                       tATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTg                            <  3:331220/149‑1 (MQ=255)                                                                                                                       tATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGACGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTg                            >  8:301872/1‑149 (MQ=255)                                                                                                                         tGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGCCGCTGGc                          >  3:52888/1‑149 (MQ=255)                                                                                                                            cTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCGCGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCtt                       >  4:102562/1‑149 (MQ=255)                                                                                                                             tGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTc                                           >  6:151000/1‑128 (MQ=255)                                                                                                                             tGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                      >  5:25712/1‑149 (MQ=255)                                                                                                                               cAACGCGCGAGTGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTgg                     <  8:87619/148‑1 (MQ=255)                                                                                                                               cAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCGCGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                      <  3:102562/147‑1 (MQ=255)                                                                                                                               cAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      <  7:218138/83‑1 (MQ=255)                                                                                                                               aaaCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                      >  8:218138/2‑83 (MQ=255)                                                                                                                                   gcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtc                 >  6:364419/1‑148 (MQ=255)                                                                                                                                    cgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTGCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcg                >  1:279941/1‑149 (MQ=255)                                                                                                                                     gcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgat              <  7:212574/149‑1 (MQ=255)                                                                                                                                       gAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgt            <  5:364419/149‑1 (MQ=255)                                                                                                                                       gAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgt            >  6:169321/1‑149 (MQ=255)                                                                                                                                         gggCGGCCTCGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgat              >  8:117972/1‑145 (MQ=255)                                                                                                                                         gggCGGCCTCGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgat              <  7:117972/145‑1 (MQ=255)                                                                                                                                           gcggcCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGAc        >  5:351250/1‑149 (MQ=255)                                                                                                                                             ggcCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGTTCTTTCTGTAGCCAAGTGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgtggaag       >  4:34188/1‑146 (MQ=255)                                                                                                                                               ccGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCGTGGTGTCGATGTCGACGTAc    >  1:286201/1‑149 (MQ=255)                                                                                                                                               ccGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGGCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAc    >  7:97303/1‑149 (MQ=255)                                                                                                                                                cGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAcc   <  7:135529/149‑1 (MQ=255)                                                                                                                                                 gAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCg                                              <  3:241248/105‑1 (MQ=255)                                                                                                                                                 gAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCg                                              >  4:241248/1‑105 (MQ=255)                                                                                                                                                 gAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCa  <  4:145816/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     ATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTG‑TGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGT‑CAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCA  >  NC_002947/593290‑593581

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