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.4% 137.8 / 12.4 41intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base . (22/13);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.85.4% 134.0 / 16.2 41intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base . (22/13);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGC  >  NC_000913/2173128‑2173610
                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                                      
aGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                         <  1:580593/250‑1 (MQ=255)
  ccgccAGTTGGGGGTCCAAACCAATATTTTGCGTCGCGACCCCAAGGGTGATGCTCAGAATGATGACCCCCGAGTTTAAGGTGCGGAACA‑AATTCCCACGTGGCCGGGCGACCGCCATCCCCACGAAGAGGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGTATAAAAATCAGTCTTGCCTATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                      <  1:254637/251‑1 (MQ=255)
           gggtggcaaaaCCAAAAGTTTGCGTCGCGATCCACAGGGTGATGCCCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGTATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCATCAACGCCGGa                                                                                                                                                                                                                             <  1:361488/243‑1 (MQ=255)
            ggggtgCAAACCAATAGTTGGCGTCGCGATCCACAGGGTGAGGCCCATAAGGATGACACCCGAGATTAAGGTGCGGAACA‑GATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                                            <  1:507804/248‑1 (MQ=255)
                tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGCTCa                                                                                                                                                                                                                          >  2:580290/1‑249 (MQ=255)
                 gCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCAAAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                                       <  1:48918/251‑1 (MQ=255)
                  gaaaCCAATAGTTTGCGTCGCGATACACAGGGTGATGCTCATAATGATTACACCCGAGATTAAGGAGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCCATTAAAATGGTGAGTGGGAAAAAAATCAGGCTTGCCGACAC‑‑CCCCGCCGTATGTCCAAGAAGCAACGCCGGATCACGGc                                                                                                                                                                                                                      >  2:664839/2‑251 (MQ=255)
                     aCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCACGGcaa                                                                                                                                                                                                                    >  2:541823/1‑248 (MQ=255)
                            gTTTGCGTCGCGACCCACAGGGTGATGCTCATAATGAGTACACCCGAGATAAAGGGGCGGAACA‑GATTTCCACGATGCGCGGCGACCGCCAGCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAACCGGCAGCACCTGATTCCCCGGCACACAAACACCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                                                            <  2:271664/251‑1 (MQ=255)
                                    cgcgATACACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCCACAACGCCGGATCAAGGCCAATCAGGAACTCCttac                                                                                                                                                                                                    >  2:349300/1‑248 (MQ=255)
                                       gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GCTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGACCCCCTGAccgc                                                                                                                                                                                                 >  2:222193/1‑251 (MQ=255)
                                           cacaGGGTGATGCTAATAATGAGGACCCCCGAGATTAAGGTGCGGAACA‑GATTCCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                                                                                                                             <  1:44566/251‑1 (MQ=255)
                                             caGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCAACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGacct                                                                                                                                                                                           >  2:455169/1‑248 (MQ=255)
                                              aGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                                          >  1:182700/1‑251 (MQ=255)
                                                ggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTcc                                                                                                                                                                                        >  2:64104/1‑250 (MQ=255)
                                                         cATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                               >  1:734910/1‑251 (MQ=255)
                                                            aaTGATGACCCCCGAAATAAAGGGGCG‑AACAAGATTTCCACGATGCAGGGCACCCGCCATCGCCACAAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                                                                                                                            <  1:675719/251‑1 (MQ=255)
                                                                            aTTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGc                                                                                                                                                            >  2:520262/1‑251 (MQ=255)
                                                                              tAAGGTGGGGAAAA‑GTTTTACACGAGGCAGGGGGACCGCCATCGCCACGAAGAAGCTGAGGGTGGAAAGATCGCAAAAGGGCAGCACCTGATTACCCGGCACAAAAACAGCAATTAAAATGTGGAGGGGGATAAAAATCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCg                                                                                                                                                           <  1:541823/250‑1 (MQ=255)
                                                                                 ggtggggacaa‑GATTTCCACGATCCAGGGCGCCCGCCTCCCCCACAAAGAAGCCGAGGGTGGCAAGACCGCCAAACGGCAGCACCGGATTCCCCGGCACACAAACAGAAATAAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCATTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAGCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATgg                                                                                                                                                       <  1:222193/241‑1 (MQ=255)
                                                                                                      aagcccgggCGCCCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                                   <  1:655142/244‑1 (MQ=255)
                                                                                                                gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                          >  1:514804/1‑250 (MQ=255)
                                                                                                                   cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                                      >  1:341723/1‑251 (MQ=255)
                                                                                                                               gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                           >  1:261928/1‑250 (MQ=255)
                                                                                                                                                          aaccggcacaccttGATTACCCGGCACAAAACCAGCAATTAAAATGGTGGGTGGAAAAAAAATCAGGCTTCCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCa                                                                               <  1:349300/238‑1 (MQ=255)
                                                                                                                                                          aaaCGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                            <  1:295168/238‑1 (MQ=255)
                                                                                                                                                          aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                            >  2:295168/1‑238 (MQ=255)
                                                                                                                                                           aaggccagcaCCTGTTTCCCCGGCACAAAAAAAGAAATAAAAATGGTGAGTGGGAAAAAAATCAGGTTTGCGGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                                                              <  1:603909/246‑1 (MQ=255)
                                                                                                                                                            aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                                             >  2:745805/1‑251 (MQ=255)
                                                                                                                                                                           tACCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACAAGCTTGCTTAGCGATGAGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCGTTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCtttt                                                              >  2:566966/1‑250 (MQ=255)
                                                                                                                                                                           accccGGAACAAAAACAGAAATAAAAATGGGGAGGGGAAAAAAAACCAGGCTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                               <  1:64104/248‑1 (MQ=255)
                                                                                                                                                                             cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                            >  2:170911/1‑251 (MQ=255)
                                                                                                                                                                                 gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGCCATCg                                                        >  1:177646/1‑251 (MQ=255)
                                                                                                                                                                                      caAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                   >  1:303744/1‑251 (MQ=255)
                                                                                                                                                                                        aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCt                                                                            >  1:349016/1‑224 (MQ=255)
                                                                                                                                                                                           aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                               >  1:242050/1‑250 (MQ=255)
                                                                                                                                                                                                  aaaaAGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCAAGCTGCAATACACCTTTGACATCGTAACCCGCGGGGTTGcc                                       >  1:25778/1‑251 (MQ=255)
                                                                                                                                                                                                             tGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGCTAATCAGcc                            >  1:500423/1‑251 (MQ=255)
                                                                                                                                                                                                                         cAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAc                <  1:408810/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAACAcc               >  2:730592/1‑251 (MQ=255)
                                                                                                                                                                                                                          aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAcc               >  1:233434/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc      >  1:411986/8‑251 (MQ=255)
                                                                                                                                                                                                                                       cccccccGCCGTATGTCCCAGCAGCAACCCCGGATCAAGGCCAGTCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAg    <  1:455169/248‑1 (MQ=255)
                                                                                                                                                                                                                                         ccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGc  >  1:784097/4‑251 (MQ=255)
                                                                                                                                                                                                                                         ccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGc  >  1:561821/4‑251 (MQ=255)
                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                                      
AGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGC  >  NC_000913/2173128‑2173610

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