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.85.9% 225.1 / 33.9 64intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/9);  new base . (39/16);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 7.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
*NC_0009132,173,3640C.85.9% 228.6 / 30.3 64intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/9);  new base . (39/16);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 7.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173122‑2173612
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                        
aGCATTGGCCGCCAGTTGGTGGGGCAAACCAAAAGTTTGGGTGGCGATCCACAGGGTGATGCTAATAATGATGACACCCGAGATTAAGGTGCGAAACAGATTTCCACGATGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTc                                                                                                                                                                                                                                                <  2:212213/251‑1 (MQ=255)
  cATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGTTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                              <  1:266261/251‑1 (MQ=255)
  cATTAGCCCCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAAGGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCACCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATAC‑‑CTCCGCCGTATGTccc                                                                                                                                                                                                                                              <  2:435756/251‑1 (MQ=255)
   gcgttcttccgatctggTGTGCAAACCAATAGTTTGCGTCGCGATCCACGGGGGGATGCTCATAATGAGGACCCCCGAGATTAAGGTGCGAAACAGATTTCCACGATGCACGGCGACCCCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCca                                                                                                                                                                                                                                             <  2:258398/236‑1 (MQ=255)
            cAGTGGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACg                                                                                                                                                                                                                                    <  1:87161/251‑1 (MQ=255)
              gTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                                  <  1:114021/251‑1 (MQ=255)
                  ggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCca                                                                                                                                                                                                                                             >  1:258398/1‑236 (MQ=255)
                       gCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                                         >  2:52958/1‑251 (MQ=255)
                          aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                       >  1:62998/1‑250 (MQ=255)
                            ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcgcat                                                                                                                                                                                                                    >  1:461787/1‑247 (MQ=255)
                              aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCCATTAAAATGGAGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGACGGATCAAGGCCAATc                                                                                                                                                                                                                   >  2:493831/1‑250 (MQ=255)
                                 aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGAGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGGCCCATCAGTAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                                                >  1:187217/1‑250 (MQ=255)
                                        gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                        <  1:8348/251‑1 (MQ=255)
                                           gcaaaccACAGGGTGATTCTCATAATGAGGACCCCCGAGATTAAGGTGGGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                     <  2:465153/246‑1 (MQ=255)
                                             gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCCCCCCGAAGCAACCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                                                                   >  1:578016/1‑251 (MQ=255)
                                                              tCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                                   >  1:210650/1‑250 (MQ=255)
                                                              tCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                                   >  1:165625/1‑250 (MQ=255)
                                                               cATAATGATGACCCCGGAGATTAAGGTGCGGAAAAGATTCCCAGGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                                 <  1:343145/251‑1 (MQ=255)
                                                                      atgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                          >  2:188613/1‑251 (MQ=255)
                                                                             ccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGc                                                                                                                                                                   >  2:169639/1‑251 (MQ=255)
                                                                                  aTTAAGGTGCGGAACAGCTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                              >  2:333475/1‑251 (MQ=255)
                                                                                  aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                             <  2:33792/188‑1 (MQ=255)
                                                                                  aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                             >  1:33792/1‑188 (MQ=255)
                                                                                    tAAGGTGCGGAACAGATTTCCCCGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAAATCGCCAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                             <  1:55560/250‑1 (MQ=255)
                                                                                    tAAGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                             >  2:156694/1‑250 (MQ=255)
                                                                                        gTGCGGAAGAGATTTCCACGATGCACGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAACCTCGCCAAACGGCAGCACCTGATAACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGAATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                        >  2:510112/1‑251 (MQ=255)
                                                                                        gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                        >  2:499730/1‑251 (MQ=255)
                                                                                            ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                                                    >  1:106301/1‑251 (MQ=255)
                                                                                            gaaaCAGATTCCCAGGAGGCACGGCGACCCCCATCGCCACGAAGAAGCCGATGTTGGCAAGATCGCCAATCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGTATAAAAA‑TCAGGCTGGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                                                    <  1:493831/249‑1 (MQ=255)
                                                                                                           gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                                     >  1:130239/1‑251 (MQ=255)
                                                                                                                        cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                                        >  2:432281/1‑251 (MQ=255)
                                                                                                                         nccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                                       >  2:282379/2‑251 (MQ=255)
                                                                                                                         cccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCTCCTGATTACCCGGCACAAAAACAGCAATTAAAATGGTGAGTTGTAAAAAAA‑CCAGGCTTGCC‑GAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTTAACCATCCATGTTGGGTTTAATCACCCg                                                                                                                       <  1:510112/250‑1 (MQ=255)
                                                                                                                          cccccGCCACGAAGAGGCCGAGGGTGGCAAGACCGCCAACCGGCAGCACTTAATTACCCGGCCCACAAACAGAAATAAAATTGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCgn                                                                                                                      <  2:278893/247‑2 (MQ=255)
                                                                                                                               gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                                 >  1:300216/1‑251 (MQ=255)
                                                                                                                               gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                                 <  1:188613/251‑1 (MQ=255)
                                                                                                                                ccgagaaaaagCCGAGGGGGGCAGGATCGCAAAAGGGAAGCCCCTGATAACCGGGACCACAAACGGAAATAAAAAGGGGGAGTGGGATAAAAA‑CCAGGCTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                <  2:130239/243‑1 (MQ=255)
                                                                                                                                  ccaaagaacCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAA‑TCAGGCTTGCC‑GAT‑CC‑CTCCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTTCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                                              <  1:245149/248‑1 (MQ=255)
                                                                                                                                    gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                             >  2:610458/1‑250 (MQ=255)
                                                                                                                                          gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                                                       >  1:460362/1‑250 (MQ=255)
                                                                                                                                            ggaggggggaaaGATCGCAAAACGGAAGCACCTGATTACCCGGGACCAAAAAAGAAATTAAAATGGGGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                                    <  2:165625/241‑1 (MQ=255)
                                                                                                                                             gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg                                                                                                   >  2:344181/1‑251 (MQ=255)
                                                                                                                                              aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGAACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCTCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgc                                                                                                  >  2:321516/1‑251 (MQ=255)
                                                                                                                                                   ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGccgccg                                                                                             <  1:156694/251‑1 (MQ=255)
                                                                                                                                                      aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                                          >  2:435420/1‑251 (MQ=255)
                                                                                                                                                               aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                  >  1:392880/1‑250 (MQ=255)
                                                                                                                                                                       gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                                         >  2:62925/1‑251 (MQ=255)
                                                                                                                                                                                 ccccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                <  1:26564/249‑1 (MQ=255)
                                                                                                                                                                                 aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                >  2:595630/1‑250 (MQ=255)
                                                                                                                                                                                   ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                                              >  1:97823/1‑186 (MQ=255)
                                                                                                                                                                                        acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                         >  2:538712/1‑250 (MQ=255)
                                                                                                                                                                                           aaaaccgcattaaaaaTGGTGAG‑GGGAAAAAAACTCAGG‑TTGCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                     <  2:401894/239‑1 (MQ=255)
                                                                                                                                                                                           aaaaCAGCAATTAAAAGGGTGAGTGGGATAAAAA‑TCAGGCTGGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTCGACATCGTAAcc                                                     <  1:432281/250‑1 (MQ=255)
                                                                                                                                                                                            aaaCAGCAATAAAAAGGGGGAGTGGGATAAAAA‑TCAGGCTGGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                                    <  1:282379/251‑1 (MQ=255)
                                                                                                                                                                                               cAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                                 >  1:207867/1‑251 (MQ=255)
                                                                                                                                                                                                    aTTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                                            >  1:622397/1‑251 (MQ=255)
                                                                                                                                                                                                             gTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACTCGCAAGGATGCCGATAAt                                   >  2:185331/1‑251 (MQ=255)
                                                                                                                                                                                                                 gTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                        >  2:39302/1‑194 (MQ=255)
                                                                                                                                                                                                                 gTGGGAAAAAAA‑CCAGGCTTGCC‑GATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                        <  1:39302/194‑1 (MQ=255)
                                                                                                                                                                                                                    ggATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                     >  2:625379/1‑194 (MQ=255)
                                                                                                                                                                                                                    ggATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                     <  1:625379/194‑1 (MQ=255)
                                                                                                                                                                                                                    ggATAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc                            >  1:411988/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTAAAAA‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCAt                          >  1:153949/1‑251 (MQ=255)
                                                                                                                                                                                                                           aa‑TCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCGGATAATCAGCCCCATcaca                      >  1:482189/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ggccgtccccaccgccgTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc              <  2:62998/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt   <  1:52958/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:109203/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  accgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  >  1:564175/2‑251 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                        
AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAA‑TCAGGCTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173122‑2173612

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