Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,667	0	G	A	37.7%	67.5 / 24.8	61	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 (27/11);  new base A (0/23);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.21e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGT  >  NC_000913/3862430‑3862897                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑AGGTGCgggg                                                                                                                                                                                                                                                >  1:1085342/1‑238 (MQ=255) tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGGTGCgggg                                                                                                                                                                                                                                                >  1:3493366/1‑238 (MQ=255)      aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTA‑AA‑GGTGGGGGGaaca                                                                                                                                                                                                                                             >  1:1975488/1‑236 (MQ=255)      aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTG‑AA‑GGGTGCggggagta                                                                                                                                                                                                                                            >  1:2886037/1‑233 (MQ=255)            gTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGAT‑AA‑GGGTTCGGGGAT‑caa‑c‑aaag                                                                                                                                                                                                                                   >  1:3426405/1‑236 (MQ=255)                 catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑TA‑AGTTGGGGGGAA‑CAA‑C‑AAACCCGcc                                                                                                                                                                                                                              >  1:1924177/1‑238 (MQ=255)                            cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTTCGGGGAA‑CAG‑C‑AATCCCCCCACATTTTCCCGa                                                                                                                                                                                                                  >  1:800336/1‑239 (MQ=255)                                       aaCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTGCGGGGAA‑CAA‑CAAAACCCCCCCTTTTTCCCCGCTTTTCCgtgg                                                                                                                                                                                                        >  1:3400569/1‑236 (MQ=255)                                                      ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAA‑GGTGGGGGGACA‑AAC‑A‑AAGCCCCCCCCTTTTCCCCATTTTTCCCGTGGGGCAG‑GGAGTAC‑a                                                                                                                                                                                        >  1:1114516/1‑238 (MQ=255)                                                      tgtggggacaaaaCCTAATAAACCAGCCCTCCGCGCGTGCTCCCCCCCACAAAAAAGAAATACACTACGCGCGAAAGGGAGCCCCCCAACACGGCGTGCACCGAAAACGGCAGCGGTAAAT‑GGAACAGGAGGTTAAATTCAAGAGTTTCGGTAATCCCCACTGGGGGGGT‑GTTT‑AA‑GGTGGGGGGAAT‑AAC‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑At                                                                                                                                                                                       <  2:2468702/230‑1 (MQ=255)                                                         tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GAGTTTGGGAAG‑AAG‑C‑ACACCCCCCACTTTTAACCCCTTTTTCGGTGCGGCAGAGAAGGAC‑Atag                                                                                                                                                                                     >  2:1769382/1‑236 (MQ=255)                                                            ggTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTT‑AA‑GG‑TGCGGGGAA‑CCG‑C‑AAACCCCCCACcttttc                                                                                                                                                                                                                      >  1:3425762/1‑201 (MQ=255)                                                                       gTCAATCACCCCCCCGCCTATGTCGCCCCCCCCACAAAAGATATACATCCGGCTCAACATTGGGCCCGCCAGCACGGCGTGTCCCAAAAACAGAAGCGGAGAAA‑TGAACGGGAAGGTAAATTCCGGGGGTTCGTTAACCCCCCCCAGCATGGC‑GGTA‑AA‑GGTGGCGGAAAA‑CAC‑C‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCa                                                                                                                                                                     <  2:500802/239‑1 (MQ=255)                                                                                                     cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGTTTCGGGAAT‑CAG‑A‑CAACCCGCCCCCTTTATCCGCTTTTCCCGTGCGGGCGGGAAGGAA‑AACGCTTAAAAAAAG‑CCCCCCCCACCAAAGC‑ATTTTTAG‑TTA‑CCCTGCa                                                                                                                                    >  2:1174966/1‑239 (MQ=255)                                                                                                       caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGGTGGGGGGAT‑CAG‑C‑AAAACCCCCCCCTTTACCCGATTTTTCCGAGCGGCAGGGAAGTACAATGGCCAAAAAAATC‑CCCACGCCCCACAAAA‑TTTTTAGT‑TTA‑CCG‑GCAgg                                                                                                                                  >  2:1227261/1‑239 (MQ=255)                                                                                                                  tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCCAACAGCAACGGCGAAA‑TGAACAGGAAGGCAAATTCCAGCGGTTCTGTAATTCCCACCAGCATGGC‑GGTT‑AA‑TGTTTCGGGGAT‑CTGAC‑AACCCGCCCCCTTTTACCC‑ATTTTCCGGGGCGGGCGGGAAGGAA‑CTCACTCAAAAAAAA‑ACCCCCGCGCCCATAA‑CTTTTTAG‑TTATCCGTGCTGGGCaaaagacg                                                                                                                        >  2:1054440/1‑235 (MQ=255)                                                                                                                    cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGGTGCGGGAAA‑CAA‑C‑AAGACCGCCCCCTTTACCCCGTTTTTCCGTGCGGCCGGGAAGAAC‑AACGCGAAAGAAAAA‑CCCACGGCACCAAAAA‑ATTTTTAGTTTA‑CCG‑GCACGGCAAAAcaaccc                                                                                                                      >  2:1794272/1‑233 (MQ=255)                                                                                                                            tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGTTTCGGGGAA‑CAG‑C‑AAAACCCCCCCCTTTAACCGCTTTTTCGGGGCGGCAGTGGAgga                                                                                                                                                                                           >  1:2519771/1‑164 (MQ=255)                                                                                                                               aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGTTGGGGGGAT‑TAG‑C‑CAGCCCGGCACtttttc                                                                                                                                                                                                                      >  1:2443745/1‑134 (MQ=255)                                                                                                                               aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTT‑AA‑GG‑TGCGGGGAA‑CAA‑CAAAACCCCCCCCTTTTACCCGGTTTTCCGG‑GCGGGCGGGAAGGAA‑ACACCCAAAAAAAAG‑CCCCCCGCCCCAAAAATTTTTTGGG‑TTC‑CCGCG‑AGGGCAAAAACG‑CCTTCCGGGgaccac                                                                                                         >  2:1757724/1‑239 (MQ=255)                                                                                                                                     aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGTTGCGGGGAT‑CAA‑C‑CAACCCCCCCCCTTTTCCCCCTTTTTCGGTGGGGCCGGGAAAGAA‑AACGCTAAAAAAAAA‑CCCCCCGCCCCCAAAA‑TTTTTTAA‑GTT‑TCGTGGCGGGCAAAAACG‑CCCTCCTGGgacaaccaattt                                                                                                   >  2:2846435/1‑233 (MQ=255)                                                                                                                                                   gtgtctACCGCAAACACCCAGGGTAAAA‑TAAAAGGGAAGGTAAATCCCGGGGTTTCGGTAATTCCCACCGCCAGGGC‑GGTT‑GG‑GGTGGGGGGAAT‑CAG‑A‑AAACCCGCCACTATTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCAc                                                                                                                                                                    <  2:368282/159‑1 (MQ=255)                                                                                                                                                         aCCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑TA‑GGGTGCGGGGAA‑CAA‑CAAACCCCCCCCCTTTTACCCCTTTTTCCGGGGGGGCGGGGAGAAAA‑ACCCCAAAAAAAATC‑CCCCCGCACCCAAAAA‑TTTTTGGG‑TTA‑‑CGCGCAGGGGAAAAACC‑CCCTCCCGGGACAAAGCTTTCTCCGGCTCGgcgccccaaa                                                                                >  2:1133718/1‑232 (MQ=255)                                                                                                                                                            gaaaaCAGCAA‑GGTGAAAGTGAAAGGGAGGGTAACTCCCGGGGGTGCGGAAACCCCCCCAGGCTGGGG‑GGTT‑AG‑GGGTGGGGGAAT‑CAC‑C‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                            <  2:1089777/237‑1 (MQ=255)                                                                                                                                                             cAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGGTGCGGGAAT‑TAG‑C‑CAACCCGCCACTTTTACCCGATTTTCCCGTGCGGCAGTGAAAGAA‑ATCGCTAAAGAAAAGCCCCCCGGCCCCaaa                                                                                                                                                           >  1:368282/1‑164 (MQ=255)                                                                                                                                                                aCAGAAACGGTAAAA‑TGAACAGGAAGGAAAATCCCAGGGGTTCGGTAATCCCCACCAGCAGGGG‑GTTT‑AG‑GGTTGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                        <  1:2331164/117‑1 (MQ=255)                                                                                                                                                                  aaaaaCGGTAAAA‑TGAACGGGAAGTTAATTCCCGGCGGTTCGGTATTCCCCACCAGCTGGGC‑GTTT‑AG‑GGTTGGGAAAAA‑AAC‑A‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCg                                                                                               <  1:437563/211‑1 (MQ=255)                                                                                                                                                                     aaCGGGGAAA‑TGAACAGGAAGGAAATTCCCGGCGTTTCGTTAATCCCCACCAGCATGGC‑GGTT‑AG‑GGTGGGGGGAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCg                                                                                          <  1:2453067/216‑1 (MQ=255)                                                                                                                                                                         ggaaaa‑aGAACAGAAAGTAAAATCCCAGGGGTCCGGTAATCCCCACCACCGGGGG‑GGTA‑AG‑GGTGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                <  1:1047547/235‑1 (MQ=255)                                                                                                                                                                         gTGAAA‑TGAACAGGAGAGTAAATGCCAGGGTTTCGGTAACCCCCACGGGCAGGGG‑GTTT‑AA‑GGTGGGGGAAAT‑CAG‑C‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                <  2:279067/238‑1 (MQ=255)                                                                                                                                                                          tGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTTCGGGGAT‑CAG‑C‑AAACCCCCCACCTTTTCCCCATTTTTCCGTGCGGGAGGGaa                                                                                                                                                                                              >  2:2331164/1‑117 (MQ=255)                                                                                                                                                                            aaa‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑AGTTGCGGGGAT‑CAG‑C‑AAGCCCCCCCCTTTTTCCCCATTTTTCGGTGCGGCAGTGAAGTAC‑ATCGCTTAAGAAATG‑CCCACGGCACCAAAGA‑tttttta                                                                                                                                                 >  1:640393/1‑155 (MQ=255)                                                                                                                                                                            aaa‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑TA‑GGGTGGGGGGAT‑CAg                                                                                                                                                                                                                                          >  2:2384981/1‑73 (MQ=255)                                                                                                                                                                            aaa‑TGAACAGGAAGGAAAATCCCGGGGTTTGGTTAATCCCCCCCAGCAGGGC‑GTTT‑AG‑GGTGGGGGAAAC‑CGG‑C‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                <  1:2194497/235‑1 (MQ=255)                                                                                                                                                                             aa‑TGAACAGGAAGGAAATTCCCGGGGGTTCGGTATTCCCCCCCAGCATGGG‑GGTT‑AG‑GGTGGGGGAAAT‑CAG‑C‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                           <  1:153194/239‑1 (MQ=255)                                                                                                                                                                                tGAAAAGGAAGGAAAATTCCGGCGTTTCGTTAATCCCCCCCAGCAGGGC‑GGTT‑AG‑GGTGGGGGAAAT‑AAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                         <  1:1133718/239‑1 (MQ=255)                                                                                                                                                                                 aaaagggAAGTTAAATCCCAGCGGTTCGTAAATTCCCCCCATGAGGGC‑GGTT‑AA‑GGTTGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGccc                                                        <  1:2376252/234‑1 (MQ=255)                                                                                                                                                                                  aaCAGGAAGTTAAATTCCGGCGGTTGGGTAATCCCCCCCAGCAGGGC‑GTTT‑AG‑GGTGGGGGAAAC‑CAC‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGAt                                                                                                      <  1:1279858/192‑1 (MQ=255)                                                                                                                                                                                         aGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGa                                                                                                                                                                                                                                                <  1:1843742/56‑1 (MQ=255)                                                                                                                                                                                          ggaaaaTCCCAGGGTTTCGTTAATCCCCCCCAGCAGGGG‑GGTT‑AA‑GGTGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGAttt                                                                                                                                                                                                               <  1:3561767/82‑1 (MQ=255)                                                                                                                                                                                           gaaaaTCCCGGCGTTTCGTTATTCCCCCCCAGCTTGGC‑GTTT‑AA‑GGTTGCGGAAAA‑CAG‑A‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                         <  1:1769382/226‑1 (MQ=255)                                                                                                                                                                                            tAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGTTGCGGGGGA‑CAG‑G‑AAGGCCGCCACTTTTTCCCCATTTTTCGGTGCGGCAGGGGAGGAC‑ATCCCCAAAGAAATG‑CCCACGGCACCCAAGA‑ATTTTTAG‑TTA‑CCGTGCAGGGCAAAACCGCCCTTCCGGGAACAA‑GATTTTAGCgggc                                                                                            >  2:1279858/1‑190 (MQ=255)                                                                                                                                                                                            aaaaCTCCGGCGGTTGAGTAATCCCCGCAAGCGTGGG‑GTTT‑GA‑GTTGGGGGAAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                         <  2:1975488/226‑1 (MQ=255)                                                                                                                                                                                                         ttCGGTAATCCCCACCAGCGTGGC‑GGTT‑AA‑GGTGGCGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                          <  1:566194/197‑1 (MQ=255)                                                                                                                                                                                                              aactcccccccAGCATGGC‑GGTT‑AA‑GGTTGCGGAAAT‑CAG‑C‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAg                                                                                                                                                <  2:2023933/120‑1 (MQ=255)                                                                                                                                                                                                                atccccACCAGCATGGC‑GGTT‑AA‑GGTGGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                         <  1:1174966/236‑1 (MQ=255)                                                                                                                                                                                                                 tcccccccGGCAGGGC‑GTTT‑AA‑GGTGGGGGGAAC‑CAG‑A‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAAt                                                                                                                                                                           <  1:1717048/95‑1 (MQ=255)                                                                                                                                                                                                                     caccaGCGGGGC‑GGTA‑AA‑TGTGGCGGAAAT‑CAG‑C‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                    <  1:1237659/239‑1 (MQ=255)                                                                                                                                                                                                                     caccaGCAGGGG‑GGTT‑AA‑GGGGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                    <  1:884707/239‑1 (MQ=255)                                                                                                                                                                                                                      cccaGCTGGGC‑GGTT‑AA‑GTTTGGGGAAAA‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                           <  1:2738076/198‑1 (MQ=255)                                                                                                                                                                                                                      cccaGCA‑GGG‑GGTT‑AAGGGTTGCGGAAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                     <  1:138546/172‑1 (MQ=255)                                                                                                                                                                                                                       ccaGCTGGGC‑GGTT‑AG‑GTTGGCGGGAAT‑AAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                  <  1:3468314/239‑1 (MQ=255)                                                                                                                                                                                                                       ccaGAAGGGC‑GGTG‑AG‑GGTGGCGGAAAC‑CAA‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                  <  1:2199944/239‑1 (MQ=255)                                                                                                                                                                                                                        caGCATGGC‑GGTT‑AA‑GGGTGCGGGGAT‑CAG‑C‑AAACCCCCCCCTTTTTCCCGCTTTTCCCGTGCGGCCGGGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTTAG‑GTT‑CCGCGCAGGGCAAAACCC‑CCTTCCGGGGACAACGATTTTTGCCGCTCCGCCGTCCACCCGCACTCTTGCTGAAGACGCGCCCAGTACCTATCACTGCCCCCCTCTACAGCTccccggcaca                 >  2:2199944/1‑232 (MQ=255)                                                                                                                                                                                                                           cATGGC‑GGTT‑AA‑GGTTGCGGGGAT‑CAG‑C‑AAGCCCCCCCCTTTTTCCCGATTTTCCGGGGCGGGAGGGGAGGAC‑a                                                                                                                                                                                        >  1:3029603/1‑74 (MQ=255)                                                                                                                                                                                                                           cATGGC‑GGTT‑AA‑GGTGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCgg                                                                                                                                                                                                          <  1:448399/57‑1 (MQ=255)                                                                                                                                                                                                                                c‑GGTT‑TA‑GGGTGCGGGAAT‑CAG‑C‑AAGCCCGCCCCTTTTTCCCCATTTTTCCGGGCGGCAGTGGAGGAC‑ATCGCTTAAGAAATG‑CCCCCGGCACCCAAAA‑TTTTTTAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTTAGCGGCTCCGCGGTCCAAACGAAc                                                                           >  2:138546/1‑173 (MQ=255)                                                                                                                                                                                                                                 ggttt‑AA‑GGTTGCGGAAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa          <  1:2142843/235‑1 (MQ=255)                                                                                                                                                                                                                                         gagtggCGGGAAACCAG‑C‑AACACCGACACTTTTACCCGATTTTCCGGTGCGACAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa          <  2:2922638/228‑1 (MQ=255)                                                                                                                                                                                                                                           gtggCGGAAAC‑CAG‑C‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt  <  1:3410793/235‑1 (MQ=255)                                                                                                                                                                                                                                            ttGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCa                                                                                                                                                             >  1:2078901/1‑86 (MQ=255)                                                                                                                                                                                                                                            tggCGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt  <  1:2089341/235‑1 (MQ=255)                                                                                                                                                                                                                                               cGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                     >  2:2274838/1‑199 (MQ=255)                                                                                                                                                                                                                                               cGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                     <  1:2274838/199‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGT  >  NC_000913/3862430‑3862897

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