Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,664	0	C	G	43.8%	19.0 / 22.2	49	pseudogene (940/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base C (16/11);  new base G (6/15);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862430‑3862902                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑a                                                                                                                                                                                                                                                <  1:1079146/239‑1 (MQ=255)      aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑TAAG‑GTGGGGG‑GATCa                                                                                                                                                                                                                                            >  1:1059983/1‑238 (MQ=255)      aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAG‑TGGCGGG‑GATCa                                                                                                                                                                                                                                            >  2:1079146/1‑238 (MQ=255)       tttGGGTCAGCATCATATCGGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑AATCAGc                                                                                                                                                                                                                                          <  1:1989095/239‑1 (MQ=255)                            cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAG‑GTGCGGG‑GATAAGC‑CAACCCGCCCCCTTTTccc                                                                                                                                                                                                                      >  1:1606364/1‑237 (MQ=255)                                           tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑GT‑TAAG‑TTGCGGG‑GAATAGCAAAGCCCCCCCTTTTTACCC‑AATTTTCggggcgggc                                                                                                                                                                                                     >  1:990989/1‑233 (MQ=255)                                                    aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑GT‑AAGG‑GTGCGGG‑GATCAGC‑AAGCCCGCCCCTTTTACCC‑GATTTTTCCGTGCGGGCGTG‑AAGTAc                                                                                                                                                                                          >  1:198988/1‑239 (MQ=255)                                                            gaaaacccatAGACAACCAGCCCCCCGCCCTTTTCCCCCACCACAA‑AAAAGAGATAAATACCGGTCGATAGAGAGCCCGCCAGCAGGGCGTGCACCAAAAAAAGAAAGGGAAAAAGAAACGGGAAGAAAAATCCCAGCGGTTC‑GATAATCCCCCCCAGCGGGGG‑GG‑TT‑AGGG‑TGGGGGA‑AACCAAA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTa                                                                                                                                                                                   <  2:2195119/229‑1 (MQ=255)                                                                                 ccccccgccATTTTGCCCCCCACCCAAAAAGATAAACATCACGGTCGACATGGAGCCCGCAAGCCCGGCGTGTCCCAAAAACAGAAAGGTAAAAATAAAAAGGAAGGTAAATTCCAGCGGTTC‑GTTAACCCCCCCCAGCAGGGC‑GG‑TT‑AGGG‑TGGCGGA‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCaa                                                                                                                                                             <  2:198988/232‑1 (MQ=255)                                                                                          aTGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGCGGG‑TT‑AGGG‑TGGCGGG‑GA‑CACC‑AAACCCCCCCCATTTACCC‑CCGTTTTCTGTGGCGGAGGG‑GAGGAAAACCCCAAAGAAAAGTCCACGGC‑GCCAAACAatttttt                                                                                                                                                   >  2:798189/1‑237 (MQ=255)                                                                                          aTGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGCGGG‑TT‑AAGG‑‑TGCGGG‑GATCAGC‑AAACCCGCCCCCTTTACCC‑CATTTTTCGGTGCGGCAGTG‑AAGTACa                                                                                                                                                                                         >  1:2124034/1‑202 (MQ=255)                                                                                                                aaaTACATTAGGGCCAAAATGGGGCCCGCCAGCGGGGCGTGCACCAAAAACACCAACGGTAAAATGAAGAGGAAGTTAAATCCCGGCGTTTG‑GGTAATCCCCCCCAGCAGGGG‑GG‑TA‑AGGG‑TGGGGGA‑AACCAAA‑AACCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTg                                                                                                                                                                                                 <  2:204924/172‑1 (MQ=255)                                                                                                                     cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGc                                                                                                                                                                                                                                                     <  1:2328574/120‑1 (MQ=255)                                                                                                                        tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGCGGG‑TA‑AAGG‑TTGGGGG‑AATAACC‑AAGCCCCCCCCTTTTTCCC‑GATTTTTCGTGGCGGCAG‑G‑GAGAAAAACCCCCAAAAAAAAACCACCGC‑ACCCAAAAATTTTTTGTTACCCTGCAGAGCAAAAACGCCCTccc                                                                                                                      >  2:834463/1‑238 (MQ=255)                                                                                                                        tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GGTTT‑AAGG‑TTG‑GGG‑GATCAGC‑AAACCCCCCACCTTTACCC‑CGTTTTTCGGGGCGGGAGGG‑aa                                                                                                                                                                                              >  1:1432316/1‑168 (MQ=255)                                                                                                                          cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGG‑GTGCGGG‑GAACAGC‑AAACCCCCCCCCTTTTCCC‑GATTTTTCGGGGCGGCAGGG‑AAggaaaacg                                                                                                                                                                                      >  1:204924/1‑167 (MQ=255)                                                                                                                          cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑AT‑AAAG‑GTGCGGG‑AAGCAGC‑AAAGCCGCCCCCTTTTCCC‑CCTTTTTCTGTGCGGGAGTG‑GAGGAAAACGCCCAAAAAAAGACCCCCGC‑ACCCAAGAATTTTTAGTTACCGTGCAGGGGAAAAACGccctaccgg                                                                                                                    >  2:878734/1‑232 (MQ=255)                                                                                                                               gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAAGG‑TGGGGGG‑aacaag                                                                                                                                                                                                                                           >  2:2328574/1‑116 (MQ=255)                                                                                                                               gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑GT‑TAGG‑GTGGGGG‑GAACACC‑CAAGCCGCCCCCTTTTCCC‑CACTTTTCGGTGCAGCCGGGAAAGAACAAC‑CCAAAAAAAAGACCCCCGC‑ACCAAAAAATTTTTTGTGTCCGTGCATGGCAAACACACCTTCCggggcaaa                                                                                                               >  2:1273341/1‑234 (MQ=255)                                                                                                                                 ctttGGGGCCGCCAGCGCGGCGTATACCAAAAACAGCAAGGGAGAAAAGAAGAGGAAGGTAAATTCCAGCGGTTC‑GGTAATCCCCCCCAGCGGGGG‑GG‑TT‑AGGG‑TGGGGGG‑AAACAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg                                                                                                             <  2:1175135/237‑1 (MQ=255)                                                                                                                                 catgggggCCGCCAGCGCGGCGTGCACCGAAAACAGCAAGGGTAAAAGGAAGAGGAAGGAAAATTCCGGGGGTTG‑GGTAATCCCCCCCAGCGGGGG‑GT‑TT‑AGGG‑TGGGGGG‑AAACAGC‑AACCCCGCCACTTTTACCC‑GATTTTCCGGTGCgg                                                                                                                                                                                                      <  2:455342/146‑1 (MQ=255)                                                                                                                                         ccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑GTGCGGG‑GATAAGCAAAGCCCCCCCCTTTTCCCC‑ATTTTTCCGG‑GCGGGGGGG‑GAGGAAAACGCCAAAAAAAAACCCCCCGCGACCAAA‑AATTTTGTGGTACCGGGCAGGGCAAAAACACCCTCCCGGGCAAAAGATTTCaac                                                                                                     >  2:327369/1‑237 (MQ=255)                                                                                                                                           gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAAG‑TTGCGGG‑GAACCGC‑AAACCCCCCCCCTTTTCCC‑GATTTTTCCGGGGGGCAGGG‑GAGTa                                                                                                                                                                                           >  1:455342/1‑152 (MQ=255)                                                                                                                                            ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGG‑GTGCGGG‑GAACAGC‑AAACCCCCCCCCTTTACCC‑CCTTTTTCGGGGGGGGCGGG‑GAGGACACCCCCAAAAAAAAAACCCCGGC‑ACCCAAAAATTTTTTGTTTACGTGCGTGGCCAAACACCCCTCCGGGgaaaacaagtttatagg                                                                                                   >  2:489857/1‑227 (MQ=255)                                                                                                                                                cACGGCGTGACCCAAAAACAGAAGGGGTAAAATAAACGGGAAGAAAAATCCCAGGGTTTG‑GTTAATTCCCACCAGCTGGGC‑GG‑TA‑GGGG‑GGGGGGA‑AATCAAC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTg                                                                                                                                                                                                         <  2:1321990/134‑1 (MQ=255)                                                                                                                                                     cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑GTGGGGG‑GAACAGC‑AAGCCCGCCACCTTTTCCC‑CATTTTTCCGGGCGGCCGTG‑AAGGACATCg                                                                                                                                                                                      >  1:926831/1‑147 (MQ=255)                                                                                                                                                          aCCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGG‑GTGGGGG‑GATAAGC‑AAACCCCCCCCTTTTTCCC‑GATTTTTCCGGGCGGCAGGG‑aa                                                                                                                                                                                              >  1:1321990/1‑134 (MQ=255)                                                                                                                                                                            aaaaTGAAAAGGAAGGAAATTTCCGGGGGTTG‑GTTAATTCCCCCCAGCGGGGC‑GG‑TT‑AGGG‑TGGGGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:489857/237‑1 (MQ=255)                                                                                                                                                                              aaTGAACAGGAAGG‑AAATTCCCAGGGTTCGGGTATTCCCCACCAGCAGGGC‑GG‑TT‑AGGG‑TGGGGGG‑AACCAGA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                <  1:1273341/239‑1 (MQ=255)                                                                                                                                                                                tGAACAGGAAGGTAAATCCCGGCGTTGC‑GGTAATCCCCACCAGCATGGC‑GT‑TA‑AGGG‑TGGCGGA‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:264189/239‑1 (MQ=255)                                                                                                                                                                                agAACGGGAGGGAAAATCCCGGGTGTTC‑GGTACTCCCCCCCGGCAGGGG‑GT‑TT‑AGGT‑TTGGGGA‑AAACAGC‑AGCCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2130806/238‑1 (MQ=255)                                                                                                                                                                                  aaCAGTAAGGTAAATCCCAGGGTTTG‑GGTAATCCCCCCCAGCAGGGG‑GG‑TA‑AAGG‑TTGGGGA‑AACCAGA‑AACCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTAc                                                                                                                                                                                          <  1:1045186/114‑1 (MQ=255)                                                                                                                                                                                  aaCAGAAAGGAAAATTCCAGGGGTTG‑GGTAATCCCCACCAGCAGGGG‑GG‑TT‑AAGG‑TGGGGGG‑AACCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCa                                                                                                                                                                                                    <  1:2060251/105‑1 (MQ=255)                                                                                                                                                                                       gAAGTTAAATTCCGGGGGTTC‑GGTAATCCCCCCAAGAAGGGC‑GG‑TT‑AAGG‑TGGGGGA‑AATAAGA‑AAGCCCGCCACTTGTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                         <  1:1850716/237‑1 (MQ=255)                                                                                                                                                                                       gAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑a                                                                                                                                                                                                                                                <  1:982943/58‑1 (MQ=255)                                                                                                                                                                                        aaGGTAAATCCCGGGGTTTC‑AGTAACTCCCCCCAGCAGGGC‑GG‑TA‑AAGG‑GTGCGGA‑AATAAAC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCa                                                            <  1:2230435/233‑1 (MQ=255)                                                                                                                                                                                            tAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑a                                                                                                                                                                                                                                                <  2:960365/53‑1 (MQ=255)                                                                                                                                                                                            tAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑a                                                                                                                                                                                                                                                >  1:960365/1‑53 (MQ=255)                                                                                                                                                                                            tAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑GAACAAC‑AAACCCGCCCCTTTTTCCC‑GGTTTTCCGGTGCGGCAGGG‑AAGTACATCGCCAAAg                                                                                                                                                                                >  2:1045186/1‑114 (MQ=255)                                                                                                                                                                                            tAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑GTGCGGG‑GAACAGC‑AAAACCGCCACTTTTACCC‑CATTTTTCGGGGCGGGAGGG‑aaagaaa                                                                                                                                                                                         >  2:2060251/1‑100 (MQ=255)                                                                                                                                                                                               tgtccggggtttC‑GATAATCCCCGCCGGTATGCC‑GG‑‑T‑AAGG‑TTGGGGGCAACCAGA‑AGGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                               <  1:1644543/229‑1 (MQ=255)                                                                                                                                                                                               aTTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑TAGG‑GTGCGGG‑GATTAGC‑AAAACCCCCACTTTTACCCGGATTTTCC‑GTGCGGCAGTG‑AAgga                                                                                                                                                                                           >  1:509552/1‑98 (MQ=255)                                                                                                                                                                                                          ttgggTAATCCCCACCAGATGGGC‑GG‑TT‑AAGG‑TTGCGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:399304/206‑1 (MQ=255)                                                                                                                                                                                                                   ccccACCAGCATGGC‑GG‑TT‑AGGG‑TGGCGGG‑AATCAGA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                               <  1:2099487/219‑1 (MQ=255)                                                                                                                                                                                                                      ccccaGGAGGGC‑GT‑TA‑AAGG‑TTGCGGA‑AATCAGC‑ACCCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  1:26792/237‑1 (MQ=255)                                                                                                                                                                                                                      ccccaGCAGGGC‑GG‑TT‑AGGT‑TGGCGAG‑AATCAGA‑ACGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  1:1211885/237‑1 (MQ=255)                                                                                                                                                                                                                             aTGGC‑GG‑TT‑AAGG‑TTGCGGG‑AAACAGC‑AAACCCCCCCCTTTTACCC‑CGTTTTTCGGTGCGGCAGAG‑AAGTACATCGCTAAAGAAAAGCCCACGGC‑ACCAAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAAAAACGATTTCACCGGCTCCGCGCCCCAACCGAACTCCTTCAAATGTGGCGCCCCGTACAACTGCATACCGCCtccta                                     >  2:2099487/1‑216 (MQ=255)                                                                                                                                                                                                                                  ggg‑TA‑AAGG‑TGGGGGG‑AACCAAC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa               <  1:2009331/237‑1 (MQ=255)                                                                                                                                                                                                                                     ttt‑AAGT‑TTGGGGG‑AACCACC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTg                                                                                                                                                                                                         <  2:948559/51‑1 (MQ=255)                                                                                                                                                                                                                                      tt‑AAGGTTTGGGGA‑AATCAAA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCgg                                                                                                                                                                                                      <  1:1596831/55‑1 (MQ=255)                                                                                                                                                                                                                                                 cGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑aa                                                                                                                                                                                              >  1:948559/1‑52 (MQ=255)                                                                                                                                                                                                                                                 cGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTa                                                                                                                                                                                           >  2:1596831/1‑55 (MQ=255)                                                                                                                                                                                                                                                 cGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCCGATGCTGCGCCCCGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCCAAGATACACTTGCCGTAGATg  >  2:2092916/1‑239 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TT‑AAGG‑TTGCGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862430‑3862902

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