Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,580,628	A→G	7.1%	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	7.1%	196.8 / 4.8	70	I255V (ATC→GTC)	benR	BenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base G (2/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

GCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGG  >  NC_002947/3580482‑3580774                                                                                                                                                    |                                                                                                                                                   gCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATc                                                                                                                                                  <  6:302110/149‑1 (MQ=255)   cTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCAc                                                                                                                                                >  8:333891/1‑149 (MQ=255)   cTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCAc                                                                                                                                                >  8:269989/1‑149 (MQ=255)     gAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGc                                                                                                                                              >  5:111549/1‑149 (MQ=255)            ccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCa                                                                                                                                       >  5:75736/1‑149 (MQ=255)            ccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCa                                                                                                                                       >  6:153445/1‑149 (MQ=255)             cGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCGTGATGCCCAAgtgc                                                                                                                                      <  1:180942/149‑4 (MQ=255)             cGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCGTGATGCCCAAgtgc                                                                                                                                      <  2:301401/149‑4 (MQ=255)               cAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTa                                                                                                                                    <  6:111549/149‑1 (MQ=255)                       gCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGc                                                                                                                                        <  7:22529/137‑1 (MQ=255)                       gCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGc                                                                                                                                        >  8:22529/1‑137 (MQ=255)                        cTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAg                                                                                                                                         <  8:183554/135‑1 (MQ=255)                        cTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAg                                                                                                                                         >  7:183554/1‑135 (MQ=255)                               gcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                    <  4:70778/149‑1 (MQ=255)                                cATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAAc                                                                                                                   <  4:61097/149‑1 (MQ=255)                                 aTTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACt                                                                                                                  <  3:311729/149‑1 (MQ=255)                                 aTTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACt                                                                                                                  <  5:153445/149‑1 (MQ=255)                                  ttCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                 >  6:117954/1‑149 (MQ=255)                                  ttCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                 >  2:252349/1‑149 (MQ=255)                                           aCATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGt                                                                                                         <  4:287328/148‑1 (MQ=255)                                             aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCa                                                                                                      >  1:101734/1‑149 (MQ=255)                                                      aaCCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                             <  6:75736/149‑1 (MQ=255)                                                        ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGa                                                                                           <  5:165791/149‑1 (MQ=255)                                                          tGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCGTCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGc                                                                                         >  1:73164/1‑149 (MQ=255)                                                              aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAATCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGAGAca                                                                                     >  1:274840/1‑148 (MQ=255)                                                              aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                     >  6:248906/1‑149 (MQ=255)                                                              aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                     >  6:365107/1‑149 (MQ=255)                                                              aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACCTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                     <  1:252349/149‑1 (MQ=255)                                                                  gAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCATGAGCGACCCTAc                                                                                 >  7:59104/1‑149 (MQ=255)                                                                   aaTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCAc                                                                                                                                                <  6:355694/85‑1 (MQ=255)                                                                   aaTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCAc                                                                                                                                                >  5:355694/1‑85 (MQ=255)                                                                       gTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTtgt                                                                             <  5:117954/148‑1 (MQ=255)                                                                          ggCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGaa                                                                                                              <  8:94935/112‑1 (MQ=255)                                                                          ggCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGaa                                                                                                              >  7:94935/1‑112 (MQ=255)                                                                            cTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGc                                                                                                       >  5:205611/1‑117 (MQ=255)                                                                            cTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGc                                                                                                       <  6:205611/117‑1 (MQ=255)                                                                                    cTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACg                                                                                                                                               >  4:188150/1‑69 (MQ=255)                                                                                    cTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACg                                                                                                                                               <  3:188150/69‑1 (MQ=255)                                                                                    cTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAgtgt                                                               >  4:12730/1‑149 (MQ=255)                                                                                       gCAGAGCAGGCTTGCATGAGTTTGCGTTCGGTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGa                                                             <  8:59104/148‑1 (MQ=255)                                                                                                gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                               >  8:203587/1‑73 (MQ=255)                                                                                                gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                               <  7:203587/73‑1 (MQ=255)                                                                                                    gCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTg                                               >  8:38021/1‑149 (MQ=255)                                                                                                     cATGAGGTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTgg                                              <  3:202186/149‑1 (MQ=255)                                                                                                       tGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGAt                                            >  7:41012/1‑149 (MQ=255)                                                                                                           tttGCGTTCGCTGTATGCGTTGTTCGACCAGCCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTACCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTAc                                         >  1:224065/1‑149 (MQ=255)                                                                                                           tttGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACg                                        <  7:38021/149‑1 (MQ=255)                                                                                                              gCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgc                                                                                                   <  6:26762/87‑1 (MQ=255)                                                                                                              gCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgc                                                                                                   >  5:26762/1‑87 (MQ=255)                                                                                                                 ttCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCGTCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTCGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTcc                                  <  2:73164/149‑1 (MQ=255)                                                                                                                 ttCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAt                                                       <  2:296732/128‑1 (MQ=255)                                                                                                                 ttCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAt                                                       >  1:296732/1‑128 (MQ=255)                                                                                                                 ttCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTcc                                  <  3:12730/149‑1 (MQ=255)                                                                                                                 ttCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTcc                                  <  7:269989/149‑1 (MQ=255)                                                                                                                 ttCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTcc                                  <  7:333891/149‑1 (MQ=255)                                                                                                                    gCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCtt                                >  1:220942/1‑148 (MQ=255)                                                                                                                                 ttCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGt                   <  5:365107/148‑1 (MQ=255)                                                                                                                                 ttCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCACTGGATTACGGTTTCCTTCATCTGGGGCGGtt                  >  8:242700/1‑149 (MQ=255)                                                                                                                                     aCCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc              >  3:35435/1‑149 (MQ=255)                                                                                                                                     aCCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTACCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc              <  2:224065/149‑1 (MQ=255)                                                                                                                                      ccAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACt                                                                                >  5:39463/1‑82 (MQ=255)                                                                                                                                      ccAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACt                                                                                <  6:39463/82‑1 (MQ=255)                                                                                                                                      ccAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCg             <  2:101734/149‑1 (MQ=255)                                                                                                                                      ccAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCg             <  5:248906/149‑1 (MQ=255)                                                                                                                                         g‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga          <  8:41012/149‑1 (MQ=255)                                                                                                                                             cTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGAAGTGCGACGGAATTGGCGCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGAtt        >  4:71545/1‑148 (MQ=255)                                                                                                                                             cTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgggg                       >  8:88367/1‑133 (MQ=255)                                                                                                                                              ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTa      >  5:118245/1‑149 (MQ=255)                                                                                                                                                 ggCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCg   >  3:160204/1‑149 (MQ=255)                                                                                                                                                   cATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCgg  >  3:21694/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAG‑CACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGG  >  NC_002947/3580482‑3580774

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