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.97.7% 181.8 / ‑4.2 43intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (28/14);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.100.0% 187.8 / NA 43intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (28/15);  total (28/15)

GGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACA  >  NC_000913/2173139‑2173607
                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                    
gggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                                                                                                                           <  1:462798/251‑1 (MQ=255)
   tgtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CCGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                         >  1:451934/1‑250 (MQ=255)
   ggtgaaaaCCAATAGTTGGGGTCGCGATCCACAGGGTGAGGCTAATAATGATGACCCCGGAGATTAAGGGGCGGAACAGATTTCCAGGAGGCGCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                        <  2:350644/246‑1 (MQ=255)
              aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                               >  1:153486/1‑249 (MQ=255)
              aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GACGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                               >  1:241459/1‑249 (MQ=255)
                   ttgcggcgcaacccACAGGGGGATGCTCATAATGATGACCCCCGAGATAAAGGTGCGGAACAGATTCCCACGATGCCCGGCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAAT‑CAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGG‑AAc                                                                                                                                                                                                        <  2:229748/238‑1 (MQ=255)
                           cgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGAccg                                                                                                                                                                                               >  2:396890/1‑251 (MQ=255)
                             aTCCACAGGGTGATGCTCATAATGAGGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGAccgcc                                                                                                                                                                                             <  1:396890/251‑1 (MQ=255)
                               ccACAGGGGGATGCCCATAATGAGGACCCCCGAGATAAAGGGGCGGAACAGATTTCCACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCCCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGGGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGa                                                                                                                                                                                           <  2:291036/251‑1 (MQ=255)
                                              cATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCAGGGCGCCCGCCATCGCCACAAAGAAGCCGATGGTGGCAAAAGCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAC‑CAGCCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                            <  1:218016/251‑1 (MQ=255)
                                                    gatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                      >  1:317545/1‑251 (MQ=255)
                                                     atgaCACCCGAGATTCAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                     >  1:418732/1‑251 (MQ=255)
                                                        cccccggatataaaGGGGCGGACCAGATTTCCGCGTTGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCCCCGGATTACCCGGCACAAAAACGGCAATTAAAATGGGGAGTGGGATAAAACT‑CAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                                  <  2:216723/239‑1 (MQ=255)
                                                        acacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                                  >  1:225249/1‑251 (MQ=255)
                                                                    aaGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCAGCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCATTTAAAATGGTGTGGGGGATAAAATT‑CAGGTTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAt                                                                                                                                                      <  1:552554/251‑1 (MQ=255)
                                                                           ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTCGCGATGGGCGtt                                                                                                                                               >  2:389607/1‑251 (MQ=255)
                                                                                    ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                      >  2:184473/1‑251 (MQ=255)
                                                                                           gtcacgtggcccccaccgcGCCAGGAAGAAGCCAATGTGGGAAAGACCGCTAAACGGCCGCCCCAGATTCCCCGCCACACAAACAGCAATTAAAGTGGTGATGGGGAAAAAACT‑CGGGCTTGCCAATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                               <  2:177325/233‑1 (MQ=255)
                                                                                           aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGAATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAACCGCCTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                               >  1:49109/1‑251 (MQ=255)
                                                                                                cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAG‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa                                                                                                                          <  1:62532/251‑1 (MQ=255)
                                                                                                    gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                       >  1:45497/1‑250 (MQ=255)
                                                                                                     aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc                                                                                                                     >  2:95596/1‑251 (MQ=255)
                                                                                                              gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACATGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                            >  2:261592/1‑251 (MQ=255)
                                                                                                                cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                           >  2:186550/1‑250 (MQ=255)
                                                                                                                       aagCCGATGGTGGCAAAATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                   <  1:223002/251‑1 (MQ=255)
                                                                                                                                    cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAAAAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                      >  2:101382/1‑251 (MQ=255)
                                                                                                                                                  ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAACAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACAATCCATGATGGGTTTAATCACACGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                        >  1:289880/1‑251 (MQ=255)
                                                                                                                                                     agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                      >  2:527627/1‑250 (MQ=255)
                                                                                                                                                                     gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                     >  1:528951/1‑251 (MQ=255)
                                                                                                                                                                     gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                     >  2:187697/1‑251 (MQ=255)
                                                                                                                                                                       acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                    >  2:112415/1‑250 (MQ=255)
                                                                                                                                                                         acaAACAGCAGTTGCAATGGAGCGTGGGCTAAAACT‑CCGGCTTGCCGATACCACC‑‑GCAGTCTGTCCCAGCAGCAACGCAGGATAAAGGACAATCAGA‑CACTACTGACCGCCGAACTTCGCCTGTAAACGACTACAAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTCGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCACTAATCCTTTGACAtagaaac                                                 >  2:263092/1‑245 (MQ=255)
                                                                                                                                                                            aaCAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcccc                                              >  2:208546/1‑250 (MQ=255)
                                                                                                                                                                             aCAGCAATTCAAATGGTGAGTGGGATAAACAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                             >  2:245518/1‑251 (MQ=255)
                                                                                                                                                                             aCAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGATCAAGTCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAACAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                             >  1:247153/1‑251 (MQ=255)
                                                                                                                                                                              acaaaattgaaagtggtggggggAAAAAAATACAGTCTTCCCGATACCACC‑‑GCCGTATGCCCCAGAAGCACCGCCGGATCAAGGCCAATCAGG‑AACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCCTCACTGCCGCCGTTTTTACCGCCAGCTTCAATACACCTTTGACATCGTAACCcgc                                             <  2:451934/231‑1 (MQ=255)
                                                                                                                                                                                        aaTGGAGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTACCAGCAGCAAAGCCGGATCAAGGCCACTCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGAGCTTGCTTCGCGATGGGCGGTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCACTACACCTTTGACATCGTAACCCGCGAGGATGCCg                                   >  2:534139/1‑250 (MQ=255)
                                                                                                                                                                                            gTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAAt                              <  1:150476/251‑1 (MQ=255)
                                                                                                                                                                                               aGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATagtc                             >  2:108504/1‑246 (MQ=255)
                                                                                                                                                                                                        aaaaT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                  >  1:460111/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaaT‑CAGACTTGCCGATACCACC‑‑GCCGTATTTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATTGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                  <  1:2505/251‑1 (MQ=255)
                                                                                                                                                                                                        aaaaC‑CAGCCTTGCCAATACCACCC‑‑CCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                  <  1:250727/251‑1 (MQ=255)
                                                                                                                                                                                                         aaaaaaccaggcttcccgttcccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGT‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                    <  2:427441/227‑1 (MQ=255)
                                                                                                                                                                                                         aaaT‑CACGCTTGCCGATACCACC‑‑GCCGCATGTCCCAGCAGCAAAGCCGGATCACGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAAAGACTACGAGCTTGCTTAGCGATGGGCGTTATACCATCAATGCTGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGACGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCTGATAATCCGACCCATcac                  >  2:357241/1‑250 (MQ=255)
                                                                                                                                                                                                                 gggcttgcctacacctccgccgTGTGTCCCAGCAGCAACGCCGGATCAAGGCCAA‑CAGGTAACTCCTTACCGCCGAATTTCGCCTGTAAACGACTCCGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTATTCGCCCTTGTCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACATCGTAACCCGCGAGGATGCCGATAATCATCCCCATCACAAAAcc            <  2:241257/235‑1 (MQ=255)
                                                                                                                                                                                                                           atcccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACa  <  1:10849/245‑1 (MQ=255)
                                                                                                                                                                                                                           ataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc   >  2:160804/7‑250 (MQ=255)
                                                                                                                                                                                                                            taccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc   >  1:440064/6‑249 (MQ=255)
                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                    
GGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACA  >  NC_000913/2173139‑2173607

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: ‑