Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,667	0	G	A	38.1%	42.6 / 15.0	42	pseudogene (937/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 (18/8);  new base A (0/16);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 4.54e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGTTGCGG‑G‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862429‑3862902                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 aTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGTTGCgg‑g                                                                                                                                                                                                                                                  <  1:372146/239‑1 (MQ=255)       aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑GTAAAGTTGGGG‑G‑GATCa                                                                                                                                                                                                                                            >  1:291438/1‑238 (MQ=255)        tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATC‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAAGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑GTAAAGGTGGGG‑G‑GATCaat                                                                                                                                                                                                                                          >  1:934835/1‑237 (MQ=255)                      atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAAGTGGCGG‑G‑AAACAGC‑AAACCCGCCCCttt                                                                                                                                                                                                                           >  1:223957/1‑239 (MQ=255)                                 gTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGGTCGGT‑AATTCCCACCAGCATGGC‑GG‑GTAAGGGTGCGG‑G‑GGTCAGC‑AAGCCCGCCCCTTTTTTCCGActtt                                                                                                                                                                                                                >  2:2228996/1‑239 (MQ=255)                                          cATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GGTTTAAGGTTGGGG‑G‑AAAAGGA‑AAGCCCGCCCCTTTTCCCCGATTTTCCgggtgg                                                                                                                                                                                                        >  1:1905859/1‑235 (MQ=255)                                          cATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAAGGTTCGG‑G‑GAACAGC‑AAACCCGCCACCTTTTTCCCCTTTTTCGCTGCg                                                                                                                                                                                                        >  1:1328351/1‑238 (MQ=255)                                           actggggtactgTTTTGCGGAAAAAACCT‑AGAAAACCCCCCCCCCGCCCTTGTCCCCCCCCAC‑CCAAAAAAAACCATT‑TCGGGTGAACATTGGGCCCGCCGGCACGGGGGTCCCCGAAAACAACAACGGAAAAAAAACAAGGAAGGTAAATCCCGGGGGTTCGGT‑AACCCCCCCCAGCAGGGC‑GT‑TTAAGGTGGGGG‑A‑AAACAGC‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCgg                                                                                                                                                                                                       <  2:316389/228‑1 (MQ=255)                                            gcgacaTCCATTTTTGCGA‑AAAACCCT‑AGACATCCACACCCCCGCCTATGTCCCCCCCCCCACCAAAGAGATAATTT‑AC‑GGGCGACAGTGGCGCCCCCAGCCAGTCTTGGACCGAAAACGGAAACGTAAAAAAGAGGAGGAA‑GAAAAGTCCCGGCGTGAGGTAAATCCCCACAAGGTGGGG‑GG‑TTATGGGTGGGG‑G‑AGAAAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                     <  2:403557/234‑1 (MQ=255)                                                atcccAGTTTTGCGG‑AAAACCAT‑AGTCACTCACACCCCCCCCTATTTCCCCCCCCACACAAAAG‑GAAACATTGGC‑GGACGACAGTGAGGCCGCCAGCGCGGCGTGTCCCGCAACCAGCAACGGAAAAATGAAAAGAAAGGAAAATCCCAGCGGTTCGGT‑ACTCCCCACCAGAATGGG‑GT‑TTAGGGTGGGGG‑A‑AACCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTg                                                                                                                                                                                                  <  2:1108314/235‑1 (MQ=255)                                                 gatcAGTTTTGCGG‑AAAAAACTAAGTCACCCAGCCCCCCCCCAATGTCCCCCCCCCCCAAAAAGAAATCAATT‑AC‑GGTCACAATGGAGGCCGCCAGCCGGGGGTGTACCGAAAACAGAAAGGGTAAAATGAACAGAAGGATAAATCCCGGCGGTTCGGT‑AATCCCCACCGGCGGGGG‑GT‑TTAGGGTGGGGG‑A‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                     <  2:1476639/232‑1 (MQ=255)                                                             gg‑TAAAACCT‑AGTAAATCAGCCCTCCCCCCATTTGGCCCACCACCCCAAAGGGTACCATT‑GC‑GGTCGCCATTGAGGCCGCCAGCACGGCGTTCACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAGGTTTGCGG‑G‑AATCAGC‑AAGCCTGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑a                                                                                                                                                                                   <  2:1654495/238‑1 (MQ=255)                                                                           tCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGCGGG‑TTAAGG‑TGCGG‑G‑GATCAGCAAAACCCCCCCCTTTTCCCCACTTTTCCGGT‑CGGCATGGAGG‑TACATCTCTCAAAAAAATTccc                                                                                                                                                                        >  2:112242/1‑238 (MQ=255)                                                                                       ccGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GGTTTAAGGGTGGGG‑G‑AA‑AAGC‑AAACCCCCCACCTTTTCCCCAtttttcc                                                                                                                                                                                                             >  1:1512463/1‑187 (MQ=255)                                                                                         ccccttgtgccccACAACACAAAAGAGATCAATG‑GC‑GGTCGACATGGAGGCCCCCAGCACGGGGGTTACCGAAAACAGCAACGGAAAAAGAAACAGGAAGGAAAATTCGAGCGTTTCGGT‑AACCCCCCCCAACAGGGG‑GT‑TTAGGGTTGCGG‑A‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGAttt                                                                                                                                                        <  2:456380/229‑1 (MQ=255)                                                                                                     caccacaccaAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGTTGCgg‑g                                                                                                                                                                                                                                                  <  1:427734/141‑1 (MQ=255)                                                                                                        cacaccaAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGTTGCGG‑G‑a                                                                                                                                                                                                                                                <  1:2117144/139‑1 (MQ=255)                                                                                                        cacaccaAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAAAGAACAGGAAGGTAAATTCCAGCGGTTCGGG‑AATTCCCACCAGCATGGG‑GG‑TATAAGTAATGG‑G‑CAGCCGC‑ACACCCGCCCCTTTTTCTCGAAATTCCGTGGCGGCGGGGAAG‑AAGAACGCT‑AAAAAAAATACCACGGCGACACAAATATTTGAATTAcacagac                                                                                                                                         >  2:379407/1‑233 (MQ=255)                                                                                                               aaGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTTAGGGTTCGG‑G‑GATCAGccaa                                                                                                                                                                                                                                       >  2:427734/1‑138 (MQ=255)                                                                                                                  agaTACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGCGGG‑TTAAGGTTGGGG‑G‑AATAAGC‑AACCCCCCCCCTTTTACCCGATTTTCCGGGGGGGCGGGGAGG‑AAAATGCCTAAAAAAAATCCCCAC‑ACACCAAAAAATTTTTAGTTACCGAGCA‑GGGCcaaa                                                                                                                                >  2:2123212/1‑238 (MQ=255)                                                                                                                  agaTACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑GTTAGGGTGCGG‑G‑GAACAGC‑aaaa                                                                                                                                                                                                                                     >  2:2117144/1‑137 (MQ=255)                                                                                                                                   gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGGGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑GTAAAGTTGGGG‑G‑GATCCCC‑CACACCGCCCCCTTTTCCCCATTTTCCCGGGGGGGCGGGGAG‑AACAACCCT‑AAAAAAAATTCCCCGCGACCAAAAATGTTTTAGTTGTCGTGCACGGGGAAAACGCCCTTCCGCGacca                                                                                                                >  2:404150/1‑239 (MQ=255)                                                                                                                                          aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGCGGG‑TTAAGGTT‑CGG‑G‑GAAAAGC‑CAGCCCCCCCCCTTTTTCCCCTTTTTCGGGGGGGGAGGGAAGAAACACGCCA‑AAAAAAAACCCCCCGGCCCCAACAAATTTT‑AGTTTCCCCGCAGGGGAAAAACCCCCTCCCGGGaaaaaca                                                                                                             >  2:2054666/1‑234 (MQ=255)                                                                                                                                                               ccGAAACCAGCACGGGTAAAATAAACAGAAAGTAAAATTCCGGCGTTTCGTT‑AATCCCCACCAGCTTGGC‑GT‑TTAAGTTTGCGG‑A‑AACCACC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTg                                                                                   <  1:365922/239‑1 (MQ=255)                                                                                                                                                                             gttaaaagcaaCAGGAGGGAGAATTCCGGCGGTTGAGA‑TGTTCCCACCGGCGTGGG‑TA‑TTAGGGTGGGGG‑G‑AATCAGC‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGGTCGGCGCTCAAACTGAACTCTTGCAg                                                                        <  2:2202672/227‑1 (MQ=255)                                                                                                                                                                                gAAATGAAGAGGAAGGTAAATCCCTGGGTTGCGTG‑AACTCCCACCAGGGTGGC‑GG‑TTAGGGTGGGGG‑A‑AACCAGA‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:442161/238‑1 (MQ=255)                                                                                                                                                                                aaattGAACAGGAAGGAAACTCCCGGCGGTTCGTT‑AATCCCCACCAGCGTGGC‑GG‑TTAGGGTGGCGG‑A‑AATCAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:305533/234‑1 (MQ=255)                                                                                                                                                                                 aaaTAAAAAGGAAGGAAAATCCCGGCGTTTCGGT‑ACTCCCCCCCAGAAGGGC‑‑G‑TTTAGGTTGGGG‑GAAACCACC‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTg                                                                                                                                                     <  1:154012/156‑1 (MQ=255)                                                                                                                                                                                      acagggaaggaaACTTCCGGGGGGTCGG‑‑AATCTCCCCCAGGATGGG‑AG‑GTAAGGGTGGGGTA‑AATCAGC‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  2:1354589/233‑1 (MQ=255)                                                                                                                                                                                       aaaggaggGTAAATGCCGGGGGTTGGAT‑AATCCCCCCCAGCATGGG‑GT‑TAAAGGTTGGGG‑A‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                           <  1:627355/237‑1 (MQ=255)                                                                                                                                                                                           gAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGTTGGGG‑G‑AATCAGC‑AAGCCCGCCACTTTTTCCCGATTTTCCCGGGCGGCAGTGAAG‑GACATCGCT‑AAAGAAATGCCCCCGGCACCCAAGATTTTTTAGTTACCGTg                                                                                                                                           >  2:154012/1‑156 (MQ=255)                                                                                                                                                                                           gAAGGAAAATCCCGGCGGTTGGGT‑AATCCCCACCAGCAGGGG‑GG‑TTAGGGTTGCGG‑G‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACAt                                                                                                                                                                                         <  1:1489754/111‑1 (MQ=255)                                                                                                                                                                                            aaGGTAATTCCCAGCGGTTCGGA‑ATTCCCCACCAGCAGGGC‑GG‑TTAAGGTGGGGG‑G‑AATAAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:1547480/226‑1 (MQ=255)                                                                                                                                                                                                aaaaTTCCGGTGGTTCGGT‑AACCCCCCCCAGCGGGGG‑GG‑TAAGGGTGGGGG‑A‑AACCAGA‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGt                                                          <  1:1615177/230‑1 (MQ=255)                                                                                                                                                                                                   ttcccggCGGTTGGGT‑AATTCCCCCCAGTGTGGG‑GG‑TTAGGGTGGGGG‑A‑AATCAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                               <  1:126137/233‑1 (MQ=255)                                                                                                                                                                                                   attcgggCGGTTCGGT‑ATTTCCCACAAGTGGGGC‑GG‑TTAAGGTGGGGG‑A‑AATCACC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                               <  1:139042/233‑1 (MQ=255)                                                                                                                                                                                                    ccccgggggTTGGGT‑AATCCCCCCCGGAGGGGG‑GG‑TTAGGGTGGGGG‑G‑AACCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgc                                              <  1:2054666/232‑1 (MQ=255)                                                                                                                                                                                                     tCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGGTGCGG‑G‑GATCAGC‑AAACCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGGGAAG‑GACAACCCT‑AAAGaa                                                                                                                                                                              >  2:1489754/1‑111 (MQ=255)                                                                                                                                                                                                              tcagT‑AATTCCCACCAGCGGGGC‑GT‑TTAAGGGTGCGG‑A‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCaa                                     <  1:2039007/235‑1 (MQ=255)                                                                                                                                                                                                                          ccctgccggggC‑GT‑TTGAGTTGGCGG‑A‑AATCAGC‑GACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  2:658647/230‑1 (MQ=255)                                                                                                                                                                                                                          ccccaGCTGGGG‑GG‑TTAAGGTTGGGG‑A‑AAACAGC‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  1:2123212/237‑1 (MQ=255)                                                                                                                                                                                                                          ccacaacgtGGC‑GG‑TTGGGGTTGGGG‑A‑ACACAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                         <  1:1406220/231‑1 (MQ=255)                                                                                                                                                                                                                                                tgcggg‑g‑AACCAGC‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt       <  1:1387541/233‑1 (MQ=255)                                                                                                                                                                                                                                                   cGG‑G‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGcc                                                                                                                                                                         >  1:808480/1‑73 (MQ=255)                                                                                                                                                                                                                                                   cGG‑G‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  >  1:1146295/1‑239 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGG‑TAAAACCT‑AGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATT‑AC‑GGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGT‑AATTCCCACCAGCATGGC‑GG‑TTAAGGTTGCGG‑G‑AATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAG‑TACATCGCT‑AAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862429‑3862902

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