Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,580,628	A→G	6.2%	I255V (ATC→GTC) ‡	benR →	BenABC operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,580,628	0	A	G	6.2%	184.0 / 2.8	64	I255V (ATC→GTC) ‡	benR	BenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base G (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

AGCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTAT  >  NC_002947/3580481‑3580771                                                                                                                                                      |                                                                                                                                                agCCTGAGTGCCCCGCAGGCTGGGCTAGTGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCa                                                                                                                                                 <  7:133711/148‑1 (MQ=255) agCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGCCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCCTTGTT‑CGACCAGCACTTGGGAAt                                                                                                                                                >  4:130858/1‑149 (MQ=255)  gCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATc                                                                                                                                               >  8:1050/1‑149 (MQ=255)     tGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACg                                                                                                                                            >  1:286465/1‑149 (MQ=255)                 aGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGCCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTAc                                                                                                                                <  1:107595/149‑1 (MQ=255)                     tGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                            <  3:57956/149‑1 (MQ=255)                          tAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCTGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAg                                                                                                                       >  6:268690/1‑149 (MQ=255)                           agagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                      >  4:242432/1‑149 (MQ=255)                            gagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCTGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                      <  5:268690/148‑1 (MQ=255)                            gagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTTCCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                      <  1:236409/149‑1 (MQ=255)                            gagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                      <  6:52134/148‑1 (MQ=255)                              gcgcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGa                                                                                                                   <  4:91684/149‑1 (MQ=255)                                gcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                 <  5:126187/149‑1 (MQ=255)                                     cTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTgg                                                                                                             <  5:238054/148‑1 (MQ=255)                                      ttGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGaa                                                                                                           >  8:257368/1‑149 (MQ=255)                                          ttACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACgg                                                                                                        <  8:223843/148‑1 (MQ=255)                                          ttACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCAGAAACTGGAACggg                                                                                                       >  6:73890/1‑149 (MQ=255)                                          ttACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGATTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACggg                                                                                                       >  3:234906/1‑149 (MQ=255)                                             cATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGc                                                                                                    <  2:286465/149‑1 (MQ=255)                                              aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCa                                                                                                   <  3:242432/149‑1 (MQ=255)                                                cGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAt                                                                                                  >  1:35600/1‑148 (MQ=255)                                                cGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAt                                                                                                  <  2:35600/148‑1 (MQ=255)                                                   aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcc                                                                                               <  1:119541/148‑1 (MQ=255)                                                     gCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATTGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctg                                                                                             >  1:275652/1‑149 (MQ=255)                                                      cAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgcc                                                                                           >  4:141315/1‑149 (MQ=255)                                                       aaCCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                      >  6:315915/1‑105 (MQ=255)                                                         ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGa                                                                                        <  7:74499/149‑1 (MQ=255)                                                          cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                              <  6:116819/126‑1 (MQ=255)                                                          cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                              >  5:116819/1‑126 (MQ=255)                                                               aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                  <  7:1050/149‑1 (MQ=255)                                                                     aTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCGTGATGCCCAAgtgct                                                                                                                                  >  8:17956/1‑91 (MQ=255)                                                                     aTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCGTGATGCCCAAgtgct                                                                                                                                  <  7:17956/95‑5 (MQ=255)                                                                        gTCGGCTGAAGTGCTTGCATTGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGTCCAGCACTTGGGCATGACGCCCAAGCTTTACGTGCGCCAGCGGAAACTGGTACGGGTGCATGCCTGCCTGAGCGACCCTACTtgt                                                                          <  3:130858/148‑1 (MQ=255)                                                                         tCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                            <  2:151909/129‑1 (MQ=255)                                                                         tCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                            >  1:151909/1‑129 (MQ=255)                                                                            gCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGt                                                                     >  6:306812/1‑149 (MQ=255)                                                                            gCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGt                                                                     >  4:320112/1‑149 (MQ=255)                                                                                       ggCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGa                                                          >  7:234841/1‑149 (MQ=255)                                                                                          agagCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACg                                                        >  2:115984/1‑148 (MQ=255)                                                                                            agCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGaa                                                     >  5:274567/1‑149 (MQ=255)                                                                                             gCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAt                                                    <  3:320112/149‑1 (MQ=255)                                                                                                      cATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCAGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTgg                                           <  5:73890/149‑1 (MQ=255)                                                                                                            tttGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGTGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACg                                     >  6:117563/1‑149 (MQ=255)                                                                                                            tttGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGAGCAGTGTGACGGAATTGGCCCTGGATTACg                                     >  8:161554/1‑149 (MQ=255)                                                                                                                  ttCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACTGTTTcc                               <  4:234906/149‑1 (MQ=255)                                                                                                                    cGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCt                              <  3:141315/148‑1 (MQ=255)                                                                                                                            tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                      >  5:328644/1‑149 (MQ=255)                                                                                                                            tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                      >  2:27771/1‑149 (MQ=255)                                                                                                                            tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                      >  3:56128/1‑149 (MQ=255)                                                                                                                            tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                      >  1:154816/1‑149 (MQ=255)                                                                                                                               gttgtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTc                                >  1:181635/1‑136 (MQ=255)                                                                                                                               gttgtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTc                                <  2:181635/136‑1 (MQ=255)                                                                                                                                  gtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGt                >  1:332380/1‑149 (MQ=255)                                                                                                                                  gtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGt                >  5:112808/1‑149 (MQ=255)                                                                                                                                       gACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGtttt             <  7:257368/148‑1 (MQ=255)                                                                                                                                        aCCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc           >  6:293812/1‑149 (MQ=255)                                                                                                                                        aCCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc           >  7:264837/1‑149 (MQ=255)                                                                                                                                          cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCga         <  6:274567/149‑1 (MQ=255)                                                                                                                                            gCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag        >  8:302741/1‑148 (MQ=255)                                                                                                                                            gCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag        <  7:302741/148‑1 (MQ=255)                                                                                                                                              gcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTa                                       <  1:295133/114‑1 (MQ=255)                                                                                                                                              gcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTa                                       >  2:295133/2‑115 (MQ=255)                                                                                                                                                ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag        <  6:327306/144‑1 (MQ=255)                                                                                                                                                ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag        >  5:327306/1‑144 (MQ=255)                                                                                                                                                 tGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTACTTCATCTGGGGCGGTTTTCCGAGATTTAt  >  6:330946/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                AGCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTAT  >  NC_002947/3580481‑3580771

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