Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,664	0	C	G	41.0%	28.5 / 12.5	40	A314P (GCA→CCA)	glvC	putative PTS enzyme II component GlvC
Reads supporting (aligned to +/- strand):  ref base C (10/13);  new base G (5/11);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTG  >  NC_000913/3862419‑3862910                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                gCACAACCCAATGGCGATTTGGGTCAGTATCATATCCGCGTAGTTGATGAACATCGGGATCCAGTTTTGCGGTAAAAACTAGTCAATCAGACCTCCGCCCATGTTGCCCAACACACCAAAGAGATCCA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTAACGCAAACAGCAACGGTGAAATGCACCGGAAGGAAAATTCAAGCGGCTCGGTAATTCCCACCCGCAAGGCGG‑TTCAGGTTGGGGGa                                                                                                                                                                                                                                                             >  1:389384/1‑250 (MQ=255)        ccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTTAAGTTGCGGGGATCA‑Gcc                                                                                                                                                                                                                                                     >  1:164304/1‑249 (MQ=255)        ccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCaa                                                                                                                                                                                                                                                    <  2:37968/251‑1 (MQ=255)        ccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCaa                                                                                                                                                                                                                                                    <  2:240648/251‑1 (MQ=255)           aTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCAAGcc                                                                                                                                                                                                                                                 <  2:455433/251‑1 (MQ=255)               ggATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTGGCGGTAAAACCGAGTAAATCAGACCGCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCAGGGCGG‑TTAAGGTTGCGGGAATCA‑GCAAGCCC‑Gcc                                                                                                                                                                                                                                            <  2:191892/250‑1 (MQ=255)               cGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCAAGCCC‑Gcc                                                                                                                                                                                                                                            <  1:231402/251‑1 (MQ=255)                gATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAAGGTGCGGGGATCA‑GCAAGccccccc                                                                                                                                                                                                                                            >  1:574099/1‑247 (MQ=255)                         cAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑GTAAGGGTGGGGGAATCA‑ACAAGCCC‑GCCCC‑CTTT‑A‑ccc                                                                                                                                                                                                                               >  2:124240/1‑251 (MQ=255)                            catcatACCCGCGTGGTTGCTGAACATCGGGATCCATTTTTGCGGAAAAACCTAGTCAATCAGCCCCCCGCCAATGTTGCCCCCCACACAAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATTCCAGCGTTTCGGTAATTCCCCCCAGCATGGGGG‑TTAAGGTGGCGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGAt                                                                                                                                                                                                                            <  1:175531/251‑1 (MQ=255)                                    ccgggtggtttctaaccatggggACCCAGTTTGGCGATAAACCCTAGTAAATCAGACCCCCGCCATTGTTGCCCCCCACCCAAAAGAGATACA‑TTCCGGTAGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAAGGGGAAAATAAACAGAAAGGTAAATCCCGGGGTTTCGTTATTCCCCCCCAGCAGGGCGT‑TTAGGGTTGGGGGAATCA‑GCAACCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGt                                                                                                                                                                                                                    <  2:239075/231‑1 (MQ=255)                                     cgcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTTAGGTTGCGGGGAACA‑GCAAGCCC‑ACCAC‑CTTT‑A‑CCCGCtttttccgtg                                                                                                                                                                                                                   >  2:481203/1‑245 (MQ=255)                                        gTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTTCGGGGATCC‑GCAAGCCC‑CCCCC‑TTTT‑T‑CCCGCTTTTCCGGTTCgg                                                                                                                                                                                                                >  1:581012/1‑251 (MQ=255)                                         tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGGTTCGGGAATCA‑GCAAACCC‑GCCAC‑TTTT‑A‑CCCGCtttttccgggcggc                                                                                                                                                                                                               >  1:416806/1‑241 (MQ=255)                                                            tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGGGGGGAA‑CC‑GCAAACCC‑GCCCCCTTTT‑A‑CCCCTTTTTCCgggggggcgggaaggaacaccccca                                                                                                                                                                                              >  2:550093/1‑227 (MQ=255)                                                             ccattttggCGGAAAACCCTAGTCAATCAGACCCCCGCCCATGTGGCCCACCACACAAAAGAAATACA‑TTACGGTCGACTTGGGGGCCGCCAGCCCGGCGTGTACCACAACCAGCAACGGTGAAATGAACAGGAAGGAAAATTCCGGGGTTTCGGTAATTCCCCCCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAg                                                                                                                                                                                         <  1:221366/243‑1 (MQ=255)                                                               aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTT‑CGGGGATCA‑GCAAACCC‑GCCCC‑CTTT‑T‑CCCCCTTTTTCGGTGCGGCAGTGAAGTACA‑ACGcctaaaa                                                                                                                                                                                         >  2:220394/1‑244 (MQ=255)                                                                                  gccaaccacccccccGCCCATGTCCCCCACCACACAAAAGAAATACA‑TAAGGGCAGACATTGAGGCCGCCGGCGGGGCGTGTCCCGAAAACAGCAACGGTAAAAGGAAAGGGAAGAAAAATTCCGGGGGTTCGGTAATCCCCCCCAGCAGGGGGG‑TAAAGGTGGGGGAAACCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGa                                                                                                                                                                   <  2:574099/238‑1 (MQ=255)                                                                                       acacccccccGCCCATGTTGCCCACCACCCAAAAGAGA‑ACAATTACGGTCGACATTGGGCCCGCCAGCGCGGGGTGTACCGAAAACAGCAAGGGTGAAATGAACAGGAAGTTAAATTCCGGGGGTTCGGTAATTCCCCCCAGCAGGGCGG‑TAAAGGTGGCGGAAATCA‑GAAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGAttttt                                                                                                                                                              <  1:6139/243‑1 (MQ=255)                                                                                                        ttGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGGGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGGGCGGCAGTGAAGGACA‑TCGC‑TAAAGAAATGCCCCCGGCACC‑AAAAATTTTTTAATTACCGTGGAgggg                                                                                                                                             >  1:57333/1‑250 (MQ=255)                                                                                                                             aaCA‑TTACGGCAGACTTGGAGGCCGCAAGCACGGCGTGAACCGCAAACAGCAACGGAAAAATGAAAAGGAAGGAAAATACCAGGGTTTCGGTATTTCCCCCCAGCAGGGCGG‑TTAAGGTGGGGGAAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGaacaa                                                                                                                       <  1:171340/250‑1 (MQ=255)                                                                                                                                      ggTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑GTAAGGGTGCGGGAAGCA‑GCAAAACC‑GCCCC‑TTTTTA‑CCCAATTTTC‑GGGGGGGGAGGGAAATACA‑TCGCCTAAA‑AAATTCCCCAGGCACCCAAAAATTTTT‑TGTTTCCGTGCAGGGCAAAAACCCCTTCCCGGGaaaacggtttt                                                                                                                 >  2:320950/1‑241 (MQ=255)                                                                                                                                      ggTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGGAATTCCCACCAGCATGGCGG‑TTAAGGGTGCGGGAATCA‑GCAAACCC‑CCCAC‑TTTT‑AACCCAATTTTC‑GGTGCGGGAGTGAAGTACAATCCC‑TAAA‑AAAAGCCCACGGCCCC‑AACGAATTTTGAGTTAACCTGCAGAGCaaaaaacccctccggggacaaaaatttc                                                                                                                 >  1:308631/1‑226 (MQ=255)                                                                                                                                          gACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAATGGTGCGGGAATCA‑GCAAGCCC‑GCCAC‑CTTT‑T‑CCCCATTTTCCGGGGCGGCCGGGAAGTACC‑TCCC‑TAAAAAAATGCCCACCGCACC‑AAAGATTTTTTAGTTACCGGGCAGGGCAAAACCCCCTTCC‑GGGCACAACGATTTCACCgg                                                                                                           >  1:114517/1‑251 (MQ=255)                                                                                                                                                      gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑GTAAAGGTGGGGGAAACA‑GCAAACCC‑CCCCC‑TTTT‑T‑CCCGATTTTTCGGGGCGGCAGTGAAGGACA‑ACGC‑TAAAAAAATGCCCCCGGCACC‑AAAAATTTTTTAGTTACCGGGCAGGGCAAAACCGCCTTCC‑GGGGAAAACGATTTCAGCGGCCCGGCGCTcc                                                                                                >  1:169101/1‑249 (MQ=255)                                                                                                                                                       ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGGTGCGGGGATCA‑GCAAACCC‑CCCCC‑TTTT‑A‑CCCCATTTTTCGCGGCGGCCGTGAAGGACA‑ACGC‑TAAAGAAAAGCCCACGGCCCC‑CAAAATTTTTTAGTTAACCTGCCGGGCCAAACCGCCTTCC‑CGGGACAACGAATTCAGCGGcgccgcgccccac                                                                                              >  2:236668/1‑239 (MQ=255)                                                                                                                                                            acggcggtcaCCGAAAACAGCAACGGGAAAATGAACAGGAAGGTAAATTCCGGCGGTTCGTTAATCCCCACCAGAAGGGCGG‑TAAAGGTGGCGGGAAGCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                         <  1:198007/242‑1 (MQ=255)                                                                                                                                                             gggCGTGTACCGAAAACGGAAAGGGAGAAAAAAACGGGAAGATAAATCCCCGGGGTTCGGTAATCCCCACCAGCAGGGGGT‑TTAAGGT‑GGGGGAAACAAGCACGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAc                                                                                        <  2:419096/250‑1 (MQ=255)                                                                                                                                                             gggCGTGTACCGAAAACAGCAGCGGTGAAAGGAAGGGGGGGGTAATTGCCAGCGGTGGGGAAATCCCCCCCAGGTGGGCGG‑TTAAGGGGGGGGGAAAGA‑GCAACGCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGAGTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAc                                                                                        <  2:189149/250‑1 (MQ=255)                                                                                                                                                                cGTGAACCGCAAACAGAAAGGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGTTAATCCCCCCCAGCAGGGCGG‑TAAAGGTGGGGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAActct                                                                                     <  1:94059/251‑1 (MQ=255)                                                                                                                                                                           accAGCAACGTTGAAATAAACAGGAAGTTAATTTCCAGCGTTTCGGTAATCCCCCCCAGGATGGTGT‑TAAAGTTTGGGGAAATCG‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgctg                                                                          <  1:514366/249‑1 (MQ=255)                                                                                                                                                                                  aCGGTAAAATGAACAGGAAGGTAAATCCCGGCGGTTCGGTATTCCCCCCCAGCATGGCGG‑TTAAGGTTGCGGAAACCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                                   <  1:96915/251‑1 (MQ=255)                                                                                                                                                                                             aaCAGAAAGGTAATTCCCGGCGTTTCGGTATTCCCCCCCAGCAGGGCTG‑TTAAGGTTGGGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAAt                                                        <  1:244187/251‑1 (MQ=255)                                                                                                                                                                                                   aaggaaaatcccggggtttCGGTAATCCCCCCCAGCATGGGGG‑TAAAGGTGGGGGGAAACA‑GAAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                                  <  1:550093/234‑1 (MQ=255)                                                                                                                                                                                                   aaGGTAAATCCCAGGGTTTCGGTAATCCCCACCAGAAGGGCGT‑TTAAGGTTGCGGGAACCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                                  <  1:432245/251‑1 (MQ=255)                                                                                                                                                                                                     gtaaaatcccgggggtgggtaaatccccACCAGCTGGGCGG‑TTAAGGTGGCGGAAATAA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCtt                                                <  1:222730/226‑1 (MQ=255)                                                                                                                                                                                                       tAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAAGGTGAGGGAATCA‑CCAAACCC‑GACCC‑TTTT‑A‑CACCATTTTCCGGTGCGGGCGGGAAGTACA‑TCGC‑TAAAGAAATTCCCCCGGCAAC‑AAAAATTTTTGAGTTACCGTGCAGGGCAAAACCCCCTTCC‑GGGGACAACGATTTCAGCGGCTCCGGGCCCAAACTGAACTTTTTCAGATTATGCGCCCCGTACCTCTTAAAGCCGCCtttc                                              >  1:263302/1‑249 (MQ=255)                                                                                                                                                                                                                gggtgtcgaaacccccccccacgagggCGG‑TTAGGGTTGGGAGAATCA‑GAAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGcc                                     <  2:553080/226‑1 (MQ=255)                                                                                                                                                                                                                  ggttgggtaacccccACCAGGAGGGCGT‑TAAAGGTGGGGGAAATCA‑G‑AAGCCCCGCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCgg                                   <  1:11294/239‑1 (MQ=255)                                                                                                                                                                                                                                                    gtggCGGGAACCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtggtg  <  1:220394/248‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                GCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACA‑TTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGG‑TTAAGGTTGCGGGAATCA‑GCAAGCCC‑GCCAC‑TTTT‑A‑CCCGATTTTCCGGTGCGGCAGTGAAGTACA‑TCGC‑TAAAGAAATGCCCACGGCACC‑AAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTG  >  NC_000913/3862419‑3862910

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