Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.85.3% 102.1 / 15.6 34intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (20/9);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 7.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01
*NC_0009132,173,3640C.85.3% 104.9 / 12.9 34intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (20/9);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 7.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173130‑2173612
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                        
ccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                        >  1:24037/1‑251 (MQ=255)
 cgccAGTTGGGTGTGCAAACCACTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGCGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagccgc                                                                                                                                                                                                                                        >  2:267629/1‑247 (MQ=255)
                      aaTAGTTTGCGTCGCGCTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                   >  1:426232/1‑250 (MQ=255)
                                  cgcgATCCAAAGGGTGATGATCATAATCATGACACCCGAGATTAAGGTGCGGAACAGACTTCCACGATGCACGGCGACCGCCCTCGCAACGAAGAAGCCGATCGTGGCAAGATAGCCGAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGACCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCGGAc                                                                                                                                                                                                      >  2:390457/1‑251 (MQ=255)
                                         cacaGGGTGATGCTCATAATGATGACAACCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGGTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                                                                                                                                                >  2:573146/1‑250 (MQ=255)
                                          acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                              >  2:144282/1‑251 (MQ=255)
                                              ggtaaTGCTCAAAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTCCCCCGATGCCCGGCCACCCCCATGCCCCCGAAGAAGCCGATGGTGGCAAGATCGCAAACCGGAAGCTCCGGATTCCCCGGCACACAAACAGCAATTAAAAGGGTGGGTGGAATAAAAACCAGGCTTGCCGTTCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                          <  2:490243/247‑1 (MQ=255)
                                                               tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                          >  1:82340/1‑250 (MQ=255)
                                                               tgaCACACGAGCTTAAGGGGCGGCACAGATTTCCACGCTGCACGGCGACAGCAATCGCAGCGAAGAAGCCGATGGTGGAAAGATCGACAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACAACC‑‑GCCGTAGGTCCCAGCAGCAAGGCAGGATCAAGGCCAATCCGGAACTCATGAACGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                           >  2:94205/1‑249 (MQ=255)
                                                                     ccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGc                                                                                                                                                                   >  1:98061/1‑251 (MQ=255)
                                                                               ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                                                            >  2:217830/1‑200 (MQ=255)
                                                                               ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                                                            <  1:217830/200‑1 (MQ=255)
                                                                                   cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGt                                                                                                                                                     >  1:268513/1‑251 (MQ=255)
                                                                                    ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGtt                                                                                                                                                    >  1:42467/1‑251 (MQ=255)
                                                                                           aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTACAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                             >  2:27499/1‑251 (MQ=255)
                                                                                                  cgagccacgcccaCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGCCAGACCCAGATTCCCCGGCCCAAAAACGGCAATCAAAATGGGGAGTGGGATAAAAACAAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCCGCTGCAACGCCGGACTACGGCCATTCAGGAACTCCTGAACGCCGAACTTCGCAGGGAAAAGAAGAAGAGCTTGCTCAGCGATGGGCGTTAAACCATCCatgat                                                                                                                                      <  1:543546/239‑1 (MQ=255)
                                                                                                         cGGCCACCGCCATCGCCACGAAGATGCCGATGGTGGTAAGATCGCCAAACTGCATCCCCTGATTCCCCGCCACACAAACATCAATTAAAATGGGGAGTGGGATAAAATTCAGGCTTGCCTATACCACC‑‑GCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAATCGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa                                                                                                                               <  1:358132/251‑1 (MQ=255)
                                                                                                           gcgccggCCACCGCCCCGAAGAAGCCGAGGGTGGCAAGACCGCCAAACGGCAGCACCTGATAACCGGGCACACAAACAGCAATAAAAATGTGGAGTGGAAAAAAAATCGGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGATTTAAt                                                                                                                             <  2:105767/245‑1 (MQ=255)
                                                                                                            cccccgccaccgccaCGAAGAAGCCGTGGTTGCAAAGAACGCAAACCGGCAGCACCGGATTACCCGGACCCAAAACAGCAATAAAAAGGGGGAGGGGGATAAAAATCAGGTTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                            <  2:310936/248‑1 (MQ=255)
                                                                                                                            gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCACCGGTACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCATTagc                                                                                                             >  1:235790/1‑250 (MQ=255)
                                                                                                                                   ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGAAGCATCAc                                                                                                     >  2:281410/1‑251 (MQ=255)
                                                                                                                                                    cGCCAAGCGGCAGCACCTGATCACCCGGCACAAAAACACCAATAAAAAGGGTGAGGGGGATAAAAACCAGGCTTGCCGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                     <  2:42467/251‑1 (MQ=255)
                                                                                                                                                       cAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                  <  1:139087/251‑1 (MQ=255)
                                                                                                                                                              cagccCCTGATTACCCGGCACACAAACACCAATTAAAATGGGGAGGGGGATAAAAACCAGGCTGGCCGATAC‑‑CGCCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                           <  1:518367/251‑1 (MQ=255)
                                                                                                                                                                     tGATCACCCGGCACAAAAACACAAATAAAAATGGGAAGGGGGATAAAAACCGGGCTTGCAGATAC‑‑CCCCACCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGACTACGAGCTTGCGTAGCGATGGGCGTTAAACCATCCATGATGGGGTTAATCACCCGGGTCATTCGCAGCAGCACTGCCGCCGTTTTGACCGCCAGTTGCAATacac                                                                     <  2:391132/250‑1 (MQ=255)
                                                                                                                                                                       aTTACCCGGCACACAAACAGCAATTAAAAGGGTGAGGGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCtt                                                                  <  1:144282/251‑1 (MQ=255)
                                                                                                                                                                                  cccgatctagaaaTTAAAAGGGTGAGGGGGATAAAAATCAGGCGTCCGGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                         <  2:454948/239‑1 (MQ=255)
                                                                                                                                                                                          aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                         >  1:454948/1‑242 (MQ=255)
                                                                                                                                                                                                  aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTGTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAATCCTGAACGCCGAACTTCGCCTCTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATACCg                                        >  1:555137/1‑251 (MQ=255)
                                                                                                                                                                                                    agggTGAGTGGAATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                      <  1:236863/249‑1 (MQ=255)
                                                                                                                                                                                                                        cAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCt                                                                                                                                                                  >  1:341178/1‑107 (MQ=255)
                                                                                                                                                                                                                        cAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAc                  >  1:480034/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCa               <  1:246513/251‑1 (MQ=255)
                                                                                                                                                                                                                             ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg             >  2:187373/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ccgatacccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCATCCCCATCACAAAACCCACGGTg          <  1:143640/244‑1 (MQ=255)
                                                                                                                                                                                                                                     atcccaccgccgTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAGGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACa       <  1:69484/245‑1 (MQ=255)
                                                                                                                                                                                                                                     ataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGCGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTgcc        >  2:234705/7‑248 (MQ=255)
                                                                                                                                                                                                                                         caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATAACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt   >  2:176393/3‑251 (MQ=255)
                                                                                                                                                                                                                                          cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:467157/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                        
CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173130‑2173612

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: