Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,580,632	C→T	7.9%	T256M (ACG→ATG)	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,632	0	C	T	7.9%	180.2 / 4.5	64	T256M (ACG→ATG)	benR	BenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (35/23);  new base T (2/3);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

CCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAG  >  NC_002947/3580491‑3580780                                                                                                                                              |                                                                                                                                                     ccccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                               <  1:234311/149‑1 (MQ=255) ccccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                               >  8:23839/1‑149 (MQ=255)  cccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                               >  6:277431/1‑148 (MQ=255)        ggCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGCCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgtgctggt                                                                                                                                        <  2:112214/149‑8 (MQ=255)           tGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                                     >  7:273573/1‑149 (MQ=255)             ggCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgcgc                                                                                                                                   <  3:99104/149‑1 (MQ=255)                         ttCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCAAgtgctg                                                                                                                                          >  4:141386/1‑125 (MQ=255)                         ttCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCAAgtgctg                                                                                                                                          <  3:141386/130‑6 (MQ=255)                          tCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                       >  4:308926/1‑148 (MQ=255)                              gATTACATCGATAGCTACCTGATATTGGAATTGTCTGCTGAAGTGGTTGCAGTGCAGGCTTGCATGAGTTTGGGTTGGCTGTATGCGTTGTTCGACCTGGACTTTGTCATCAGGCCCATGCTTTACTTGCGCCAGCGGAAACTGTAAc                                                                                                                   <  4:174923/148‑1 (MQ=255)                              gATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAAc                                                                                                                   <  8:273573/148‑1 (MQ=255)                               aTTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACgg                                                                                                                 <  7:23839/149‑1 (MQ=255)                                 tACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                          <  8:207150/138‑1 (MQ=255)                                 tACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                          >  7:207150/1‑138 (MQ=255)                                    aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCa                                                                                                            >  2:3407/1‑149 (MQ=255)                                    aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAACAGTACGTGCGCCAGCGGAAACTGGACCGGGTGCa                                                                                                            >  8:249297/1‑149 (MQ=255)                                       gAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgc                                                                                                         >  6:64901/1‑149 (MQ=255)                                         aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                       >  4:10214/1‑133 (MQ=255)                                         aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                       <  3:10214/133‑1 (MQ=255)                                         aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcct                                                                                                       >  1:243297/1‑149 (MQ=255)                                          cGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctg                                                                                                      >  1:68303/1‑149 (MQ=255)                                               ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCAAGAGTTTGCGTTCGCTGTATGCGATGTTCGACCAGCACTTGGGCAACACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGCTGAATgcctgcct                                                                                                   >  8:30601/1‑147 (MQ=255)                                                cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                                >  2:39025/1‑149 (MQ=255)                                                          aaTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGTGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACt                                                                                      >  6:59592/1‑149 (MQ=255)                                                            ttGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGAGCATGCCTGCCTGAGCGACCCTACTtg                                                                                    >  4:150021/1‑149 (MQ=255)                                                                 ggCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGc                                                                                <  5:277431/148‑1 (MQ=255)                                                                 ggCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGc                                                                                <  3:308926/148‑1 (MQ=255)                                                                  gCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGt                                                                              <  2:243297/149‑1 (MQ=255)                                                                   cTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTg                                                                             >  3:257305/1‑149 (MQ=255)                                                                   cTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTg                                                                             <  2:221451/149‑1 (MQ=255)                                                                           cTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAgtgt                                                                     >  8:245058/1‑149 (MQ=255)                                                                              gCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                                  <  7:245058/149‑1 (MQ=255)                                                                              gCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                                  <  5:64901/149‑1 (MQ=255)                                                                                    cAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAtt                                                            >  6:103652/1‑149 (MQ=255)                                                                                     aGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCTGAAACTGGAACGGGTGCATGCCTGCATGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTg                                                           >  4:301969/1‑149 (MQ=255)                                                                                       gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                                   >  5:27530/1‑107 (MQ=255)                                                                                       gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                                   <  6:27530/107‑1 (MQ=255)                                                                                       gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGaa                                                              <  1:9642/144‑1 (MQ=255)                                                                                       gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGaa                                                              >  2:9642/1‑144 (MQ=255)                                                                                         ttGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGccc                                                       >  8:23362/1‑149 (MQ=255)                                                                                            cATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTg                                                     >  2:38328/1‑148 (MQ=255)                                                                                             aTGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                                   <  1:39025/149‑1 (MQ=255)                                                                                             aTGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                                   >  8:41363/1‑149 (MQ=255)                                                                                              tgagtttgagttCGCTGTATGCGTTGTTAGACCAGCACTTGGGCATCACGCACAAGCATTACGTGCTCCAGCGGAAACTGGAACGTGTGCATGCCTGCCTGAGCGACCCTACTTGTGGAGTGCGCAGTGTGACGGAATTGGCCCTGGAt                                                  >  6:218497/1‑149 (MQ=255)                                                                                              tGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGAt                                                  >  5:134285/1‑149 (MQ=255)                                                                                                    tGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGACAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGt                                            >  1:250026/1‑149 (MQ=255)                                                                                                        ttCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                                     >  7:104217/1‑88 (MQ=255)                                                                                                        ttCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                                     <  8:104217/88‑1 (MQ=255)                                                                                                           gCTGTATGTGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTc                                     <  5:59592/149‑1 (MQ=255)                                                                                                             tGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGAGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCAt                                   >  6:145104/1‑149 (MQ=255)                                                                                                               tATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCt                                 <  2:68303/149‑1 (MQ=255)                                                                                                                 tGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                               >  3:98999/1‑149 (MQ=255)                                                                                                                   cGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGAGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTggg                              <  3:150021/148‑1 (MQ=255)                                                                                                                    gttgttCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGc                            >  2:9032/1‑149 (MQ=255)                                                                                                                           gACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTc                     >  1:316680/1‑149 (MQ=255)                                                                                                                             ccAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                    >  2:15269/1‑148 (MQ=255)                                                                                                                              cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGACCGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                    >  3:289524/1‑147 (MQ=255)                                                                                                                              cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCACTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCga                  <  7:187512/149‑1 (MQ=255)                                                                                                                                gCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga                <  7:23362/149‑1 (MQ=255)                                                                                                                                 tgcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCGGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga                <  2:195520/146‑1 (MQ=255)                                                                                                                                 tgcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCGGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga                >  1:195520/3‑148 (MQ=255)                                                                                                                                       ggCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCg         <  3:191820/149‑1 (MQ=255)                                                                                                                                             aCGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTc                     >  3:178479/1‑131 (MQ=255)                                                                                                                                             aCGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTc                     <  4:178479/131‑1 (MQ=255)                                                                                                                                              cGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGcagcag  >  2:309554/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     CCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAGCAG  >  NC_002947/3580491‑3580780

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