Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,468	0	C	G	43.3%	97.4 / 20.1	60	G330G (GGC→GGG)	ghxP	guanine/hypoxanthine permease, high affinity; guanine/hypoxanthine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base C (5/29);  new base G (26/0);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 2.48e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATC  >  NC_000913/4279323‑4279618                                                                                                                                                      |                                                                                                                                                            aGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGC‑GGC‑gg‑g‑gggc‑g                                                                                                                                                           >  2:91961/1‑147 (MQ=255)   gCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GT‑C‑GGGC‑GGT‑aaaa                                                                                                                                                    <  2:302798/151‑1 (MQ=255)       aCCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑Cgg                                                                                                                                               <  1:551399/151‑1 (MQ=255)         ctgctgGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAT‑GGG‑Gg‑c‑gggg‑ggggaaaaaggg                                                                                                                                                >  2:690460/1‑140 (MQ=255)             tgGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAG‑GGG‑GG‑G‑GGGG‑GGG‑GAAAA‑CAGttttgc                                                                                                                                         >  1:673025/1‑148 (MQ=255)             tgGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGCGGG‑G‑GGGG‑GGG‑AAAAC‑CGGtt                                                                                                                                             >  1:649980/1‑148 (MQ=255)                aTAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGG‑GG‑G‑GGGG‑GGG‑aaaaac                                                                                                                                                  >  1:966778/1‑138 (MQ=255)                 tAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑CGGGGG‑GGA‑AAAAC‑CGGttggcgg                                                                                                                                        >  1:608720/1‑144 (MQ=255)                       tGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGAAC‑GGG‑GG‑G‑CGGC‑GGT‑AAAAA‑CGCGTTGAcaggctttac                                                                                                                                >  1:846845/1‑141 (MQ=255)                          gCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTGC‑GGC‑GG‑G‑GGGG‑GGG‑AAAAA‑CGGGTTGGCGGCCATctacgat                                                                                                                            >  2:480849/1‑145 (MQ=255)                            aGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑CGG‑GG‑C‑GGGG‑GGG‑AAAAA‑CCGGTTGGCGGGCATCAcccgt                                                                                                                            >  1:804948/1‑146 (MQ=255)                              atcatcGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAA‑CGC‑GGCC‑GGGG‑GGA‑AAAAC‑AGGTTTCA‑GGCTATCAcc                                                                                                                               >  1:855399/1‑144 (MQ=255)                                     aCGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGtt                                                                                                                 <  1:31657/151‑1 (MQ=255)                                     aCGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGACGGGG‑GG‑G‑GGGG‑GGG‑AAAAC‑CGGttttcgggtctc                                                                                                                                   >  2:909845/1‑125 (MQ=255)                                       ggTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑C‑GGGG‑GGG‑AAAAA‑CCGGTTGACGGGTATCAACCTTgggggggtggt                                                                                                                 >  1:577873/1‑140 (MQ=255)                                        gTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑C‑GGGG‑GGG‑AAAAA‑CGGGTTGGCGGGCATCACCGGTGGGGGGCTGtttc                                                                                                               >  1:754983/1‑148 (MQ=255)                                         tGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGG‑GGG‑GAAAA‑CGGTTTTACGGCTATCACCCGTTGCGGGCTGtttc                                                                                                               >  1:527558/1‑147 (MQ=255)                                           ggAAAGCGCTGACCACTGACTCCATGAGCAGCGTTTTGTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GTGC‑GGT‑AAAAC‑CGGTTTGACGGGTATCACCGTTGGCGTGCTGTTcctcc‑t                                                                                                          <  2:236531/151‑1 (MQ=255)                                              aaGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑C‑GGGG‑GGG‑AAAAA‑CGGGTTTGCGGCTATCACCCCTTGCGTGCTGTTTCCCC‑TGAt                                                                                                       >  1:413215/1‑151 (MQ=255)                                               aGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑CGGGGG‑GGG‑AAAAC‑CGGTTT‑TCGGCGCTCAACGCTGGCGGGGTGTttctcccgctt                                                                                                       >  2:185994/1‑146 (MQ=255)                                                 cactgaccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑CGGGGG‑GGG‑AAAAC‑CGG‑TTGGGGGCGATCCCCCTTGGGGGGGTGTtcccccccatttt                                                                                                     >  1:64453/1‑144 (MQ=255)                                                  actgaccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTg                                                                                                   <  1:427382/151‑1 (MQ=255)                                                   ctgaccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGCGTCTGCGGCGGGGAA‑GGC‑GG‑C‑GGGG‑GGTGAAAAC‑CGGTTT‑ACGGGGTTCACCCTTGGCGGGGTGGTcctcc‑tcgtt                                                                                                      >  2:173116/1‑143 (MQ=255)                                                     gaccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCATCGTTGGCGTGCTGTTCCTCC‑TGATTCTGtt                                                                                                 <  2:258550/150‑1 (MQ=255)                                                       ccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCct                                                                                              <  1:612071/151‑1 (MQ=255)                                                       ccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCct                                                                                              <  2:909729/151‑1 (MQ=255)                                                       ccactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCct                                                                                              <  1:338685/151‑1 (MQ=255)                                                        cactgacTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑G‑GGGG‑GGG‑GAAAA‑CGGTTTGGCGGCCATCCCCGCTGGGGGGCCGTTCCTCC‑TGTTTCTGtt                                                                                                 >  2:1044984/1‑147 (MQ=255)                                                             acTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑G‑GGGG‑GGG‑GAAAC‑CGGGTTTGCGGCTATCACCCTTGGCGGGCTGTTCCTCC‑TTGTTTTTTTTCTCTCTcc                                                                                        >  2:483945/1‑151 (MQ=255)                                                               tCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGG‑GG‑CGGGGG‑GGG‑AAAAC‑CGGTTTTCCGGC‑GTTACCCGTGGGGGGGTGCTCCTCC‑TGAttttttttccctttcc                                                                                        >  2:753289/1‑135 (MQ=255)                                                                ccATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCg                                                                                       <  1:350685/149‑1 (MQ=255)                                                                  aTGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGctct                                                                                   <  1:625692/151‑1 (MQ=255)                                                                    gAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGctctct                                                                                 <  1:801861/151‑1 (MQ=255)                                                                    gAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGctctc                                                                                  <  2:48927/150‑1 (MQ=255)                                                                    gAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGctctc                                                                                  <  1:164500/150‑1 (MQ=255)                                                                      gcagcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGGC‑GGG‑GG‑G‑GGGC‑GGG‑GAAAA‑CGGGGTGGCGGCGCTTTCCGCTTGGGGGCTGTTCGTCC‑TTATTCTGTTTTTccctcccctcccct                                                                                >  2:609505/1‑146 (MQ=255)                                                                       cagcGTTTTCTCTGGCCTGGTTGGTGCGGCTCCGGCAGCAGTATACATCGAATCTGCGGCAGTTAC‑GGC‑GT‑C‑GTGC‑GGT‑AAAGC‑CTGTTTGACTGCTATCACCGTTGGCGTGCTTTTCCTCC‑GGAGTCTGTTCCTCTCTCCGctctc                                                                                  <  2:223607/147‑1 (MQ=255)                                                                        agcGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTAcc                                                                             <  2:949158/151‑1 (MQ=255)                                                                           gTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTGC‑GGC‑GG‑C‑GGGG‑GGG‑AAAAC‑CGGTTTGGCGGGTATAACCCCTGGCGTGCTGTTCCTCCTTGTTTTTGTT‑CTCTCTCCGCCCTCttt                                                                               >  2:289710/1‑145 (MQ=255)                                                                             tttCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGG‑GGT‑AAAAC‑CCGGTTTACGGCCATCACCCGTGGCGGGGGGCTCCTCC‑TTATTCTGTTCCTCTCCCCCccctccttcccccct                                                                        >  1:55781/1‑141 (MQ=255)                                                                                    ggCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGt                                                                <  1:1027850/151‑1 (MQ=255)                                                                                      ccTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTGC‑GGC‑GG‑C‑GGGCGGGG‑AAAAA‑CGG‑TTGACGGCCATCACCCGTTGCGGGCTGGTGCTCC‑TTGTTCTTGTCCTCTCTCCGCTCTCCTTCCTC‑CTTCCcggg                                                                 >  2:472321/1‑148 (MQ=255)                                                                                      ccTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑G‑GGGG‑GGG‑GAAACAGGGTTTGA‑GGCCTTCCCCCGTGGCGGGGTGTTCCTCC‑CTCGTTTTTTTCTCTCTCCCCCCTCTTTCCTCGGTTCCgggg                                                                 >  2:1056931/1‑149 (MQ=255)                                                                                         ggtgggtgCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCa                                                           <  2:90705/151‑1 (MQ=255)                                                                                         ggtgggtgCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGGAC‑GGC‑GG‑CGGGGG‑GGA‑AAAAC‑CGGTTT‑TCGGCTATCACCGGTGGGGGGCTGTTTCTCC‑TGGTTTTTTTTCTCTCTCCCCCCTCTTACCTC‑CTTCCcgggg                                                                >  1:318188/1‑145 (MQ=255)                                                                                                gCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑CGGGGG‑GGG‑AAAAC‑CGGTTTG‑CGGCTATTACCCTTTGGGGGCTGTTTCTCC‑CGATTTTTTTTCTCTCTCCCCCCTCTTTCCCC‑CTTCCCGGGGGTGCAACgg                                                       >  1:741086/1‑148 (MQ=255)                                                                                                  ggCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGG‑GGG‑GAAAC‑CCGGTTTTCGGGTATTCCCCTTGGGGGGGTGTTGCTCC‑CCGTTTTGTTTCTCTTTCCCCCCCCTTTCCTC‑GTTTCGGGGGGTGGAACGGGTcc                                                   >  2:347517/1‑150 (MQ=255)                                                                                                  ggCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGG‑GGG‑AAAAA‑CGGGTTTGCGGCTATCACCGTTGGCGGGCTGTTCCTCC‑TTATTTTTTTTCTCTCTCCCCTCTCTTTCCTC‑GTTCCCGGGGGTGCAACCGcccc                                                   >  1:654284/1‑150 (MQ=255)                                                                                                    cTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑G‑GGGG‑GGT‑AAAAA‑AGGGTTGACGGCTATTAACGGTGGCGGGGTGTTTCTCC‑CGATTTTTGTGCTCTTTCCCCCCTCTTTCCTC‑GTTCCGGGGGATGCAACGGGTcc                                                   >  1:892170/1‑148 (MQ=255)                                                                                                      ccGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGgcgc                                              <  2:966431/151‑1 (MQ=255)                                                                                                          cAGCGGTATACATCGAGTCTGCGGCGGGGGC‑GGG‑GG‑A‑GGGC‑GGG‑AAAAA‑CGGTTTTAACGCCTTCATCCTTGGAGGGCCGTTTCTTC‑TCATCATTTTCCTCTTTCCGCCCCCTTTTCTC‑TTTCCTGGGGATGCCACACCACCGgcgccgg                                           >  1:791974/1‑147 (MQ=255)                                                                                                            gCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATg                                         <  2:176046/150‑1 (MQ=255)                                                                                                              ggTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTCCCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTAc                                      <  1:15470/151‑1 (MQ=255)                                                                                                                       cGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCct                              <  2:754689/150‑1 (MQ=255)                                                                                                                         aGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCctgct                           <  1:853430/151‑1 (MQ=255)                                                                                                                          gTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCctgctg                          <  1:832734/151‑1 (MQ=255)                                                                                                                              gcggcgGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCtgctgatg                       <  2:148105/150‑1 (MQ=255)                                                                                                                                ggcggGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCtgctgatgct                     >  1:825076/1‑150 (MQ=255)                                                                                                                                 gcggGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCtgctgatgctga                   <  1:310769/151‑1 (MQ=255)                                                                                                                                     gTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCaa               <  2:989550/151‑1 (MQ=255)                                                                                                                                       aC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAAca             <  2:432695/151‑2 (MQ=255)                                                                                                                                       aC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAAc              <  2:1009498/150‑1 (MQ=255)                                                                                                                                              gg‑c‑gggc‑gGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCg        <  2:608449/151‑1 (MQ=255)                                                                                                                                                     gc‑ggT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAAt   <  1:753583/151‑1 (MQ=255)                                                                                                                                                      c‑ggT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATc  >  1:832293/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                            AGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAC‑GGC‑GG‑C‑GGGC‑GGT‑AAAAC‑CGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCC‑TGATTCTGTTCCTCTCTCCGCTCTCTTACCTC‑GTTCCGGGGTATGCAACGGCTCCGGCGCTGATGTACGTTGGCCTGCTGATGCTGAGCAACGTGGCGAAAATC  >  NC_000913/4279323‑4279618

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