Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% pseudogene (1‑2/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.92.6% 262.6 / 7.0 67pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (40/22);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 2.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.89.6% 247.9 / 14.4 67pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base . (40/20);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 1.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAC  >  NC_000913/2173123‑2173591
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                     
acATTAGCCGCCAGTTGGGTGTGAAAACCAATAGTTGGCGTAGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGAGGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                                                                                                                        <  1:919093/238‑1 (MQ=255)
  aTTAGCCGCCGGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑CCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATcccac                                                                                                                                                                                                                                        <  1:374510/237‑1 (MQ=255)
      gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTACGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                                                                                                                  >  1:743154/1‑234 (MQ=255)
           catgtGGGTGTGCAAACCAATAGTTTGCGTCGCGCTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTACCACAATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAATT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                                                                                                                <  1:880787/232‑8 (MQ=255)
           cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGCTTTCCACGATGCACGGCG‑CCCGCCATCGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAACACCTGATTACCAGGCACACAAACAGCAATTAAAAT‑GGTGCGTGGGCTAAAAATCAGGCTTGCCGATaacaccgccctatgt                                                                                                                                                                                                                              >  2:1084551/1‑229 (MQ=255)
           cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                            >  2:675130/1‑238 (MQ=255)
           cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                            >  1:30443/1‑238 (MQ=255)
                   tgtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                                                                                                                   >  1:1201744/1‑239 (MQ=255)
                        aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                              <  1:308107/239‑1 (MQ=255)
                        aaaCCAATAGTTGGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAAATTAAGGTGCGGACCAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                              <  1:108066/239‑1 (MQ=255)
                                   ttGTGTCGCGATCCACAGGGTGATGCTCAAAATGATGCCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                    <  1:727065/238‑1 (MQ=255)
                                    tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                  >  1:962536/1‑239 (MQ=255)
                                    tGCGTCGCGATCCACAGGGAGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAATGCCGAATCAAggac                                                                                                                                                                                                  >  1:1045904/1‑237 (MQ=255)
                                      cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCaa                                                                                                                                                                                                >  2:697520/1‑239 (MQ=255)
                                      cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGAGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                >  2:1088605/1‑239 (MQ=255)
                                                         gtgcgaagaaggattaCGCCCGCGATGAAGGTGGGAAACAGATTCCCACGATCCAGGGCC‑ACCGCCATCCACAGGAAGAAGCCGATGGTGCCAAGATCGACAACCGGCAGCACCGGATTCCCCGGCACACAAACAGCAATTAAAGT‑GGGGAGTGTGATAAAAATCAGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCGGAccgc                                                                                                                                                                              <  2:777406/227‑1 (MQ=255)
                                                         aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAACCCGGATCAAGGCCACTCAGGAACCCCTGAccccc                                                                                                                                                                             >  2:967194/1‑239 (MQ=255)
                                                               aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                                        >  1:763899/1‑238 (MQ=255)
                                                                aactaaggaCCCCGGAGATTAGGAGGGGAAACAGATGCCCCGGATCCACGGCC‑ACGGCAACCGCCACGAAAAGGCCGATGGTGGCAACACCGCCAAAGGGCAGCACCTGATTCCCCGGCACACAAACAGCAATAAAAAG‑GGTGAGTGGAATAAAAATCAGGCTGGCGGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                      <  2:922289/232‑1 (MQ=255)
                                                                 aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTc                                                                                                                                                                      >  1:1045702/1‑238 (MQ=255)
                                                                      tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                >  1:707449/1‑239 (MQ=255)
                                                                         cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGCAGAAGCCGATGGTGGCACGATCGCCAAACGGCATCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGCATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCGGCAGCAACGCCGGATAAAGGCCAATCAGCCACTCCTGACCCCCGAACTTCGCCTGTaa                                                                                                                                                              >  2:1157273/1‑238 (MQ=255)
                                                                                    agggTGCGGAAGATACGTCCGCCATGCACGGCA‑CCCGCCGTCGCCCCGAAGAAGCCGGTGGGGGAAAAACCTCCAAACGGCACCCGTGGATTACCCGGCACACAACCACCAATAAAAAG‑GGTGAGTGGGATAAAAACCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:645526/237‑1 (MQ=255)
                                                                                    aagtgggggACCATATTCCCACGAGCCCCGCCA‑ACCCCCACCGCCACGAAGAACCCGAGGGTGGCAAGATCGCAAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAAT‑GGTGAGTGGAATAAAATCCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:468671/232‑1 (MQ=255)
                                                                                    aaGGGGCGGACCAGATTTCACCGTGGCGCGCGG‑ACCGCCACCGCCACAAAGAAGCAGGTGGTGGAAAAAGCGCAAACCGGCAGCACGGGATTACCCGACACAAAAACAGAAATTAAAAT‑GGGGAGGGGAAAAAAAATCAGGTTGGCCGAT‑CC‑CCCCGCCGTATGCCCCAGCAGCAACCCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:864333/239‑1 (MQ=255)
                                                                                      ggTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                                   >  1:274562/1‑236 (MQ=255)
                                                                                      ggTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                                   >  1:926590/1‑236 (MQ=255)
                                                                                          cGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                            >  1:574060/1‑239 (MQ=255)
                                                                                          cGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                            >  1:206302/1‑239 (MQ=255)
                                                                                          cGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                            >  1:821240/1‑239 (MQ=255)
                                                                                           gaaaCAGATTTACACGATGCACGGCG‑ACCGCCATCGCGACGAAGAAGCCGATGGCGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCCTGCCGATACCA‑‑CTGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCAAACTTCGCCTGTAAAAGACAACAACCCAGCtca                                                                                                                                           >  2:189241/3‑237 (MQ=255)
                                                                                                agtttcccACGATCCACCGCG‑CCCGCCAGCCCCACAAAGAAGCCGATGGTGAAAAGAATCCCAAACGGCAGCACCTGATTCCCCGGCCCACAAACCGCAATTAAAAT‑GGGGAGTGGAATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                      <  2:962536/233‑1 (MQ=255)
                                                                                                 ctattgccggATGCACCGCGCACCGCCACCGCCAAGAAGAAGCCGAGGGGGGCAAACGCGCCAAATGGCGGCCCCTGTTCCCCCGGCCCAAAAACAGCAATTAAAAT‑GGGGAGTGGAATAAAAATCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                      <  2:743154/230‑1 (MQ=255)
                                                                                                  aTTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCACCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                    >  1:957885/1‑239 (MQ=255)
                                                                                                   tttCCACGATGCACGGCG‑CCCGCCATCGCCACGAAGAAGCCGATGGTGGAAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAT‑GGTGGGTGGGATAAAAACCGGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                      <  1:303917/236‑1 (MQ=255)
                                                                                                       cACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                               >  1:964414/1‑239 (MQ=255)
                                                                                                         cGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                                                                                                             >  1:969060/1‑239 (MQ=255)
                                                                                                          gATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                            <  1:1088605/239‑1 (MQ=255)
                                                                                                          gATGCACGGCG‑ACCGCCACCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                            <  1:192360/239‑1 (MQ=255)
                                                                                                            tGCACGGCG‑CCCGCCATCGCCAGGAAGAAGCCGATGGTGGCAAGAACGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGAGAGAGGGATAAAAATCAGGCTTCCCGATACCA‑‑CTGCCATATGTCCCAGCAGCAACGCCAGATCAAGTCCAATCAGGAACTCCTGACCGCAGAACTTCGCCTGTAACCGACTCACAACCTGCTAAGCGACGGGCGTCAAGcc                                                                                                                          >  2:433049/1‑239 (MQ=255)
                                                                                                            tGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATAAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                          >  2:767922/1‑239 (MQ=255)
                                                                                                            tGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                          >  2:858007/1‑239 (MQ=255)
                                                                                                                    g‑ACCGCCATCGCCCCGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                   <  1:723762/238‑1 (MQ=255)
                                                                                                                     caCCCCCACCGCACCGAAGAACCCGATGTTGGCAAAATCGCAAAACGGCAGCCCCTAATTCCCCGCCACACAAACAGCAATTAAAAT‑GGTGAGTGGGAAAAAAACCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                  <  2:926590/238‑1 (MQ=255)
                                                                                                                               gCCACGAAGAAGCGGATGGTGACAAGACCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAAG‑GGTGAGTGGAATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                        <  1:1151755/239‑1 (MQ=255)
                                                                                                                               gCCACAAAGAAGCCGATGGTGAAAAGACCGCAAAACGGCAGCCCCTGATTACCCGGCACCCAACAAGAAATTAAAAT‑GTGGATGGGAAAAAAAACCGGGCTTGCCAAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                        <  2:904749/239‑1 (MQ=255)
                                                                                                                                    gaagaaGCCGATGGTGGAAAGATCGCCAACCGGCAACACCTAATTACCCGGACCCCAAACAGCAATTAAAAG‑GGGGAGTGGAATAAAAATCAGGCTTCCCGATA‑CA‑CCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  1:1035369/239‑1 (MQ=255)
                                                                                                                                     aagaagCCGAGGTTGGCAAGATCGCCAACGGGCAGCCCCTGATTACCCGGCACACAAACAGCAATAAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                            <  2:286451/165‑1 (MQ=255)
                                                                                                                                               tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACTCAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                         >  2:562881/1‑238 (MQ=255)
                                                                                                                                               tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                                                                                                  >  1:286451/1‑165 (MQ=255)
                                                                                                                                               ggggggCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACATAACCAGCAATAAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                         <  1:957268/234‑1 (MQ=255)
                                                                                                                                                     cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                                  >  1:334927/1‑239 (MQ=255)
                                                                                                                                                           cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                            >  2:1182079/1‑239 (MQ=255)
                                                                                                                                                                       gcaCCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACACGTGGCATTACCAGCATCACTGCCGCCGATTTTAcc                                                                >  1:916833/1‑239 (MQ=255)
                                                                                                                                                                          ccTAATTACCCGGCCACCAAAAAGAAAATAAAAT‑GGTGAGTGGAATAAAAATCGGGCTTGCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                             <  1:451856/239‑1 (MQ=255)
                                                                                                                                                                          acgaaTTACCCGCCACA‑AAAAAACAAAGAAAATGGGCGACGGGAAAAAAAACCAGGCTTG‑CGAAAAC‑GCCCGCGGTAGGCCCCAGCAGCAACGCGGGATCAAGGCCAATGAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                             <  2:1201744/235‑1 (MQ=255)
                                                                                                                                                                             gATTACCGGGCACAAAAACAGCAATTAAAAT‑GGTGGGTGGGATAAAAATCAGGCGTGCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                          <  1:967194/239‑1 (MQ=255)
                                                                                                                                                                               ttccccGGAACACAAACAGAAATTAAAAG‑GGGAAGTGGAATAAAAACCAGGTTTCCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                        <  1:675130/236‑1 (MQ=255)
                                                                                                                                                                               ttACCCGGCACACAAACAGAAATAAAAAG‑GGGGAGTGGAAAAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                        <  1:311232/239‑1 (MQ=255)
                                                                                                                                                                                  cccGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCaa                                                     >  2:477353/1‑239 (MQ=255)
                                                                                                                                                                                     ggCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                  >  1:837151/1‑239 (MQ=255)
                                                                                                                                                                                     ggCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAACCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCAGTTTTTACCGCAAGCTGCAATac                                                  >  2:594756/1‑239 (MQ=255)
                                                                                                                                                                                       cacacaAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                 >  1:8023/1‑238 (MQ=255)
                                                                                                                                                                                        acacaAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTACGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAACGCCAGATGCAATacac                                                >  2:38570/1‑238 (MQ=255)
                                                                                                                                                                                                       aaaaT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                 >  1:770280/1‑238 (MQ=255)
                                                                                                                                                                                                                 ggTGGGATAAAAATCAGCCGAGCAGACCCCACCC‑‑CCGTATGTCCCAGCAGCAACCCCGGAACAAGGCCAATCAGGATCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                        <  2:1092465/237‑1 (MQ=255)
                                                                                                                                                                                                                 aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                       >  1:1058254/1‑239 (MQ=255)
                                                                                                                                                                                                                       cataaaaCCAGGGTTGCCGAATACA‑CCCGCCGTTGTTCCAAGCAGCACCGCCGGATCAGGGCCAATAAGGAACTCCTGACTGCCGAACTTCTCCTGGAAACGACTACGAGCCTTCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAACCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACGCCTTTGACATCGTAACCCGCGAGGATGCCg                   <  2:299240/235‑1 (MQ=255)
                                                                                                                                                                                                                             aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCATGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATAACTGCCGCCGTTTTTACCGCCAGCTGCAATACACTTTTGACATCGTACCACGCGCGTATGCCGCTAATCa            >  2:785354/1‑238 (MQ=255)
                                                                                                                                                                                                                                aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcc         >  1:1051923/1‑238 (MQ=255)
                                                                                                                                                                                                                                aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACACGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGCCATCGTAACCCGCGAGGATGCCGATAATCAGcc         >  1:614596/1‑238 (MQ=255)
                                                                                                                                                                                                                                      gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac  >  2:435996/1‑239 (MQ=255)
                                                                                                                                                                                                                                        cccatcccccccccccTGTGTCCCAGCAGCAACGCCGAATAAAGGCAAATCAGGAACTCCTGCCCGCCGAACTTCCCCTGTAACCGACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac  <  2:206302/231‑1 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                     
GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAC  >  NC_000913/2173123‑2173591

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