Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,578,290	0	A	C	42.3%	88.1 / 12.6	53	V146G (GTG→GGG)	gntR	d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (22/8);  new base C (0/22);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 2.20e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCAT  >  NC_000913/3578141‑3578438                                                                                                                                                           |                                                                                                                                                       tgatAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑ccac                                                                                                                                                       >  2:812373/1‑151 (MQ=255) tgatAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑ccac                                                                                                                                                       >  2:700994/1‑151 (MQ=255)     aGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑Ccacca‑c                                                                                                                                                   <  1:876876/150‑1 (MQ=255)        aCGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑Cggg                                                                                                                                                <  1:218924/150‑1 (MQ=255)        aCGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑Cggg                                                                                                                                                <  2:641992/150‑1 (MQ=255)          gTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AAt                                                                                                                                            >  2:99904/1‑151 (MQ=255)          gTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AAt                                                                                                                                            <  1:316445/151‑1 (MQ=255)             cgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATAcc                                                                                                                                         <  2:1032336/151‑1 (MQ=255)              gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATAccc                                                                                                                                        >  1:4528/1‑151 (MQ=255)              gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATAccc                                                                                                                                        <  2:1114903/151‑1 (MQ=255)                 gagaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGcc                                                                                                                                     >  2:149913/1‑151 (MQ=255)                   gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCa                                                                                                                                    <  1:633439/150‑1 (MQ=255)                                   ggaattggggggATGCCCGCGAGAAATAAGGGCGGTGT‑TCATCTGGCGTC‑GTGTTTAAAATTTATAAAAC‑CCGACGGCGATATCAGGGCAGGGCGCCTTGCTGTCC‑ATCAGTC‑CCACCA‑GGGA‑AATCCCCGCCACTTCAATCATCttt                                                                                                                      <  1:559619/137‑1 (MQ=255)                                   ggCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTa                                                                                                                     >  1:485358/1‑149 (MQ=255)                                    gCAATGTGGCGATGCCCGCGAGCAAT‑ATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGt                                                                                                                 >  1:384749/1‑151 (MQ=255)                                     aaagggggggATCCCCGGGAAAAATAAGGGCGGTGG‑TCTTGGGGGGGG‑CGGCTACAAATTTATAAAAC‑CCGACGGCGATATCAAGGCATGGCTATTTGTTGTCC‑ATCGGTT‑CCCCCA‑CGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGt                                                                                                                 <  2:1057720/142‑1 (MQ=255)                                      aaTGGGGCGATGCCCGCAAAAAAAAAGGGCGGTGG‑TCACCGGGCGTG‑CTGCTCAAAATTTAAAAAAC‑CCGCCGGCAATAAAAAGGCAGGGCTCTTTGCTGCCC‑ATCAGTC‑CCCCCC‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                                  <  2:469651/149‑1 (MQ=255)                                        gggggCGATGCCGGCGAGAAATAGTGGAGGTGG‑TCATC‑GGGGTGCCTGCTAAAAAGTTACAAAAC‑CCGACGGCGATATCAGGCCGTGGCGTCTTGCTGTCC‑ATCATTC‑CCCCCA‑GGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                               <  2:404166/147‑1 (MQ=255)                                           ggggATGCCCGGGAGAAATAAGGGCGGTGG‑TCATCTGGCGTT‑TTGCTCAAAAGTTAAAAAAC‑CCGACGGCGATATCAAGGCAGGGCGCCTTCCTGTCC‑ACCAGTC‑CCCCCA‑CGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                               <  1:756390/144‑1 (MQ=255)                                                        aGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGt                                                                                              >  1:253293/1‑151 (MQ=255)                                                            ataatGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                          >  1:451009/1‑151 (MQ=255)                                                              aatGGCGGTGG‑TCTTCGGGCGTC‑CT‑CTTTAAAATTAACAAAACCCGCCGGCAATACAAAGGCAGGGCTCTTTGCTGCCC‑ATCAGTC‑CCCCCA‑CGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg                                                                                          <  2:708115/149‑1 (MQ=255)                                                               atGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                                        >  1:148000/1‑150 (MQ=255)                                                               agggCGGTGG‑TCATCTGGTTTG‑CTTATTTAAGGTAACAAACG‑CCCGCGGCGATAACAAGGCAGGGCGTCTTGCTGTCC‑ACCAGTT‑CCCCCA‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                         <  1:726612/147‑1 (MQ=255)                                                                  ggaGTGGTTCATCTGGGGTG‑CTG‑TTCAAAATTTACAAAACCCGCCGGCGATATAAAGGCGTGGCGCTTTGCTGTCC‑ATCAGTC‑CCCCCA‑GGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                      <  2:1007465/148‑1 (MQ=255)                                                                     gTGG‑TCTTCGGGCGTG‑CTCCTCAAAATTTAACAAAC‑CCGACGAGAATACCAGGGCATGGCGTCTTGTTCCCC‑ATCATTC‑CCCCCA‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                 <  2:607016/151‑1 (MQ=255)                                                                        gttcttCTGGCGTG‑CTGCTTAAAATTTATAAAAC‑CCGCCGGCAATATAAAGGCATGGCGCCTTGTTGCCC‑ACCAGTC‑CCCCCA‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                                                                                <  1:653334/145‑1 (MQ=255)                                                                        ggtcttggggCGTG‑CTGTTTAAAAGTTATAAAAC‑CCGAGGGCAATACCAAGGCAGGGCGACTTGCTGCCC‑ACCAGTC‑CCCCCG‑CGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                  <  1:429843/140‑1 (MQ=255)                                                                             tcgggggTG‑CTGTTTAAAATTTAAAAAAC‑CCGGGGGCGATAACAAGGCAGGGCGATTTGCTTCCC‑ATCAGTC‑CCCCCG‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGt                                                                           <  1:464592/144‑1 (MQ=255)                                                                                ggggTG‑CTGCTCAAATTTTAAAAAAC‑CCGGCGGCGAAACAAGGGCAGGGCGTTTTCCTGTCC‑ACCAGTC‑CCCCCG‑GGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                            <  1:13839/143‑1 (MQ=255)                                                                                ggCGTG‑CTTCTAAAAATTTACAAACC‑CCGAGGGCAATATCAGGGCATGGCGTCTTGCTGTCC‑ATCA‑TT‑CCCCCC‑CGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                  <  1:742458/140‑1 (MQ=255)                                                                                ggCGTG‑CTGCTTAAAATTAAAAAAAC‑CCGCCGGCGATAACAAGGCATGGCGCTTTTCTGTCC‑ATCGGTC‑CCCCCG‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                                          <  1:251141/148‑1 (MQ=255)                                                                                   gTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                    >  1:424612/1‑151 (MQ=255)                                                                                     g‑ctgctTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                  >  1:652089/1‑151 (MQ=255)                                                                                        tgctaaaaaTTTACCAAAC‑CCGACGGCGATATCAAGGCAGGGCGTCTTGCTGTCC‑ATCAGTC‑CCCCCA‑CGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                  <  1:277430/143‑1 (MQ=255)                                                                                         gcataaaattaaacAAAC‑CCGAAGGCGAAATCAGGGCAGGGCGGTTTGCTCACC‑‑TTAGTTCCCCCCA‑CGGA‑AAGCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCa                                                                 <  2:1109013/136‑1 (MQ=255)                                                                                           ttCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                               >  1:574857/1‑149 (MQ=255)                                                                                             cAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                             >  2:1080674/1‑149 (MQ=255)                                                                                              aaatttaaaaaac‑CCGAGGGCAATATCAAGGCTGGGCGCCTGGCTGTCC‑ACCA‑TT‑CCCCCC‑CGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGga                                                              <  1:398417/138‑1 (MQ=255)                                                                                              aaatttaaaaaac‑CCGACGGCGATATCAGGGCAGGGCGTTTTGCTGTCC‑ATCAGTC‑CCACCA‑GGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                                   <  1:50606/133‑1 (MQ=255)                                                                                                   tATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                      >  1:443155/1‑150 (MQ=255)                                                                                                    aacAAAA‑CCGCCGGCAATATCAAGGCAGGGCGATTTGC‑GTCCTATCAGTC‑CCGCGA‑GGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                      <  2:630797/147‑1 (MQ=255)                                                                                                              gACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                           >  2:809026/1‑151 (MQ=255)                                                                                                                 ggCGATATCAAGGCATGGCGACTTGCTGTCT‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                        >  2:570392/1‑151 (MQ=255)                                                                                                                 ggCGATATCAAGGCATGGCGACTTGCTGCCC‑ATCAGTC‑CCCCCG‑CGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                        <  1:993670/151‑1 (MQ=255)                                                                                                                  gCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                       >  2:159523/1‑151 (MQ=255)                                                                                                                    gATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                     >  2:100875/1‑151 (MQ=255)                                                                                                                        tCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTcc                                  >  2:122503/1‑150 (MQ=255)                                                                                                                         cAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                                >  1:17959/1‑151 (MQ=255)                                                                                                                          aaGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt                               >  1:323218/1‑151 (MQ=255)                                                                                                                            ggcggggCGTCTTGCTGCCC‑ATCAGTT‑CCCCCA‑GGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAttt                              <  1:548127/145‑1 (MQ=255)                                                                                                                                 ggCGTCTTCCTCTCC‑ATCATTC‑CCCCCA‑GGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                         <  2:452712/150‑1 (MQ=255)                                                                                                                                 ggCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt                        >  2:353421/1‑151 (MQ=255)                                                                                                                                  gcccttctctcccc‑ATCA‑TT‑CCCCCATCGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                         <  1:771552/137‑1 (MQ=255)                                                                                                                                   ccttttGCTGCCC‑ATCAGTC‑CCCCCC‑CGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGttt                       <  1:920442/146‑1 (MQ=255)                                                                                                                                            gTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt               >  2:1131524/1‑149 (MQ=255)                                                                                                                                                         ccacca‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCAt  >  1:64793/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                       TGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGG‑TCATCTGGCGTG‑CTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCC‑ATCAGTT‑CCACCA‑CGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCAT  >  NC_000913/3578141‑3578438

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