Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.92.5% 156.6 / 3.4 40pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (23/14);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.92.3% 149.5 / 3.5 40pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (23/13);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAA  >  NC_000913/2173131‑2173593
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                       
cgcgAGTTGGGTGTACAAACGAATAGTTGGGGTGGGAGTACACGGTGTGATGCTAAAAATGAGGACACCCGAGATTAAGGGGGGAAACAGTTTCCCAGGATGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCAAAACGGCAGCACCTGATTC‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATaccaccgccgtat                                                                                                                                                                                                                                  <  2:244565/239‑8 (MQ=255)
   cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGCTACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                              >  2:226221/1‑238 (MQ=255)
             tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                                                                                                                    >  1:317162/1‑238 (MQ=255)
                     aaTGGTTGGCGGCGCGATCCACAGGGTGATGCTCATAATGATGACCCCGGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCAAGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                           <  1:71597/239‑1 (MQ=255)
                           ttGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACTGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                      >  2:698336/1‑238 (MQ=255)
                                       ccACGGGGTGATGCTCAAAAGGATGCCCCCCGAAATTAGGGTGGGGAACAGTTTCCCACGAGGCCCGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATaccaccg                                                                                                                                                                                                                                        <  2:117132/194‑2 (MQ=255)
                                                 aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                                                                                                                            >  1:117132/1‑194 (MQ=255)
                                                  tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                           >  1:617003/1‑210 (MQ=255)
                                                              tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                       >  2:72284/1‑234 (MQ=255)
                                                                                tggggAAAAGATTTCCAGGATGCGCGGCGACCCCCATCGCCACGAAGACGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTC‑CCCGGCCCACAAACACCAATAAAAATGGTGGGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                       <  2:40264/229‑1 (MQ=255)
                                                                                  cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATaccac                                                                                                                                                                                                                                          >  2:365806/1‑149 (MQ=255)
                                                                                  cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATaccac                                                                                                                                                                                                                                          <  1:365806/149‑1 (MQ=255)
                                                                                          aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                             >  1:40264/1‑232 (MQ=255)
                                                                                           gtcccACGGAGCGCGGCGCCCGCAGTCGCCACAGAGAACCGGGAGGGTGCAGGATCGACAAGCTGCA‑CACCTCGTTAGCCGGGCAGACAACCAGCATATAAAGTGTGGAGTGGTATAAAACCCCGCCTG‑GCGGA‑ACCA‑ACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACCCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCGGCTTAGCGAt                                                                                                                                        <  2:7039/233‑1 (MQ=255)
                                                                                              ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                                                                                                                       >  1:76398/1‑170 (MQ=255)
                                                                                                 cGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATGA‑CCCGGCATACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                       >  1:849709/1‑215 (MQ=255)
                                                                                                             cccccattccccAGGAAGAGCCGGATGGTGGCAAGACCGCAAACCGGAGGCCCTTGTTAC‑CCCGGCA‑ACAAACAAGCATTAAAAAGGTGGGGGGGAAAAAAAATAGGGCT‑TCCGGTTCC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                   <  2:156647/229‑1 (MQ=255)
                                                                                                              ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAAATTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                  >  2:583596/1‑239 (MQ=255)
                                                                                                              ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                  >  2:84335/1‑239 (MQ=255)
                                                                                                                      accggaaaaaagCCGAGGGTGGAAAAACCCCCAAACGGCGCCCCCTGATTC‑CCCGGCACACAAACAGCAA‑TAAAATAGTGGGGGGGAAAAAAAACCGGCTTGGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                          <  2:316544/230‑1 (MQ=255)
                                                                                                                                     aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                >  1:316544/1‑234 (MQ=255)
                                                                                                                                      tggtggCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                          >  1:211631/1‑239 (MQ=255)
                                                                                                                                        gtggCAAGATCGCCAAACGGCAGCACCTGATAA‑CCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                      <  2:837423/129‑1 (MQ=255)
                                                                                                                                                  cccaaaaCGGCAGCCCCTGATTA‑CCCGGCACACAACCAGCAATTAAAAGGGTGGGTGGAAAAAAAACCGGCCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                     <  2:432881/116‑1 (MQ=255)
                                                                                                                                                  cGCCAAACGGCTGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                                                                                                                            >  1:837423/1‑129 (MQ=255)
                                                                                                                                                   gCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                                                                             <  1:84335/239‑1 (MQ=255)
                                                                                                                                                            cagcaCCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                                           >  1:432881/1‑120 (MQ=255)
                                                                                                                                                               caCCTGATTA‑CCCGGCACACAAACAGCAATAAAATTGGTGAGTGGGATAAAAAGCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                                      <  2:344632/90‑1 (MQ=255)
                                                                                                                                                                    gATTC‑CCCGACACAAAAACAGCAATTAAAATGGGGAGGGGAAAAAAAACCAGGCTT‑GCGAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                            <  1:544246/239‑1 (MQ=255)
                                                                                                                                                                          cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                                                                                                            >  1:344632/1‑90 (MQ=255)
                                                                                                                                                                           ccGGCACCAAAACAGAAATTAAAATGGTGGGTGGGATAAAAAGCCGGCTT‑GCCGGTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                      <  1:226221/239‑1 (MQ=255)
                                                                                                                                                                               cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATCCCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                                      <  2:109038/75‑1 (MQ=255)
                                                                                                                                                                               cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                                      >  1:109038/1‑75 (MQ=255)
                                                                                                                                                                                          cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGTCGCCGTTTTTACTGCCAGCTGCACTACACCCTTGACATCg                                       >  2:758589/1‑239 (MQ=255)
                                                                                                                                                                                                      tGAGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                           >  1:182370/1‑239 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                          >  2:248538/1‑142 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAATCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                          <  1:248538/142‑1 (MQ=255)
                                                                                                                                                                                                             gaaaaaaaTCAGGCTT‑CCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                     <  1:319878/235‑1 (MQ=255)
                                                                                                                                                                                                                aaaaaTCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                >  2:3435/1‑96 (MQ=255)
                                                                                                                                                                                                                aaaaaTCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                <  1:3435/96‑1 (MQ=255)
                                                                                                                                                                                                                  aaaTCAGGCTT‑GCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc               <  1:698336/239‑1 (MQ=255)
                                                                                                                                                                                                                        ggCTG‑GGCGATACC‑CCCCGCCTTATGCCCCAACAGCAACGCCGGATCAAGCCCAATGAGGACCCCCGGCCCGCCGAACTTCGCCTTTAAACGACTGCGAGCCTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCGCCCGTGGCATTGGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACGTCGTAACCCGCAAGGATGCCGATAATCAGccc          <  2:211631/239‑1 (MQ=255)
                                                                                                                                                                                                                              gCCGATACCA‑‑CTGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                                                                             >  1:824006/1‑166 (MQ=255)
                                                                                                                                                                                                                              gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                              >  1:704004/1‑117 (MQ=255)
                                                                                                                                                                                                                                  gaatccccccGCCGTCTGACCAGGCAACAGCGCGGGAGCAAGGCCAACAAGAACCTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACTAGCCTGCTTAACGGTGGACGTTAAACCATCCATGATGGGTTTAATCCCCCATAACAATAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaa  <  2:182370/233‑1 (MQ=255)
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                       
CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTT‑GCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAA  >  NC_000913/2173131‑2173593

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