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.83.8% 111.4 / 19.5 37intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base . (27/4);  total (29/8)
Fisher's exact test for biased strand distribution p-value = 1.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01
*NC_0009132,173,3640C.83.8% 113.1 / 17.7 37intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base . (27/4);  total (29/8)
Fisher's exact test for biased strand distribution p-value = 1.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173130‑2173612
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                        
ccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCCTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAACCGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCagcagc                                                                                                                                                                                                                                        >  1:545955/1‑251 (MQ=255)
       ttGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCg                                                                                                                                                                                                                                 >  1:116436/1‑251 (MQ=255)
           gtgggaaaaCAAATAGTTTGCGCCGCGACCCAAAGGGTGATGCTCATAAGGAGGACACCCGAGATTAAGGGGCGGAACAGATTTCCACAATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                             <  2:47932/245‑1 (MQ=255)
                  aaCCAATATTTTGCGTCGCGATCCACAGGGTGATGCTCATAAGGATGACACCCGAGATTAAGGTGCGGAACAAATTCCCACGATGCCCGGCGACCGCCCTCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTCCCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                                      <  1:64234/251‑1 (MQ=255)
                   aCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                                      >  1:337668/1‑250 (MQ=255)
                          gTTTGCGTCGCGATCCACAGGGTGAGTCTCATAATAATGACACCCGAGATTAAGGGGCGGAAAAGATTTCCCCGATGCACGGCGACCGCTATCGCCACGAAGAAGCCGATGGGGTCAAGA‑TCGCCAAACGC‑CAGCACCTGATTACCCGGCATACAAATA‑GCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTGGCC‑GAGACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                                                              <  1:238286/251‑1 (MQ=255)
                                gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTACGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                        >  2:249989/1‑251 (MQ=255)
                                      aTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                  >  1:437836/1‑251 (MQ=255)
                                                ggATGCCCATAATGATGACACCCGAGATTAAGGTGGGGAACAGATTCCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCCCCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                                                                                                                                        <  1:33718/250‑1 (MQ=255)
                                                             gatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                           >  1:513286/1‑251 (MQ=255)
                                                                gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACg                                                                                                                                                                          >  2:171390/1‑249 (MQ=255)
                                                                               ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                                         >  1:118902/1‑251 (MQ=255)
                                                                                gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACCAGCCTGCTTAGCGATggg                                                                                                                                                        >  2:507603/1‑251 (MQ=255)
                                                                                   cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                                     >  1:218471/1‑251 (MQ=255)
                                                                                               ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                         >  2:558164/1‑251 (MQ=255)
                                                                                                   gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                     >  1:27729/1‑187 (MQ=255)
                                                                                                             gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGATTTAATc                                                                                                                            >  1:288837/1‑250 (MQ=255)
                                                                                                                cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACATGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAATCCTGACCGCCGAACTTCGCCTGTAAAAGACTACGAGCCTGCTTAGCGATGGGCGTTAAACAATCCATGATGGGTTTAATCAccc                                                                                                                        >  1:527560/1‑251 (MQ=255)
                                                                                                                           cGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                             >  1:217186/1‑251 (MQ=255)
                                                                                                                           cGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                             >  1:328341/1‑251 (MQ=255)
                                                                                                                                 agCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                                                       >  1:525390/1‑251 (MQ=255)
                                                                                                                                     ggagggggCAAGAATCGC‑AAACGGGCAGCACTTGATT‑CCCGGGCCACAAACAAGAAATAAAAAGGGTGAGGGGAAAAAAAACCAGG‑TTGCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                                    <  2:116436/245‑1 (MQ=255)
                                                                                                                                      aTGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGAAGCATCACTgc                                                                                                  >  2:550801/1‑251 (MQ=255)
                                                                                                                                                      ccAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAACCGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTAATCACCCGTGGCCTTAGCAGCATCACTGCCGCCGTTTTTTCCGc                                                                                   >  2:132346/1‑251 (MQ=255)
                                                                                                                                                                    ccTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                       >  1:353020/1‑250 (MQ=255)
                                                                                                                                                                       taTTACCCGGCACATAAATA‑CAAATGAAAATGGGTAGGGGGATAAAAACGAGGCTGGCA‑GAT‑CC‑CCCCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGCCCAACCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCTTCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACAcc                                                                    <  1:103943/249‑1 (MQ=255)
                                                                                                                                                                              cGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCAGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGGCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTCACAt                                                            >  2:159535/1‑251 (MQ=255)
                                                                                                                                                                               ggCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GACGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                           >  2:329381/1‑251 (MQ=255)
                                                                                                                                                                                   cacaAACA‑GCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCACCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                        >  1:325066/1‑250 (MQ=255)
                                                                                                                                                                                      aaaaa‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGACATTAGCAGCATCACTGCCGCCGTTTTTACCGCCACCTGCAATCCACCTTTGACATCGTAAccc                                                    >  2:313396/5‑251 (MQ=255)
                                                                                                                                                                                        aCA‑GAAATTAAATGGGTGATGGGGATAAAAATCGGGCTTGCT‑GTTCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                                  <  1:515663/251‑1 (MQ=255)
                                                                                                                                                                                                   aaatgggtggtggataaaaaaTCAGGCTGGCC‑AAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                        <  1:521562/235‑1 (MQ=255)
                                                                                                                                                                                                   aaaTGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCg                                        >  2:71321/1‑251 (MQ=255)
                                                                                                                                                                                                   aaaTGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                        >  1:539924/1‑251 (MQ=255)
                                                                                                                                                                                                              gggATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc                             >  2:309549/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaTCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAACACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca                      >  1:466034/1‑250 (MQ=255)
                                                                                                                                                                                                                              ttGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGCCATCGTAACCCGCGAGGCTGCCGCTAATCAGCCCCATCACAAAACCCACg             >  2:119863/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggtttgccaacccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc                <  1:171390/238‑1 (MQ=255)
                                                                                                                                                                                                                              gcctgccagatccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc                <  1:414030/238‑1 (MQ=255)
                                                                                                                                                                                                                                           caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt   >  2:300799/3‑251 (MQ=255)
                                                                                                                                                                                                                                            cccGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:309549/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                        
CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGA‑TCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACA‑GCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173130‑2173612

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