Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,586,352	G→A	5.9%	A182A (GCG→GCA)	benK →	benzoate MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,586,352	0	G	A	5.9%	202.3 / 2.0	68	A182A (GCG→GCA)	benK	benzoate MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base A (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

AACGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCGATGATGTGCTTC  >  NC_002947/3586209‑3586497                                                                                                                                                |                                                                                                                                                  aaCGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGtcc                                                                                                                                               >  4:260758/1‑146 (MQ=255) aaCGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGc                                                                                                                                               <  8:348448/148‑1 (MQ=255) aaCGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGc                                                                                                                                               <  3:292730/148‑1 (MQ=255) aaCGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGc                                                                                                                                               <  3:100378/148‑1 (MQ=255)  aCGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCAGTGGCa                                                                                                                                             >  6:112942/1‑149 (MQ=255)      aTACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCActgc                                                                                                                                         >  7:30209/1‑149 (MQ=255)          gCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCActgctgct                                                                                                                                     <  2:23092/149‑1 (MQ=255)           cACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCActgctgctg                                                                                                                                    >  8:277198/1‑149 (MQ=255)           cACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCAGTGGCAccgctgcgg                                                                                                                                    <  5:112942/149‑3 (MQ=255)              ccAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTa                                                                                                                                 >  4:243457/1‑149 (MQ=255)              ccAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTa                                                                                                                                 >  2:134629/1‑149 (MQ=255)              ccAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTa                                                                                                                                 >  3:26612/1‑149 (MQ=255)              ccAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCAGTGGCAccgctgcgg                                                                                                                                    <  4:34523/146‑3 (MQ=255)              ccAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCAGTGGCAccgctgcgg                                                                                                                                    >  3:34523/1‑144 (MQ=255)               cAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTAc                                                                                                                                >  6:55509/1‑149 (MQ=255)               cAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTAc                                                                                                                                <  3:243457/149‑1 (MQ=255)                aaGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTAcc                                                                                                                               <  7:198673/149‑1 (MQ=255)                    gcCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGt                                                                                                                            >  8:198687/1‑148 (MQ=255)                     ccTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGa                                                                                                                          <  5:291009/149‑1 (MQ=255)                       tGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGAtt                                                                                                                        >  1:181713/1‑149 (MQ=255)                         cgcAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATtct                                                                                                                      <  2:240443/149‑1 (MQ=255)                          gcagcaCGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATtctc                                                                                                                     <  7:277198/149‑1 (MQ=255)                          gcagcaCGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATtct                                                                                                                      >  1:128744/1‑148 (MQ=255)                           cagcaCGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTctct                                                                                                                    <  1:235708/149‑1 (MQ=255)                            agcaCGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTcta                                                                                                                   <  2:92552/149‑1 (MQ=255)                                              aTGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCAt                                                                                                 <  7:198687/149‑1 (MQ=255)                                              aTGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCAt                                                                                                 <  8:212302/149‑1 (MQ=255)                                              aTGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCAt                                                                                                 <  5:55509/149‑1 (MQ=255)                                              aTGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCAt                                                                                                 <  4:26612/149‑1 (MQ=255)                                               tGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTTGCTTGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATc                                                                                                <  1:284714/149‑1 (MQ=255)                                               tGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCTCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATc                                                                                                <  8:42306/149‑1 (MQ=255)                                                    aGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCActgctg                                                                                                                                       <  6:78518/105‑1 (MQ=255)                                                    aGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCActgctg                                                                                                                                       >  5:78518/1‑105 (MQ=255)                                                    aGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCtt                                                                                           >  3:149135/1‑149 (MQ=255)                                                     gTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGGGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTc                                                                                          >  3:227139/1‑149 (MQ=255)                                                     gTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTc                                                                                          >  8:109152/1‑149 (MQ=255)                                                      tGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTcc                                                                                         >  8:242275/1‑149 (MQ=255)                                                         ttATTCGCTGGGCGGCATGCTGTTGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTgg                                                                                      >  2:107498/1‑149 (MQ=255)                                                         ttATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTgg                                                                                      <  8:201541/149‑1 (MQ=255)                                                            ttCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCtcgata                                                                                   >  7:339921/1‑144 (MQ=255)                                                            ttCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTa                                                                                    <  4:149135/148‑1 (MQ=255)                                                                 tGGGCGGAATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGAAGTTCTTCGCAGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCaa                                                                              >  5:17549/1‑149 (MQ=255)                                                                         aTGCTGTCTGGAGGTGTCGGCATGTTCATGCTGCCGCGTTTTGGCTGCGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTAt                                                                      <  5:113137/149‑1 (MQ=255)                                                                            cTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGa                                                                   <  1:341390/149‑1 (MQ=255)                                                                              gTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAg                                                                  >  6:188328/1‑148 (MQ=255)                                                                                 ggCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGaa                                                               <  5:162422/148‑1 (MQ=255)                                                                                 ggCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAg                                                              >  4:178009/1‑149 (MQ=255)                                                                                    aGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAgcac                                                           <  2:128744/149‑1 (MQ=255)                                                                                         tCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGcaca                                                      <  8:30209/149‑1 (MQ=255)                                                                                             cATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGc                                                  >  6:49856/1‑149 (MQ=255)                                                                                                       ctgcCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCt                                        <  2:181713/149‑1 (MQ=255)                                                                                                       ctgcCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCt                                        >  4:218823/1‑149 (MQ=255)                                                                                                          ccGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGa                                     >  6:289250/1‑149 (MQ=255)                                                                                                           cgcgTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGAc                                    >  5:157471/1‑149 (MQ=255)                                                                                                            gcgTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGAcc                                   <  5:188328/149‑1 (MQ=255)                                                                                                                            tCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCa                   <  1:134629/149‑1 (MQ=255)                                                                                                                              gATTTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGCCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAgg                 <  6:66676/149‑1 (MQ=255)                                                                                                                              gATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAgg                 >  5:230828/1‑149 (MQ=255)                                                                                                                              gATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGACAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAg                  <  7:242275/148‑1 (MQ=255)                                                                                                                               aTGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGc                >  4:3299/1‑149 (MQ=255)                                                                                                                               aTGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGc                <  3:178009/149‑1 (MQ=255)                                                                                                                               aTGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGc                >  6:194550/1‑149 (MQ=255)                                                                                                                                tGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGcc               <  3:260758/149‑1 (MQ=255)                                                                                                                                  ttcttcGCCGCAGCGGTGCCACTGCCGCTGGAACCGGTGATTCTCTACTACCTCCCTGAATCCATCGGCTTCCTCGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCga             <  8:240337/149‑1 (MQ=255)                                                                                                                                        gccgcAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCGATGAtgt       >  4:175113/1‑149 (MQ=255)                                                                                                                                           gcAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCgat            <  5:268548/141‑1 (MQ=255)                                                                                                                                           gcAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCgat            >  6:268548/1‑141 (MQ=255)                                                                                                                                             aGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCGATGATGTGCTTc  <  5:289250/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  AACGAATACGCACCCAAGCGCCTGCGCAGCACGCTGGTGGCGATCATGTTCAGTGGTTATTCGCTGGGCGGCATGCTGTCGGCAGGTGTCGGCATCTTCATGCTGCCGCGTTTTGGCTGGGAGTCGATGTTCTTCGCCGCAGCGGTGCCACTGCTGCTGTTACCGGTGATTCTCTACTACCTGCCTGAATCCATCGGCTTCCTGGTACGCCAAGGCCGTATCGAGGAAGCACGCACATTGCTCAAACGGCTGGACCCGGACTGCGATGTCCAGGCCGATGATGTGCTTC  >  NC_002947/3586209‑3586497

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