Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.95.2% 273.3 / 4.9 62intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (39/20);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 4.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.95.2% 272.8 / 4.1 62intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (39/20);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 4.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173127‑2173613
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                         
tAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                            >  1:117559/1‑251 (MQ=255)
      ccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                                      >  1:243582/1‑251 (MQ=255)
            gggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                                 <  1:144966/250‑1 (MQ=255)
              gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTCCCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                              <  1:183275/251‑1 (MQ=255)
                  gCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                                                                                                                                           <  1:521247/250‑1 (MQ=255)
                     aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                        >  2:522554/1‑250 (MQ=255)
                          aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGACGGATCAAGGCCAATc                                                                                                                                                                                                                    >  2:476802/1‑249 (MQ=255)
                                 gCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTc                                                                                                                                                                                                                                                 <  1:397412/213‑1 (MQ=255)
                                 gCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTc                                                                                                                                                                                                                                                 >  2:397412/1‑213 (MQ=255)
                                         aTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                   <  1:199415/251‑1 (MQ=255)
                                                 ggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAACCAGCAATT‑AAAATGGGGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                           >  2:294660/1‑251 (MQ=255)
                                                         tCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                                    >  1:73906/1‑250 (MQ=255)
                                                                 atgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                      >  1:486532/1‑224 (MQ=255)
                                                                 aggaCACCCGAAATTAAGGTTCGAAACAGATTTCCGCGATGCACGGCTACC‑GCCATCGCCACGAAGAAG‑CCGAGGGTGGCAAAATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑TCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                      <  2:486532/222‑1 (MQ=255)
                                                                     cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcc                                                                                                                                                                       >  2:330773/1‑251 (MQ=255)
                                                                       cccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                                                                                                                                                     >  2:114769/1‑251 (MQ=255)
                                                                         gaaattagagtggcggacaaGATTCCCACGAGGCACGG‑GACCCGCCATCGCCGGGAAGAAG‑CCGATGGGGGCAAGACCGCCAAAC‑GGAAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAAGGGTGAGGGGGAAAAAAATCAGG‑CTGGCC‑GAGACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                                                                                                                                   <  2:222062/232‑1 (MQ=255)
                                                                         cGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGGGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGGAAACGACTACGAACCTGCt                                                                                                                                                                   >  2:213473/1‑251 (MQ=255)
                                                                             aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                               >  1:58287/1‑251 (MQ=255)
                                                                               tAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                              >  1:426919/1‑250 (MQ=255)
                                                                                   gTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                         >  1:221964/1‑251 (MQ=255)
                                                                                       ggAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                                                     >  1:509043/1‑251 (MQ=255)
                                                                                          aCAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                                                                                                                                   >  2:379050/1‑250 (MQ=255)
                                                                                            aGATTCCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                                <  1:227362/251‑1 (MQ=255)
                                                                                              attcccccGATCCACGGCGCCC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGAAAAAAATCAGG‑CTTGCC‑GAAACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                              <  1:476802/244‑1 (MQ=255)
                                                                                               tttCCACGATGCCCGGCGCCC‑GCCATCGCCCCGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                              <  1:183487/250‑1 (MQ=255)
                                                                                                 cccaagatccccGGCAACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTCCCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                            <  1:84573/239‑1 (MQ=255)
                                                                                                  ccACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                          >  2:476944/1‑251 (MQ=255)
                                                                                                             ggCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                                   <  2:183680/151‑1 (MQ=255)
                                                                                                             ggCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                                   >  1:183680/1‑151 (MQ=255)
                                                                                                             ggCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTa                                                                                                                                >  2:539908/1‑250 (MQ=255)
                                                                                                                cgccgcccaccGCCACAAAAAAG‑CCGTGGGTGGAAAGATCCCAAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAAGGGTGAGGGGGAAAAAAATCGGG‑TTGGCC‑GATACCACC‑‑GCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                             <  2:11443/241‑1 (MQ=255)
                                                                                                                     gCCATCGCCACGAAGAAG‑CAGATGGGGGAAAGATCGCCAAAC‑GGCACCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAATACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                                        >  1:476387/1‑251 (MQ=255)
                                                                                                                      ccATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTACTCACCCGt                                                                                                                       >  1:399549/1‑251 (MQ=255)
                                                                                                                          cGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                                    >  1:355574/1‑250 (MQ=255)
                                                                                                                          cGCCACGA‑GAAGGCCGATGGTGGCAAGA‑CCCCAAACCGGAAGCACCTGATTCCCCGGCACAAAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGC‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCa                                                                                                                   <  2:307465/251‑1 (MQ=255)
                                                                                                                               cGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                              >  1:139395/1‑251 (MQ=255)
                                                                                                                                gaagaaG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGAGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                              >  1:26382/1‑250 (MQ=255)
                                                                                                                                 aagaag‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcag                                                                                                            <  1:379050/251‑1 (MQ=255)
                                                                                                                                 aagaag‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                              >  1:452621/1‑249 (MQ=255)
                                                                                                                                     ag‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                                                        <  1:330773/251‑1 (MQ=255)
                                                                                                                                          gATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg                                                                                                    >  2:225453/1‑251 (MQ=255)
                                                                                                                                                   aagagcgcaaacc‑GGCAGCACCGGATTACCCGGCACAAAAACAGCATTT‑AAAAGGGTGAGTGGGAAAAAAACCAGG‑CTTGCG‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                                           <  1:225453/239‑1 (MQ=255)
                                                                                                                                                   aaGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                                           >  1:187463/1‑251 (MQ=255)
                                                                                                                                                         gCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTAATCACCCGAGGCATTAGCAGCATCACTGCCGCCGTTTTTCCCg                                                                                     >  2:209531/1‑251 (MQ=255)
                                                                                                                                                         gCCAAAC‑GCCACCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGAATAAAAATCAGG‑CTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCg                                                                                     <  1:539908/251‑1 (MQ=255)
                                                                                                                                                           cAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCTGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  1:85683/1‑251 (MQ=255)
                                                                                                                                                               c‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCt                                                                               >  1:99834/1‑251 (MQ=255)
                                                                                                                                                                 ggCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                              >  1:260156/1‑251 (MQ=255)
                                                                                                                                                                       aCCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                        >  1:374749/1‑251 (MQ=255)
                                                                                                                                                                            aTTACCCGGCCCACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCtt                                                                   <  1:346190/251‑1 (MQ=255)
                                                                                                                                                                                  cGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                             >  1:256807/1‑251 (MQ=255)
                                                                                                                                                                                   ggCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                            >  2:558951/1‑251 (MQ=255)
                                                                                                                                                                                   cgaccaaaaaCAGCAATTTAAAAGGGTGGGGGGGATAAAAA‑CAGGGCTGGCC‑G‑TAAC‑CCCCGCCTTATGTCCCAGCAGCACCGCCGGATAAAGGCCAATCGGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATGAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACaa                                                             <  2:169579/244‑2 (MQ=255)
                                                                                                                                                                                        acaAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                         >  1:290794/1‑249 (MQ=255)
                                                                                                                                                                                          aacaaGCAATA‑AAAAGGGGGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                                     <  1:476944/247‑1 (MQ=255)
                                                                                                                                                                                                 aaTT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                                               >  1:465834/1‑250 (MQ=255)
                                                                                                                                                                                                           ggTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATa                                      >  2:340023/1‑250 (MQ=255)
                                                                                                                                                                                                             ggagggggAAAAAAATCAGG‑‑TTGCCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                   <  2:73906/246‑1 (MQ=255)
                                                                                                                                                                                                                gTGGGATAAAAATCAGG‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGc                                >  2:435153/1‑251 (MQ=255)
                                                                                                                                                                                                                                g‑CTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATCCCGATAATCAGCCCCATCACAAAAccc                 >  1:676/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATTACAAAACCCACg              <  1:387667/250‑1 (MQ=255)
                                                                                                                                                                                                                                             taccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCAATCACAAAACCCACGGTGAc         >  2:449566/6‑249 (MQ=255)
                                                                                                                                                                                                                                                  ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  <  1:522554/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                         
TAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACC‑GCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATT‑AAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173127‑2173613

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