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.90.7% 162.2 / 10.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/4);  new base . (25/14);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01
*NC_0009132,173,3640C.90.7% 164.2 / 8.2 43intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (25/14);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGC  >  NC_000913/2173124‑2173610
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                      
cATTAGCCGCCAGTTGGGTGAGCAAAACAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCAGAGCTTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACAGCCATCGCCACGAAGAAGCCGATGGTGGCCAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                            >  1:201102/1‑251 (MQ=255)
 aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                            >  1:304107/1‑250 (MQ=255)
             gggggggtgCAAACAAATAGGTGGCGTCGCGATCCACAGGGTGTTTCTCAAAATGATGACCCCCGAGATAAGGGTGGAGAACATATTCCCACAATGCGCGG‑CGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAAAAACAGCAATAAAAAT‑GG‑T‑GAGTGGGAGAAAAATCAGG‑CGTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCg                                                                                                                                                                                                                               <  2:483491/245‑1 (MQ=255)
                     aaaaaCCAATAGTTTGCGGCGCGATCCACAGGGGGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                                       <  2:157288/249‑1 (MQ=255)
                      cAAACCAATAGTTTGCGCCGGGATCCAAGGGGTGATGCCCATAATGATGACCCCCGAGATTAAGGGGCGGAACAGATGTCCAGGA‑GGACGGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                      <  2:391231/251‑1 (MQ=255)
                            aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAATCGGCAGCACCTGATTGCCCGTCACACAAGCAGCAAGTAAAAT‑GT‑T‑GAGGGGTGTAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                      <  1:490282/229‑1 (MQ=255)
                            aaTAGTTTGAGTCGCGATCCCCAGGGCGATGCTAATAATGATCACACCCGAGCTTAAGGTCCGGCACAGCTTTCCACGATGCACGG‑CGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAACT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                      >  2:490282/1‑229 (MQ=255)
                                      gTCGCGCTCCACAGGGTGATGATCATCATGATGACACCCGAGCTTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGCGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGCGGGATAAAAATCAGG‑CTTGC‑CGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGTACTCCTg                                                                                                                                                                                                      >  2:232379/1‑251 (MQ=255)
                                        cgcgATCCACAGGGTGATGCTCAAAAGGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                                    <  1:276317/251‑1 (MQ=255)
                                                  aGGGTGATGCTCATAATGATGCCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATAGCCACGAAGAAGCCGATGGTGGCAAGATCGCCCAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                                          >  2:471784/1‑251 (MQ=255)
                                                                aaTGAGGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCCCGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                                                                                                                            <  1:335992/251‑1 (MQ=255)
                                                                      gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                        >  2:64625/1‑249 (MQ=255)
                                                                                   aaGGGGCGAAACAGATTCCCAGGGTGCCCGG‑CGCCCCCCACCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGGACCGGATTACCCGTCACACAAACAGCAATTAAAAG‑GG‑T‑GAGGGGGATAAAAATCAGG‑CTGGC‑CGATTCCACC‑‑GCCGTATGTCCCAGGAGCAACGCGGGATCAGGGCCAATCAGGAACTCCTTACCGCCGATCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                                         <  1:42569/251‑1 (MQ=255)
                                                                                         cGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                                   >  1:247149/1‑251 (MQ=255)
                                                                                                            ccACGG‑CGCCCGCCATCGCCACGAAGAGGCCGAGGGTGGCAAGACCGCCAAACGGCGGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAAT‑GT‑T‑GAGTGGGAAAAAAATCAGG‑CTTGC‑CGAGCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACGACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                                <  1:471784/250‑1 (MQ=255)
                                                                                                               ccggccaccgccctcccaccGAAGAACCCGATGGGAGCAAGATCGCCAACCGGAAGCCCCGGATTACCCGCCCCACAAACAGCA‑TTAAAAT‑GGGT‑GGGGGGAAAAAAAACCAGG‑CTTGC‑CAATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                              <  2:367229/232‑1 (MQ=255)
                                                                                                                       cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGCCACCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGAAAAAACTCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                                      >  1:34146/1‑251 (MQ=255)
                                                                                                                           aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                         >  2:60789/1‑180 (MQ=255)
                                                                                                                           aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                         <  1:60789/180‑1 (MQ=255)
                                                                                                                           aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg                                                                                                                  >  2:110741/1‑251 (MQ=255)
                                                                                                                             cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACCAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                                    >  1:188533/1‑183 (MQ=255)
                                                                                                                             cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                                 >  2:294041/1‑250 (MQ=255)
                                                                                                                              gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCAACTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                               >  1:348214/1‑251 (MQ=255)
                                                                                                                                          ccGATGGGGGAAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGAAATAAAAAG‑GG‑G‑GAGTGGGATAAAAACCAGG‑TTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                   <  1:218844/251‑1 (MQ=255)
                                                                                                                                             aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgc                                                                                                >  2:4975/1‑251 (MQ=255)
                                                                                                                                                 gggaaaGATCCCCAAACGGCAGCACCTGATTACCCGGCTCACAAACAGCATTTAAAAG‑GG‑T‑GAGTGGGAAAAAAACCGGG‑CC‑GCGCGATACCACC‑‑GCCGTATGACCCAGAAGCAACGCCGGACCAAGGCAAAACAGTAACTCCTGACCGCCGAACTTCGCCGGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCTTCCATTATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                            <  1:496304/247‑1 (MQ=255)
                                                                                                                                                             cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCCTCACTGCCGCCGTTTTTACCGcc                                                                                >  2:422600/1‑251 (MQ=255)
                                                                                                                                                                              ttACCCGGCACACAATCAGCAATTACAAT‑GG‑T‑GAGCGGGCTAACACTCCGG‑CCTGC‑AGATACCACC‑‑GCCGTATGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCATGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                               >  2:515082/1‑251 (MQ=255)
                                                                                                                                                                               tACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                              >  1:163183/1‑251 (MQ=255)
                                                                                                                                                                                 cccGGCACACAAACAG‑AATTAAAATTGG‑T‑GGGTGGAATAAAAA‑CAGGGCTGGC‑CG‑‑ACCCCACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGAGGGGCGTTAAACCATCCATGATGGGTTTAATCTCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                            <  1:232308/251‑1 (MQ=255)
                                                                                                                                                                                     gCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                                >  2:72872/1‑179 (MQ=255)
                                                                                                                                                                                     ccacacaAACAGCAATTAAAAT‑GG‑T‑GATTGGGAAAAAAATCAGT‑CTTGC‑CG‑‑ACCCCACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                                <  1:72872/178‑1 (MQ=255)
                                                                                                                                                                                       ccaaaaaaaggcaatAAAAT‑GG‑TTGGGTGGAAAAAAAACCAGG‑CTTAC‑CCATCCCACC‑‑GCCGTATGACCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACGTCGCCTGGAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGAGGGGTTTAATCACCCGTTGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                       <  1:494709/237‑1 (MQ=255)
                                                                                                                                                                                        aaaaaaCAGAAATTAAATG‑GG‑G‑GAGGGGGAAAAAAATCAGG‑CTGGC‑CGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTaa                                                     <  1:513532/248‑1 (MQ=255)
                                                                                                                                                                                         acaAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                    >  2:228592/1‑251 (MQ=255)
                                                                                                                                                                                           aaaCAGAAATTAAAAT‑GG‑T‑GAGTGGGAAAAAAACAAGG‑CTGCC‑CGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                                  <  1:3638/251‑1 (MQ=255)
                                                                                                                                                                                             aCAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CATGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGCTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTACCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGCCATCGTAACCcgc                                                >  1:314551/1‑251 (MQ=255)
                                                                                                                                                                                               agaaataaaaaG‑GG‑T‑GAGGGGGATAAAAACCAGG‑CTTGC‑CGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGa                                              <  1:64625/244‑1 (MQ=255)
                                                                                                                                                                                                        aaT‑GG‑T‑GAGTGGGATAAAAATAAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTCACCAGCGCGGAGGCCg                                      >  1:174332/1‑250 (MQ=255)
                                                                                                                                                                                                            gg‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATa                                   >  1:455691/1‑250 (MQ=255)
                                                                                                                                                                                                                  aGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                >  1:17773/1‑249 (MQ=255)
                                                                                                                                                                                                                              aaaaacaggcctcgccatcccccccccCGAATGTCCCAGCAGCAACGCGGGACAAAGGCCAACAAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATca                      <  2:247149/230‑1 (MQ=255)
                                                                                                                                                                                                                                  gg‑CTTGC‑CGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGGTGCCGATACTCAGCCCCATCACAAAAccc              >  2:284106/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                      >  1:475263/1‑156 (MQ=255)
                                                                                                                                                                                                                                               ataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTaa                                                     >  1:53574/7‑203 (MQ=255)
                                                                                                                                                                                                                                                   caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCATGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATAAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCCTCACAAAACCCACGGTGACAGGc  >  1:354125/3‑250 (MQ=255)
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                      
CATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGG‑CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GG‑T‑GAGTGGGATAAAAATCAGG‑CTTGC‑CGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGC  >  NC_000913/2173124‑2173610

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