Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,677	0	G	A	32.6%	29.7 / 13.4	48	pseudogene (927/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 G (11/18);  new base A (14/0);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 5.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAGG‑TTGCGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGG  >  NC_000913/3862457‑3862908                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAAG‑GTGCGGG‑GATAAGC‑CAACCCGCCCCCTTTTccc                                                                                                                                                                                                                            >  1:1606364/1‑237 (MQ=255)                tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑GTTAAG‑TTGCGGG‑GAATAGCAAAGCCCCCCCTTTTTACCC‑AATTTTCggggcgggc                                                                                                                                                                                                           >  1:990989/1‑233 (MQ=255)                         aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑GTAAGG‑GTGCGGG‑GATCAGC‑AAGCCCGCCCCTTTTACCC‑GATTTTTCCGTGCGGGCGTG‑AAGTAc                                                                                                                                                                                                >  1:198988/1‑239 (MQ=255)                                 gaaaacccatAGACAACCAGCCCCCCGCCCTTTTCCCCCACCACAA‑AAAAGAGATAAATACCGGTCGATAGAGAGCCCGCCAGCAGGGCGTGCACCAAAAAAAGAAAGGGAAAAAGAAACGGGAAGAAAAATCCCAGCGGTTC‑GATAATCCCCCCCAGCGGGGG‑GG‑TTAGGG‑TGGGGGA‑AACCAAA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTa                                                                                                                                                                                         <  2:2195119/229‑1 (MQ=255)                                                      ccccccgccATTTTGCCCCCCACCCAAAAAGATAAACATCACGGTCGACATGGAGCCCGCAAGCCCGGCGTGTCCCAAAAACAGAAAGGTAAAAATAAAAAGGAAGGTAAATTCCAGCGGTTC‑GTTAACCCCCCCCAGCAGGGC‑GG‑TTAGGG‑TGGCGGA‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCaa                                                                                                                                                                   <  2:198988/232‑1 (MQ=255)                                                               aTGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGCGGG‑TTAGGG‑TGGCGGG‑GA‑CACC‑AAACCCCCCCCATTTACCC‑CCGTTTTCTGTGGCGGAGGG‑GAGGAAAACCCCAAAGAAAAGTCCACGGC‑GCCAAACAatttttt                                                                                                                                                         >  2:798189/1‑237 (MQ=255)                                                               aTGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGCGGG‑TTAAGG‑‑TGCGGG‑GATCAGC‑AAACCCGCCCCCTTTACCC‑CATTTTTCGGTGCGGCAGTG‑AAGTACa                                                                                                                                                                                               >  1:2124034/1‑202 (MQ=255)                                                                                     aaaTACATTAGGGCCAAAATGGGGCCCGCCAGCGGGGCGTGCACCAAAAACACCAACGGTAAAATGAAGAGGAAGTTAAATCCCGGCGTTTG‑GGTAATCCCCCCCAGCAGGGG‑GG‑TAAGGG‑TGGGGGA‑AACCAAA‑AACCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTg                                                                                                                                                                                                       <  2:204924/172‑1 (MQ=255)                                                                                             tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGCGGG‑TAAAGG‑TTGGGGG‑AATAACC‑AAGCCCCCCCCTTTTTCCC‑GATTTTTCGTGGCGGCAG‑G‑GAGAAAAACCCCCAAAAAAAAACCACCGC‑ACCCAAAAATTTTTTGTTACCCTGCAGAGCAAAAACGCCCTccc                                                                                                                            >  2:834463/1‑238 (MQ=255)                                                                                             tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GGTTTAAGG‑TTG‑GGG‑GATCAGC‑AAACCCCCCACCTTTACCC‑CGTTTTTCGGGGCGGGAGGG‑aa                                                                                                                                                                                                    >  1:1432316/1‑168 (MQ=255)                                                                                               cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTTAGG‑GTGCGGG‑GAACAGC‑AAACCCCCCCCCTTTTCCC‑GATTTTTCGGGGCGGCAGGG‑AAggaaaacg                                                                                                                                                                                            >  1:204924/1‑167 (MQ=255)                                                                                               cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑ATAAAG‑GTGCGGG‑AAGCAGC‑AAAGCCGCCCCCTTTTCCC‑CCTTTTTCTGTGCGGGAGTG‑GAGGAAAACGCCCAAAAAAAGACCCCCGC‑ACCCAAGAATTTTTAGTTACCGTGCAGGGGAAAAACGccctaccgg                                                                                                                          >  2:878734/1‑232 (MQ=255)                                                                                                    gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑GTTAGG‑GTGGGGG‑GAACACC‑CAAGCCGCCCCCTTTTCCC‑CACTTTTCGGTGCAGCCGGGAAAGAACAAC‑CCAAAAAAAAGACCCCCGC‑ACCAAAAAATTTTTTGTGTCCGTGCATGGCAAACACACCTTCCggggcaaa                                                                                                                     >  2:1273341/1‑234 (MQ=255)                                                                                                      ctttGGGGCCGCCAGCGCGGCGTATACCAAAAACAGCAAGGGAGAAAAGAAGAGGAAGGTAAATTCCAGCGGTTC‑GGTAATCCCCCCCAGCGGGGG‑GG‑TTAGGG‑TGGGGGG‑AAACAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg                                                                                                                   <  2:1175135/237‑1 (MQ=255)                                                                                                      catgggggCCGCCAGCGCGGCGTGCACCGAAAACAGCAAGGGTAAAAGGAAGAGGAAGGAAAATTCCGGGGGTTG‑GGTAATCCCCCCCAGCGGGGG‑GT‑TTAGGG‑TGGGGGG‑AAACAGC‑AACCCCGCCACTTTTACCC‑GATTTTCCGGTGCgg                                                                                                                                                                                                            <  2:455342/146‑1 (MQ=255)                                                                                                              ccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAGG‑GTGCGGG‑GATAAGCAAAGCCCCCCCCTTTTCCCC‑ATTTTTCCGG‑GCGGGGGGG‑GAGGAAAACGCCAAAAAAAAACCCCCCGCGACCAAA‑AATTTTGTGGTACCGGGCAGGGCAAAAACACCCTCCCGGGCAAAAGATTTCaac                                                                                                           >  2:327369/1‑237 (MQ=255)                                                                                                                gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAAG‑TTGCGGG‑GAACCGC‑AAACCCCCCCCCTTTTCCC‑GATTTTTCCGGGGGGCAGGG‑GAGTa                                                                                                                                                                                                 >  1:455342/1‑152 (MQ=255)                                                                                                                 ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTTAGG‑GTGCGGG‑GAACAGC‑AAACCCCCCCCCTTTACCC‑CCTTTTTCGGGGGGGGCGGG‑GAGGACACCCCCAAAAAAAAAACCCCGGC‑ACCCAAAAATTTTTTGTTTACGTGCGTGGCCAAACACCCCTCCGGGgaaaacaagtttatagg                                                                                                         >  2:489857/1‑227 (MQ=255)                                                                                                                     cACGGCGTGACCCAAAAACAGAAGGGGTAAAATAAACGGGAAGAAAAATCCCAGGGTTTG‑GTTAATTCCCACCAGCTGGGC‑GG‑TAGGGG‑GGGGGGA‑AATCAAC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTg                                                                                                                                                                                                               <  2:1321990/134‑1 (MQ=255)                                                                                                                          cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAGG‑GTGGGGG‑GAACAGC‑AAGCCCGCCACCTTTTCCC‑CATTTTTCCGGGCGGCCGTG‑AAGGACATCg                                                                                                                                                                                            >  1:926831/1‑147 (MQ=255)                                                                                                                               aCCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTTAGG‑GTGGGGG‑GATAAGC‑AAACCCCCCCCTTTTTCCC‑GATTTTTCCGGGCGGCAGGG‑aa                                                                                                                                                                                                    >  1:1321990/1‑134 (MQ=255)                                                                                                                                                 aaaaTGAAAAGGAAGGAAATTTCCGGGGGTTG‑GTTAATTCCCCCCAGCGGGGC‑GG‑TTAGGG‑TGGGGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                         <  1:489857/237‑1 (MQ=255)                                                                                                                                                   aaTGAACAGGAAGG‑AAATTCCCAGGGTTCGGGTATTCCCCACCAGCAGGGC‑GG‑TTAGGG‑TGGGGGG‑AACCAGA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                      <  1:1273341/239‑1 (MQ=255)                                                                                                                                                     tGAACAGGAAGGTAAATCCCGGCGTTGC‑GGTAATCCCCACCAGCATGGC‑GT‑TAAGGG‑TGGCGGA‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:264189/239‑1 (MQ=255)                                                                                                                                                     agAACGGGAGGGAAAATCCCGGGTGTTC‑GGTACTCCCCCCCGGCAGGGG‑GT‑TTAGGT‑TTGGGGA‑AAACAGC‑AGCCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:2130806/238‑1 (MQ=255)                                                                                                                                                       aaCAGTAAGGTAAATCCCAGGGTTTG‑GGTAATCCCCCCCAGCAGGGG‑GG‑TAAAGG‑TTGGGGA‑AACCAGA‑AACCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTAc                                                                                                                                                                                                <  1:1045186/114‑1 (MQ=255)                                                                                                                                                       aaCAGAAAGGAAAATTCCAGGGGTTG‑GGTAATCCCCACCAGCAGGGG‑GG‑TTAAGG‑TGGGGGG‑AACCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCa                                                                                                                                                                                                          <  1:2060251/105‑1 (MQ=255)                                                                                                                                                            gAAGTTAAATTCCGGGGGTTC‑GGTAATCCCCCCAAGAAGGGC‑GG‑TTAAGG‑TGGGGGA‑AATAAGA‑AAGCCCGCCACTTGTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                               <  1:1850716/237‑1 (MQ=255)                                                                                                                                                             aaGGTAAATCCCGGGGTTTC‑AGTAACTCCCCCCAGCAGGGC‑GG‑TAAAGG‑GTGCGGA‑AATAAAC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCa                                                                  <  1:2230435/233‑1 (MQ=255)                                                                                                                                                                 tAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAGG‑TTGCGGG‑GAACAAC‑AAACCCGCCCCTTTTTCCC‑GGTTTTCCGGTGCGGCAGGG‑AAGTACATCGCCAAAg                                                                                                                                                                                      >  2:1045186/1‑114 (MQ=255)                                                                                                                                                                 tAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAGG‑GTGCGGG‑GAACAGC‑AAAACCGCCACTTTTACCC‑CATTTTTCGGGGCGGGAGGG‑aaagaaa                                                                                                                                                                                               >  2:2060251/1‑100 (MQ=255)                                                                                                                                                                    tgtccggggtttC‑GATAATCCCCGCCGGTATGCC‑GG‑‑TAAGG‑TTGGGGGCAACCAGA‑AGGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                     <  1:1644543/229‑1 (MQ=255)                                                                                                                                                                    aTTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTTAGG‑GTGCGGG‑GATTAGC‑AAAACCCCCACTTTTACCCGGATTTTCC‑GTGCGGCAGTG‑AAgga                                                                                                                                                                                                 >  1:509552/1‑98 (MQ=255)                                                                                                                                                                               ttgggTAATCCCCACCAGATGGGC‑GG‑TTAAGG‑TTGCGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                         <  1:399304/206‑1 (MQ=255)                                                                                                                                                                                        ccccACCAGCATGGC‑GG‑TTAGGG‑TGGCGGG‑AATCAGA‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                     <  1:2099487/219‑1 (MQ=255)                                                                                                                                                                                           ccccaGGAGGGC‑GT‑TAAAGG‑TTGCGGA‑AATCAGC‑ACCCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                               <  1:26792/237‑1 (MQ=255)                                                                                                                                                                                           ccccaGCAGGGC‑GG‑TTAGGT‑TGGCGAG‑AATCAGA‑ACGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                               <  1:1211885/237‑1 (MQ=255)                                                                                                                                                                                                  aTGGC‑GG‑TTAAGG‑TTGCGGG‑AAACAGC‑AAACCCCCCCCTTTTACCC‑CGTTTTTCGGTGCGGCAGAG‑AAGTACATCGCTAAAGAAAAGCCCACGGC‑ACCAAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAAAAACGATTTCACCGGCTCCGCGCCCCAACCGAACTCCTTCAAATGTGGCGCCCCGTACAACTGCATACCGCCtccta                                           >  2:2099487/1‑216 (MQ=255)                                                                                                                                                                                                       ggg‑TAAAGG‑TGGGGGG‑AACCAAC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                     <  1:2009331/237‑1 (MQ=255)                                                                                                                                                                                                          tttAAGT‑TTGGGGG‑AACCACC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTg                                                                                                                                                                                                               <  2:948559/51‑1 (MQ=255)                                                                                                                                                                                                           ttAAGGTTTGGGGA‑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)                                                                                                                                                                                                                      ggg‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGa       >  1:490128/1‑239 (MQ=255)                                                                                                                                                                                                                          aaTCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg     <  1:392238/238‑1 (MQ=255)                                                                                                                                                                                                                            tCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtgg  >  2:429475/1‑239 (MQ=255)                                                                                                                                                                                                                            tCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCGCAAACTGAACTCCTGCATATGCGGCGCCCAGTACATCTGAATGCCGCCTTCAACACGTGCCCGGCCAAAGATAAACTGCTCGTAGATGAAgtgg  >  2:266503/1‑239 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           CGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACAC‑CAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTC‑GGTAATTCCCACCAGCATGGC‑GG‑TTAAGG‑TTGCGGG‑AATCAGC‑AAGCCCGCCACTTTTACCC‑GATTTTCCGGTGCGGCAGTG‑AAGTACATCGCTAAAGAAATGCCCACGGC‑ACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGG  >  NC_000913/3862457‑3862908

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