Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,578,287	0	T	C	41.0%	105.8 / 11.7	61	E147G (GAA→GGA)	gntR	d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (31/5);  new base C (0/25);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGC  >  NC_000913/3578140‑3578436                                                                                                                                                       |                                                                                                                                                      atgatAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTcc                                                                                                                                                     >  1:195862/1‑150 (MQ=255)     tAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCcaccac                                                                                                                                                >  1:403859/1‑151 (MQ=255)      aGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACg                                                                                                                                               >  1:22622/1‑151 (MQ=255)        tACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACggg                                                                                                                                             >  2:953834/1‑151 (MQ=255)          cGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGaa                                                                                                                                           <  2:1057364/151‑1 (MQ=255)           gTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAAt                                                                                                                                          <  2:436409/151‑1 (MQ=255)            ttCGTCGAGACGTGCGCCGCGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAAt                                                                                                                                          >  1:884475/1‑150 (MQ=255)              cgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATAcc                                                                                                                                       >  2:946/1‑151 (MQ=255)                tcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCg                                                                                                                                     >  1:458058/1‑151 (MQ=255)                       gTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTc                                                                                                                              >  1:382915/1‑151 (MQ=255)                       gTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTc                                                                                                                              >  1:1123100/1‑151 (MQ=255)                                 aTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCttt                                                                                                                    >  1:105845/1‑151 (MQ=255)                                    ggataTGGGGCGTTCCCGGCGAAAAATAAGGGTAGTGGTCATCGGGCGTGCTCCTTCAAATTTAAAAAAC‑CC‑GCGGGGAATATCAAGGCCGGGCTTCTTGC‑TGTCCATCATTCCCACCGCGGGAATCCCCGCCACTTCAATCATCTTTaa                                                                                                                  <  2:440649/146‑1 (MQ=255)                                    ggCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTa                                                                                                                   >  1:543830/1‑149 (MQ=255)                                      cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                                >  2:432824/1‑150 (MQ=255)                                       agtgtgGCGTCGCCGGCGAAAAAAAAGGGCAGGGGTCTTCGGCCGTCCTCCTCAAAAGTTAAAAAAC‑CC‑GCCGGCAATACCAAGCCAGGGC‑ACTTTCTTGCCCATCATTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAg                                                                                                                 <  2:616665/146‑1 (MQ=255)                                         tgtgGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTgcg                                                                                                            >  2:59142/1‑151 (MQ=255)                                          ggggCGATGCCCGCGAAAAAAAGGGGAAGTGGTCATCGGGCGTGCTGCTAAAAATTTAAAAAAA‑CC‑GGCGGCGATAACAAGGCATGGCGATTTGC‑TCCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCttt                                                                                                                    <  1:309207/140‑1 (MQ=255)                                               gATCCCCGCGAGAAATAAGGGCAGTGGTCATCGGGCGTGCTGCTTAAAAGTAAAAAAAC‑CC‑GCCGGCAATATAAAGGCAGGGCGACTTGC‑TGCCCATCAGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgt                                                                                                      <  1:832057/151‑1 (MQ=255)                                                     cgcgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                                                                >  1:53466/1‑151 (MQ=255)                                                         aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGt                                                                                            >  1:786101/1‑151 (MQ=255)                                                            aatagtGGCGGTGTTCATCTGGGGTGCTGCTTAAAATTTATAAAAC‑CC‑GACGGCAATATCAAGGCAGGGCGACTTGT‑TGCCCATCATTCCCCCCGGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                          <  1:991707/150‑1 (MQ=255)                                                            aataagGGTGGTGTTCTTGTGGCGTGCTCCTT‑AAAATTAACAAAACCC‑GCCGGCGATACCAGGGCATGGCGTTTTTC‑TTTCCACCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGt                                                                                            <  2:768358/148‑1 (MQ=255)                                                            aataagGGGTGTGGTCATCTGGCGTGCTTCTTAAATTTTAACAAAC‑CC‑GACGGCGATACCAAGGCATGGCTACTTGC‑TTTCCATCACTTCCCCCGCGGAAAGCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgt                                                                                                      <  2:1057050/138‑1 (MQ=255)                                                             ataatGGCAGTGGTCATCTGGCGTCCTGCTTAAAATTTATAAAAC‑CC‑GGGGGCGATATCAAGGCATGGCGCTTTGC‑TGTCCTTCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                        <  1:969398/151‑1 (MQ=255)                                                               aagggCGGTGGTCATCTGGCGTGCTGCTTAAAATTTAAAAAAC‑CC‑GCCGGCGATATAAGGGCATGGCGTCTTGC‑TGTCCATCAGTCCCCCCGGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                       <  1:106091/147‑1 (MQ=255)                                                                  ggCAGTGGTCATCTGGCGTGCTGCTTAAAAGTTACCAAAC‑CC‑GACGGCGATATCAAGGCGTGGCGACTTGC‑TGTCCATCAGTTCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGa                                                                                     <  2:721135/149‑1 (MQ=255)                                                                   gCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAgg                                                                                   >  1:996067/1‑150 (MQ=255)                                                                    cAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGAt                                                                                 >  1:12001/1‑151 (MQ=255)                                                                      gTGGTCTTCGGGGGTGCTGCTCAAAATTTACAAAAC‑CC‑GCGGGCAATAAAAAGGCAGGGCGCTTTGT‑TGCCCATCATTCCCCCCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGAt                                                                                 <  1:187265/149‑1 (MQ=255)                                                                            aTCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                          >  1:899340/1‑150 (MQ=255)                                                                                        gctTAAAATTTAAAAAAC‑CC‑GACGGCAATATCAAGGCTTGGCTCCTTGC‑TGTCCACCAGTCCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                              <  2:588595/150‑1 (MQ=255)                                                                                        cctcaaaaGTTACAAAAC‑CC‑CCCGGCGATATCAAGGCAGGGCGCTTTGT‑TCTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCa                                                               <  2:276051/144‑1 (MQ=255)                                                                                            aaaaGTTAACAAAC‑CC‑GACGGAAATATCAAGGCATGGCGATTTGC‑TGTCCATCGGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                           <  2:813105/148‑1 (MQ=255)                                                                                             aatattATAAAAC‑CC‑GCCGGCAATATCAAGGCATGGCGACTTGC‑TGTCCATCTGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                  <  2:960777/137‑1 (MQ=255)                                                                                             aaattaataaaac‑CC‑GACGGCGATATCAGGGCGGGGCGCCTGGC‑TGTCCATCGGTCCCACCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgaga                                                          <  2:46358/140‑1 (MQ=255)                                                                                             aaaTTTATCAAAC‑CC‑GCCGGCGATATCAAGGCATGGCGCCTTGC‑TGTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                             <  1:112530/146‑1 (MQ=255)                                                                                              aaatttaaaaacccc‑GAGGCCGATATAAAGGTAGGGCTTCTTGT‑TTCCCATCGGTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                             <  1:598724/133‑1 (MQ=255)                                                                                              aaatttaaaaaa‑CCGGACGGCGATATCAAGGCATGGCGCTTTCC‑TGTCCATCAGTCCCCCCCGGGAAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgaga                                                          <  1:154775/141‑1 (MQ=255)                                                                                              aaGTTATCAAAC‑CC‑GGGGGCGATAGCACGGAAGGGCGACTTTC‑TTTCCATCATTTCCCCCAGGGAAATTCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                                        <  2:869559/134‑1 (MQ=255)                                                                                               aGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCAt                                                      >  2:555320/1‑151 (MQ=255)                                                                                                    tCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                 >  1:929025/1‑151 (MQ=255)                                                                                                      aaac‑CC‑GACGGAGATACGAAGGCAGGGCGACTTGT‑TTTCCATCGTTCCCACCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTc                                                <  2:917214/150‑1 (MQ=255)                                                                                                      aaaa‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCg                                               >  2:421908/1‑151 (MQ=255)                                                                                                        aacccccacgacgATAACAAGGCATGCCCTCTTGC‑TGCCCATCATTCCCCCCGGGAGAATCACCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                 <  1:882779/142‑1 (MQ=255)                                                                                                        aaccc‑gccgGCGATATCAGGGCGTGGCGTCTTGC‑TGTCCATCATTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGa                                              <  1:944654/148‑1 (MQ=255)                                                                                                            c‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCg                                          >  1:893390/1‑151 (MQ=255)                                                                                                                    gATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                   >  2:84553/1‑151 (MQ=255)                                                                                                                    gATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                   >  1:134576/1‑151 (MQ=255)                                                                                                                         cAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                              >  1:503688/1‑151 (MQ=255)                                                                                                                            ggcgtggcgtCTTGC‑TTTCCATCGTTCCCCCCGCGGAAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAttt                            <  2:47084/150‑1 (MQ=255)                                                                                                                             gCATGGCTTCTTGC‑TCTCCATCATTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTcc                          <  1:749887/151‑1 (MQ=255)                                                                                                                                 gccGACTTCC‑TGCCCATCAGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                       <  1:713541/148‑1 (MQ=255)                                                                                                                                     aCTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                   >  1:80597/1‑150 (MQ=255)                                                                                                                                      cTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtata                  >  1:212708/1‑150 (MQ=255)                                                                                                                                       ttGC‑TGTCCATCAGTTCCACCACGGGAATGCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGCTTATAAc                >  2:928414/1‑151 (MQ=255)                                                                                                                                       ttGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc                >  2:958901/1‑151 (MQ=255)                                                                                                                                             gTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt             >  1:636020/1‑149 (MQ=255)                                                                                                                                               ccATCAGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGt          <  1:113312/150‑1 (MQ=255)                                                                                                                                                 aTCATTCCCCCCACGGGAAACCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTg         <  2:884257/149‑1 (MQ=255)                                                                                                                                                  tCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgcgc      >  2:806864/1‑151 (MQ=255)                                                                                                                                                      ttCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGc  >  1:872538/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CC‑GACGGCGATATCAAGGCATGGCGACTTGC‑TGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGC  >  NC_000913/3578140‑3578436

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑