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.81.3% 253.7 / 40.2 80pseudogene (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 (1/14);  new base . (44/21);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.62e-05
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.80.1% 239.1 / 36.1 80pseudogene (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/14);  new base . (44/20);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.36e-06
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.

CAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCC  >  NC_000913/2173134‑2173573
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                   
cAGTTGGATGTAAAAACAAAT‑GTTTGCCGTGCGATTCCCAGGGGTG‑AGCCCAATATTGATGACCCCAGAGATTAGGGGGGGAAACAGATTCCCCGGAGGCCCGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATaccac                                                                                                                                                                                                                      <  2:1036528/228‑1 (MQ=255)
          tgCAAACCAATAGTTTGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAAA‑TCAGGCTTGCCGAAacaaccgccgaagga                                                                                                                                                                                                            >  1:937192/1‑219 (MQ=255)
            cAAACCAATAGTTTGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACAA‑‑CCGCCGTATGtg                                                                                                                                                                                                         >  2:1557199/1‑226 (MQ=255)
            cAAACCAATAGTTTGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAA‑TCAGGCTTGCCGATaccaccgcccaaagtc                                                                                                                                                                                                           >  1:739522/1‑217 (MQ=255)
             aaccaaaTAGTTGGGGTCGGAATCCCCG‑GGTTG‑ATGCCCAAAATGATGACACCCGAGATAAGGGTGCGGACCAGTTTCCCAGGATGCCCGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                                                                                                       <  2:734197/223‑1 (MQ=255)
             aaaCCAAGTGTTCGCGTCGCGACCCA‑A‑GGGTGCAGTCTAATAATGATGACCCCCGAGATAAAGGCGCGGAACAGATTCCCACGATGCCCGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                                                                                                       <  2:10660/228‑1 (MQ=255)
                         tGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                           >  1:651411/1‑228 (MQ=255)
                         tGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATACAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCACGCAGCAACTcc                                                                                                                                                                                           >  1:504053/1‑228 (MQ=255)
                           cGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                         >  1:1330837/1‑228 (MQ=255)
                           cGTCACGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                         >  1:779673/1‑228 (MQ=255)
                            gTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                         >  1:1162616/1‑227 (MQ=255)
                             tCGAGATCCAAA‑GGGTG‑ATCCTCATAATGATGACCCCGGAGATTAGGGTGCGGAACAGTTTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                       <  2:189626/228‑1 (MQ=255)
                                   tcacca‑GGGTG‑ATGCTCATAATGATGACACCGGAGATTAAGGTGCGGAAAAAATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                 <  1:149244/224‑1 (MQ=255)
                                           ggTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCGAGGCCAATCa                                                                                                                                                                           >  1:1295649/1‑227 (MQ=255)
                                             tG‑ATCCTAATAATGATGCCCCCCGAGATTAAGGTGCGGAACAGATTCCCCCGATGCCCGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAA‑CCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                         <  1:669938/227‑1 (MQ=255)
                                                 tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCACAC‑‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                     >  1:1499576/1‑228 (MQ=255)
                                                 tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                  >  2:301204/1‑199 (MQ=255)
                                                 tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                     >  1:744353/1‑228 (MQ=255)
                                                 tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAATGCCAATCAGGAAc                                                                                                                                                                      >  1:787626/1‑227 (MQ=255)
                                                  gCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                                                                                                 <  1:301718/199‑1 (MQ=255)
                                                   cTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAGGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                    >  1:46863/1‑227 (MQ=255)
                                                      aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACACCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                 >  1:729972/1‑227 (MQ=255)
                                                       aaagctagccccccGATATTAAGGGCGGGAACGGATTTCCACGATGCAGGGCG‑ACCGCCATCGCCACGAAGAA‑CCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                  <  2:581226/214‑1 (MQ=255)
                                                           gacgccacGCGAGATAAGGGTGCGAAACAGATTCCCACGATGCACGGCG‑ACCGCCACCGCCAGGAAGAA‑GCCGAGGGTGGAAAGACCGCCAAACGGAACCCCCTGATTCCCCGGCACACAACCAGCAATTAAAATGGTGAGGGGGATAAAAA‑CCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                            <  2:743093/222‑1 (MQ=255)
                                                                  cccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                  >  1:1800005/1‑198 (MQ=255)
                                                                      agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                 >  1:161609/1‑211 (MQ=255)
                                                                      aaaTTAGGGTCCGAAACAGTTTCCCACGATGCCCGGCG‑CCCGCAAGCGCCACGAAGAA‑GCCGATGGTGGAAAGATCCCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                  <  2:365008/224‑1 (MQ=255)
                                                                      aaaTTAAGTTGCGGAACAGATTCCCAGGAGCCACGCCA‑CCCCCCACCCCCCGAAAAAA‑CCCGATGGGGGCAAGACCGCAAAACGGCAGCACCTAATTCCCCGGCACACAAACACCAATAAAAATGGTGAGTGGGAAAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                                                                                <  2:738267/226‑1 (MQ=255)
                                                                                 cGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACgccca                                                                                                                                     >  1:698393/1‑224 (MQ=255)
                                                                                   gaaacggtttCCCAGATGCCCGGCGCGCCCCCATCGCCACGAAGAA‑CCCAATGGGGCAAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAA‑TCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCGAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                                                                                     <  2:1293333/220‑1 (MQ=255)
                                                                                      cAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGACCTTCGCCTGTAAACGACTACgcac                                                                                                                                >  1:1403210/1‑225 (MQ=255)
                                                                                      aagttcaccACGTCGCCCGGAG‑ACCCCCATCGCCACGAAGAC‑GCCGATGGGCGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAAGGGTGAGTGGAATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                <  2:744888/221‑1 (MQ=255)
                                                                                        gATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                              >  1:472253/1‑228 (MQ=255)
                                                                                                   agccagggacaCCGCCAGCCCCA‑GAAGAACGCCGAGGGTAGAAAGACCACCAAGCGGCAGCCCCTAATTACCCGGCACACAAACAGCAATTAAAAGGGGGAGGGGAAAAAAAA‑CCGGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                    <  2:226513/218‑1 (MQ=255)
                                                                                                              ccGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATaccaccgcc                                                                                                                                                                                                                  >  1:1752090/1‑122 (MQ=255)
                                                                                                              ccGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                                                                                         >  1:1314962/1‑228 (MQ=255)
                                                                                                               cGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATaccaccgccg                                                                                                                                                                                                                 <  2:1748997/125‑5 (MQ=255)
                                                                                                                gCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                                                                                                       >  2:454971/1‑228 (MQ=255)
                                                                                                                       ccAGGAAGAC‑ACCGATGGCGCCAAGATCCACAAACGGCATCAGCTAATACCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAC‑TCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                <  2:1507290/228‑1 (MQ=255)
                                                                                                                       ccACAAAGAA‑GCAGATGGTGGAAAGATCGCAAAACGGCGACACCTGATACCCCGGCAAACAAACACAAATTAAAATGGTGAGTGGAAAAAAAT‑TCAGCCTT‑CCCAAACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                <  2:512535/228‑1 (MQ=255)
                                                                                                                         aCGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                 >  1:700519/1‑193 (MQ=255)
                                                                                                                           gaagac‑GCAGATGGTGGAAAGATCCCCAACCGGACGCCCCTAATTCCCCGGCACACAAACAGAAATTAAAAGGGTGGGTGGAAAAAAAC‑TCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                            <  2:1264123/228‑1 (MQ=255)
                                                                                                                                  gCCGATGGTGGAAAGATCGCCAACGGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGGTGAGGGGAAAAAAAA‑TCAGGCTGGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                      <  2:697222/228‑1 (MQ=255)
                                                                                                                                       tggtggCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                 <  1:201355/228‑1 (MQ=255)
                                                                                                                                         gtgcaaaGGTCGCAAACCGGCAGCACCTGATTCCCAGCAACAAAAACACCAATTAAAAGGGTGAGTGGAATAAAAA‑TCAGGCTGGCCAAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                <  2:371359/222‑1 (MQ=255)
                                                                                                                                         gaaacaagtccGCCAATCGGCAGCACATGATTACCCGCCACAAAAACAGAAATAAAAAGGGTGAGGGGAATAAAAC‑TCAGCTTTGCCGACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                <  2:503234/217‑1 (MQ=255)
                                                                                                                                                  gccacacaCGGCAGCCCCAGATTCCCCGCCA‑AAAAACCACAAATAAAATGGGGAGTGGGAAAAAAACTCGGCCTTCCCG‑‑ACCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                              <  2:46784/213‑1 (MQ=255)
                                                                                                                                                     ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATGAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTtgcagc                                                                   >  2:1832651/1‑228 (MQ=255)
                                                                                                                                                                     gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                                                                                                           >  1:1580559/1‑108 (MQ=255)
                                                                                                                                                                     gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTGGCAGCTTCACTGCCGCCGgtt                                                   >  2:1561818/1‑228 (MQ=255)
                                                                                                                                                                      aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                          <  2:1577756/108‑1 (MQ=255)
                                                                                                                                                                          cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                              >  1:455093/1‑228 (MQ=255)
                                                                                                                                                                           ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                                                                                                     >  1:2035260/1‑172 (MQ=255)
                                                                                                                                                                             ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                                                                                                           >  1:977060/1‑132 (MQ=255)
                                                                                                                                                                              gCACACAAACAGCAATAAAAATGGTGAGGGGGATAAAAA‑CCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                          <  2:975312/132‑1 (MQ=255)
                                                                                                                                                                                 ccaaaacacgcaTTAAAATGGGGAGGGGAAAAAAAA‑GCCGTCTCGCC‑CGACCCCACCGCGGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                        <  2:1328458/217‑1 (MQ=255)
                                                                                                                                                                                      acggaagtaaaatTGGTGAGTGGGAT‑AAAA‑TCCGGCCTGCC‑CCATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                  <  1:482524/219‑1 (MQ=255)
                                                                                                                                                                                                  aTGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACa                      >  1:1212425/1‑228 (MQ=255)
                                                                                                                                                                                                    ggTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                                                                                                     >  1:970869/1‑131 (MQ=255)
                                                                                                                                                                                                    ggTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGACGTTTTTACCGCCAGCTGCAATACACCCTTGACATc                    >  2:1355835/1‑228 (MQ=255)
                                                                                                                                                                                                    gaccagggggAAAAAAA‑ACCGGCCCTCCCATACC‑CCCCGCCGTAGGTCCCCGCAGCACCCCCGGATCAGGCCCAATCGGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                         <  2:205044/217‑1 (MQ=255)
                                                                                                                                                                                                     gggggTGGAAAAAAAA‑TAAGCCTTCCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                    <  2:969135/127‑1 (MQ=255)
                                                                                                                                                                                                     gTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                 >  2:1065218/1‑118 (MQ=255)
                                                                                                                                                                                                      tGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                <  1:1067132/118‑1 (MQ=255)
                                                                                                                                                                                                       gAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTGTGTCCCAGCAGCAACGCCGGATCaa                                                                                                                                                                                    >  1:1885599/1‑65 (MQ=255)
                                                                                                                                                                                                       gAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                 >  2:440661/1‑228 (MQ=255)
                                                                                                                                                                                                       gAGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                 >  1:1309969/1‑228 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTGTGTCCCAGCAGCAACGCCGGATCAAg                                                                                                                                                                                   <  2:1882262/65‑1 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                              >  1:290935/1‑150 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                  >  2:1261713/1‑162 (MQ=255)
                                                                                                                                                                                                         gTGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                 <  1:1263969/162‑1 (MQ=255)
                                                                                                                                                                                                          tGGGATAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCCTTGACATCGTAAc               >  1:124997/1‑227 (MQ=255)
                                                                                                                                                                                                               tAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                      >  2:1004703/1‑103 (MQ=255)
                                                                                                                                                                                                               tAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                >  2:1118064/1‑141 (MQ=255)
                                                                                                                                                                                                               tAAAAA‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATCAACCTTTGACATCGTAACCcgcac         >  1:1800918/1‑226 (MQ=255)
                                                                                                                                                                                                                aaaaa‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                                                                                     <  1:1006529/103‑1 (MQ=255)
                                                                                                                                                                                                                aaaaa‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                               <  1:1120073/141‑1 (MQ=255)
                                                                                                                                                                                                                    a‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                      >  1:751522/1‑98 (MQ=255)
                                                                                                                                                                                                                    a‑TCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                          >  1:205396/1‑206 (MQ=255)
                                                                                                                                                                                                                      tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                                                                                     <  2:750252/98‑1 (MQ=255)
                                                                                                                                                                                                                        ggccTTGCCGAT‑CCA‑CCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc  <  1:890221/224‑1 (MQ=255)
                                                                                                                                                                                                                        aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                >  2:1850199/1‑53 (MQ=255)
                                                                                                                                                                                                                         ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                               <  1:1853481/53‑1 (MQ=255)
                                                                                                                                                                                                                              gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                      >  2:804396/1‑89 (MQ=255)
                                                                                                                                                                                                                               ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                          >  1:1371800/1‑196 (MQ=255)
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                   
CAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACA‑GGGTG‑ATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAA‑GCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCC  >  NC_000913/2173134‑2173573

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