Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/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.7% 238.7 / 36.2 71pseudogene (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/13);  new base . (42/16);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.17e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.81.7% 235.3 / 29.6 71pseudogene (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/13);  new base . (42/16);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.17e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATG  >  NC_000913/2173135‑2173571
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                 
tgTTGGGTGGACACACCATAGGTTTGCGTCGCGCTCCACGGGTTATTGCTATAAATGGTGACCCCCGAAGTTAGGTGGAGGAACAGTTTCCCCCGAGCCCGGGGGACCCCCA‑TCGCC‑ACGAAGAAGCCGATGGGGGCAAGATCGCCAAAC‑GGCAGCACCTGATTCCCCGGCACACAAA‑CA‑GCAATTAAAAG‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATcccaccg                                                                                                                                                                                                                  <  2:1007702/228‑2 (MQ=255)
  ttGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccgcc                                                                                                                                                                                                                >  1:1945157/1‑226 (MQ=255)
         tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccgccgcatgt                                                                                                                                                                                                          >  1:529488/1‑219 (MQ=255)
                        tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGA‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGAATGTCCCAGCAGCAACGCCg                                                                                                                                                                                        >  1:450054/1‑229 (MQ=255)
                        tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCAAACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGGGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGCATGTCCCAGCAACAACGCCg                                                                                                                                                                                        >  1:215883/1‑229 (MQ=255)
                                 ctCCACGGGGTAATGCTAAAAATATGGCCGCCAGAAATTAGGGGGCGAACCAGTTTCCCACGATGCCCGGCGCCCGCCA‑TCGCC‑ACGAAGAAGCGGATGGGGGCAAGATCGCCAAAC‑GGCAGCCCCTGATTACCCGGCACACAAA‑CA‑GCAATAAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCGGG‑CT‑TGCCG‑ATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                               <  2:1244134/228‑1 (MQ=255)
                                     acaGGGTGAGCCTAAGAATGAGGACCCCAAAGCTTAAGGAGAGGAACAGATTTCCACGATGAAAGCCGACGGCCA‑GCGCC‑ACGAAGAAGCCGATGGTGGAAAGATCGCAAAAC‑GGCAGCACTTGATAACCCGGCACACAAA‑CA‑CCAATTAAAAG‑GGT‑G‑AGTGGAATAAAAA‑GCAGG‑CT‑TGCAG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                            <  2:449254/228‑1 (MQ=255)
                                          gtaatgccaaaaaTGATGACGCCGGGAATTAGGGTGCGGACAAGATTCCAACGACGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCTATGGTGGCAAGATCGCCAAAC‑GGCAGCCCCTGATTCCCCGGCACCCAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGAAAAAAAA‑GCAGG‑CT‑CGCCG‑ATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                        <  2:1652872/216‑1 (MQ=255)
                                            gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCCGCAGCAACCACGGATCAAGGCCAATCAGGaa                                                                                                                                                                     >  1:1655689/1‑228 (MQ=255)
                                            gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                     >  1:1490414/1‑228 (MQ=255)
                                             aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                   >  1:213871/1‑229 (MQ=255)
                                             aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                   >  1:1366520/1‑229 (MQ=255)
                                              tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTc                                                                                                                                                                  >  1:338370/1‑229 (MQ=255)
                                              tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGG‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGCCCAATCAGGAACTc                                                                                                                                                                  >  2:485147/1‑229 (MQ=255)
                                                   aTAATGATGACACCCGAGCTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGCATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                                                                                              >  1:1828966/1‑228 (MQ=255)
                                                   aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGa                                                                                                                                                              >  1:678746/1‑228 (MQ=255)
                                                     aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGAcc                                                                                                                                                            >  1:1122626/1‑228 (MQ=255)
                                                     aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑AAGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGG‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                            >  1:442188/1‑228 (MQ=255)
                                                           gaCACCGGAAATTAGGGGGGGGAAAAAATTGCCCCAATGACCGGCACCCGCCA‑TCGCC‑ACAAAGAAGCGGATGGTGGCAAGATCGCAAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATAAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATcccac                                                                                                                                                                                                                    <  2:830525/168‑1 (MQ=255)
                                                           gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccgccgtat                                                                                                                                                                                                            >  1:100380/1‑169 (MQ=255)
                                                           gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccac                                                                                                                                                                                                                    >  1:831976/1‑168 (MQ=255)
                                                             cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACt                                                                                                                                                    >  1:926103/1‑228 (MQ=255)
                                                             cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                    >  1:278650/1‑228 (MQ=255)
                                                             cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                    >  1:1141309/1‑228 (MQ=255)
                                                              acCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑CAACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCGCACGGACAATCACGCAATCCTGACCGACAAACCTc                                                                                                                                                  >  2:773559/1‑229 (MQ=255)
                                                                              cGGAACAGTTTCCCGCAAGGCACGCCACCCCCCA‑TCGCC‑ACGAAGAAGCCGAGGGTGGCAAGATCGCCAAAC‑GGCAGCCCCTGATTCCCCGGCACAAAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGAAAAAAA‑CCAGG‑CT‑TGCCG‑ATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                     <  2:74345/178‑1 (MQ=255)
                                                                              cGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                     >  1:74467/1‑178 (MQ=255)
                                                                              cGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                  >  1:495123/1‑229 (MQ=255)
                                                                              cGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                  >  1:672173/1‑229 (MQ=255)
                                                                               ggAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACAACg                                                                                                                                 >  1:673498/1‑229 (MQ=255)
                                                                               ggAACAGATTTCCACGATCCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAG‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                       <  1:591002/191‑1 (MQ=255)
                                                                                              ggtGCACGGC‑ACCGCCAGTCGAC‑AAGAAGGAGCGAGTGGGGCAAAGATCCCAAAGC‑GCCAGCACATAATACCCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGG‑G‑GGGGGAAAAAAAG‑TCAGG‑CT‑TCCCG‑ATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                  <  2:441404/227‑1 (MQ=255)
                                                                                                tGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                               >  2:580485/1‑166 (MQ=255)
                                                                                                tGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                               <  1:581501/166‑1 (MQ=255)
                                                                                                 gCACGGCGCCCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCGGATTACCCGGCACACAAA‑CA‑GAAATTAAAAG‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                  <  1:310545/226‑1 (MQ=255)
                                                                                                  cACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                               >  1:1264169/1‑228 (MQ=255)
                                                                                                           ccccA‑TCGCCAACGAAGAGGCCGAGGGTGGCAAAACCACCAAGC‑GGCCACCCCTATTTCCCAGGCACA‑AAA‑CAGGCAATAAAATG‑GGT‑G‑GGGGGGAAAAAAA‑CCAGG‑CT‑TGCCG‑AT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                                                                                       <  2:28774/226‑1 (MQ=255)
                                                                                                                    gccagaagtacgCCGATGGTGAAAAGA‑CCCAAAACGGGAAGCCCCTGATAACCCGCCAACAAAA‑CA‑GAAATAAAAAT‑GTG‑G‑AGGGGAAAAAAAC‑CCAGT‑CT‑GGCC‑‑AT‑CACCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                              <  2:1261978/218‑1 (MQ=255)
                                                                                                                          aagaagCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                            >  1:546328/1‑225 (MQ=255)
                                                                                                                              agCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                    >  1:1053464/1‑229 (MQ=255)
                                                                                                                                   agggtgaCAAGCTCCACAACC‑GGAGGCCCCTGTTTCCCCGGCACAAAAA‑CA‑GCAATTAAAAG‑GGT‑G‑AGGGGGATAAAAA‑TCAGG‑CT‑GGCCG‑ATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                               <  2:545363/227‑1 (MQ=255)
                                                                                                                                          cAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                        >  1:1812170/1‑165 (MQ=255)
                                                                                                                                           agataccccAATC‑GGCAGCACCTGATTCCCCGGC‑CACAAAGCA‑GAAATGAAAAT‑GGT‑A‑AGGTGAAAAAAAA‑TCGGG‑CT‑TGCCG‑ATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                               <  2:33586/214‑1 (MQ=255)
                                                                                                                                              aTCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                        <  2:983120/129‑1 (MQ=255)
                                                                                                                                              aTCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                        >  1:984855/1‑129 (MQ=255)
                                                                                                                                                 gCCAACG‑GGCGGCACCTGATGACCCGGCACACAAA‑CA‑GAAATTAAAAG‑GGT‑G‑AGGGGGATAAAAA‑TCGGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                                                                                                 <  2:1602631/181‑1 (MQ=255)
                                                                                                                                                 gCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                                                                                                 >  1:1605363/1‑181 (MQ=255)
                                                                                                                                                  ccAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                  >  1:1863521/1‑179 (MQ=255)
                                                                                                                                                   cAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCATGCATAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                             >  1:1415751/1‑215 (MQ=255)
                                                                                                                                                   aaagc‑GGCAGCCCCT‑AATACCCTGGACAAAAA‑AA‑AGAAATAAAATGGGG‑G‑TGGGGGAAAAAAA‑GCGGG‑TT‑GGCCG‑TTAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                             <  2:1413320/211‑1 (MQ=255)
                                                                                                                                                         gccagcaCTTAATTCCCGCGCACACAAA‑CA‑‑CAAATAAAAT‑GGTGG‑GGGGGGAAAAAAA‑CCAGG‑TT‑TGCCA‑AT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                               <  2:677557/222‑1 (MQ=255)
                                                                                                                                                                         ccGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                                                                                               <  1:1918732/80‑1 (MQ=255)
                                                                                                                                                                         ccGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                                                                                               >  2:1915449/1‑80 (MQ=255)
                                                                                                                                                                          cgccacacaAA‑CA‑CAAATAAAAGT‑GGT‑G‑AGTGAGAAAAAAC‑ACAGG‑‑T‑TTCCC‑ACTCC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAGGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGTTGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                         <  2:1139319/226‑1 (MQ=255)
                                                                                                                                                                           ggacacaaaa‑CA‑GGAAATAAAAT‑GGG‑G‑GGGGGAAAAAAAATTCAGG‑CG‑TCCCA‑AT‑CC‑CCCCCCCTTATCTCCAGGCAGAACCGCCGGATCAAGGCCAATCAGGACCCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                         <  2:670994/222‑1 (MQ=255)
                                                                                                                                                                           ggCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccg                                                                                                                                                                                                                  <  1:773553/61‑2 (MQ=255)
                                                                                                                                                                           ggCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccg                                                                                                                                                                                                                  >  2:772199/1‑60 (MQ=255)
                                                                                                                                                                           ggCACACAAA‑CA‑GCAATTAAAAT‑GGA‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGCTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGCTTTTACCGCCa                                        >  2:33144/1‑229 (MQ=255)
                                                                                                                                                                             cacacaAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                       >  1:1607338/1‑228 (MQ=255)
                                                                                                                                                                              gcacaAA‑AC‑GCAAATAAAAA‑GGG‑GAAGGGGGAAAAAAA‑CCA‑G‑CT‑TGCCGCGCAGCACCCC‑CCGTGTGTCCGAGCGGCACCGCGGGATAAAGCCAAATCAGGAACTCTTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                       <  2:173820/228‑1 (MQ=255)
                                                                                                                                                                                acaAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccg                                                                                                                                                                                                                  <  1:87946/56‑2 (MQ=255)
                                                                                                                                                                                acaAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATaccaccg                                                                                                                                                                                                                  >  2:87802/1‑55 (MQ=255)
                                                                                                                                                                                  aaa‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCaa                                 >  2:802471/1‑229 (MQ=255)
                                                                                                                                                                                         gCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACAATCCATGATGGGTTTAATCACCCGTGGCATTAGCCGCATCACTGCCGCCGTTTTTCCCGCCAACTGCAGAacac                            >  2:379339/1‑229 (MQ=255)
                                                                                                                                                                                               taaaa‑GGT‑G‑AGTGGGAAAAAAA‑ACCGG‑GTGCGCCC‑AT‑CCACCCC‑CCGTTTGCCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGGGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                       <  2:672318/228‑1 (MQ=255)
                                                                                                                                                                                                  aT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTccc                                                                                                                                                                                                     >  1:398379/1‑51 (MQ=255)
                                                                                                                                                                                                  aT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTccc                                                                                                                                                                                                     <  2:397660/51‑1 (MQ=255)
                                                                                                                                                                                                         g‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                          >  1:1175965/1‑89 (MQ=255)
                                                                                                                                                                                                         g‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                              >  1:1312844/1‑117 (MQ=255)
                                                                                                                                                                                                         g‑AGGGGGAAAAAAA‑TCAGG‑CT‑TGCCG‑AT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGAccgc                                                                                                                                                          <  2:1173921/89‑1 (MQ=255)
                                                                                                                                                                                                          gggggggAAAAAAA‑TCAGA‑CG‑GCCGG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCGGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                              <  2:1310577/113‑1 (MQ=255)
                                                                                                                                                                                                           agggggAAAAAAA‑TCGGG‑TT‑TGCCG‑ACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                               <  2:63998/160‑1 (MQ=255)
                                                                                                                                                                                                           aGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                          >  1:1158172/1‑56 (MQ=255)
                                                                                                                                                                                                           aGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                          <  2:1156152/56‑1 (MQ=255)
                                                                                                                                                                                                           aGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                               >  1:64104/1‑163 (MQ=255)
                                                                                                                                                                                                             tGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTcc                                                                                                                                                                                                      <  1:155279/42‑1 (MQ=255)
                                                                                                                                                                                                             tGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTcc                                                                                                                                                                                                      >  2:155019/1‑42 (MQ=255)
                                                                                                                                                                                                               ggATAAAAA‑TCCGGTCT‑TGCCG‑AT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGACCTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCCCTGCCGCCGTTTTTACCGCCAGCTGCAATACCCCTTTGACATCGTAACCcg          <  1:33198/229‑1 (MQ=255)
                                                                                                                                                                                                                       a‑TCAGG‑CT‑TGCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAAGATg  >  1:1633696/1‑228 (MQ=255)
                                                                                                                                                                                                                                   gCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                         >  1:1050412/1‑36 (MQ=255)
                                                                                                                                                                                                                                   gCCG‑ATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                         <  2:1048561/36‑1 (MQ=255)
                                                                                                                                                                                                                                             cccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                              <  1:316589/103‑1 (MQ=255)
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                 
AGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCC‑ACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAA‑CA‑GCAATTAAAAT‑GGT‑G‑AGTGGGATAAAAA‑TCAGG‑CT‑TGCCG‑ATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATG  >  NC_000913/2173135‑2173571

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