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.91.1% 160.8 / 7.7 45intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base . (29/12);  total (31/14)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
*NC_0009132,173,3640C.93.3% 168.3 / 4.4 45intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (30/12);  total (31/14)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGG  >  NC_000913/2173130‑2173609
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                      
ccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCag‑cagc                                                                                                                                                                                                                                     >  2:219415/1‑251 (MQ=255)
        tGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCgg                                                                                                                                                                                                                             >  1:622733/1‑251 (MQ=255)
             gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCa                                                                                                                                                                                                                         >  1:545653/1‑250 (MQ=255)
             gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCa                                                                                                                                                                                                                         >  1:659475/1‑250 (MQ=255)
                caaacaaaTAGTTTGCGTCGCGATCAACAGGGTGATGCTCAAAATGATGACCCCCGAGATTAAGGGGCGGAACAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                     <  1:438001/251‑1 (MQ=255)
                 aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                    <  1:580850/251‑1 (MQ=255)
                      aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAAAAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCa                                                                                                                                                                                                               <  1:506708/251‑1 (MQ=255)
                         aGTTTGCGCCGCGATCCACAGGGTGATGCTCATAATGAGGACCCCCGAGATTAAGGTGCGGACCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                                                                                                                                            <  1:103597/251‑1 (MQ=255)
                          gTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                                                           <  1:87886/251‑1 (MQ=255)
                                  cgcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCATGAc                                                                                                                                                                                                   >  1:232478/1‑251 (MQ=255)
                                            aGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCACGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCACCGCCGGAACAAGGCCAATCAGGGACTCCTGACCGCCGCACtt                                                                                                                                                                                         >  2:12577/1‑251 (MQ=255)
                                                   tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGCTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                      >  1:336252/1‑231 (MQ=255)
                                                      tcatcaTGATGACACCCGCGATTACGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGACGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACCAACAGCAATTAAAATGGTGAATGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGACACGCCGGATCAAGGCCAATCCGGAACTCCTAACCGCCGACCTTCGCCTGTaa                                                                                                                                                                                >  2:105954/1‑250 (MQ=255)
                                                           atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                                   >  1:197643/1‑242 (MQ=255)
                                                           atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                          >  2:594335/1‑251 (MQ=255)
                                                            tgatgaAACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAACCGGAACACAAACCGCAATTAAAATGGAGAGTGGGATAAAACTCAGGCTTGCCGATACCACC‑‑GCCGTGTGTCCCGT‑CAGCAACGCCGGCTCAAGGCCATTCAGGACATCCAGACCGCCGCACTTCGCCTGCAAACGAcg                                                                                                                                                                          >  2:49128/1‑249 (MQ=255)
                                                                     ccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCAAGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGTTAC‑ACCC‑TCCGTATGTCCCAA‑CTGAAAAGCCGGACCAAGGCCAATCAGAAACTCCGGACCGCCCAACTACGcatggaaaagacgaagagattgc                                                                                                                                                                >  2:473556/1‑229 (MQ=255)
                                                                        agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                              >  1:630346/1‑250 (MQ=255)
                                                                           ttAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCg                                                                                                                                                          >  1:396484/1‑251 (MQ=255)
                                                                                gggcgaaaCAGATTTCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATggg                                                                                                                                                     <  1:17920/245‑1 (MQ=255)
                                                                                         aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑ACCGTATGTCCCAG‑CAGAAACGCCGGATCAAGGCCAATCAGGAACACCTGACCGCCGAACTTCGCCTGTAAACGCCTACTAACATGCTTAACGATGGGCGTTAAAcc                                                                                                                                            >  2:374757/1‑251 (MQ=255)
                                                                                               ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                      >  1:124139/1‑251 (MQ=255)
                                                                                                        aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                             >  1:567982/1‑251 (MQ=255)
                                                                                                          gccGACCGCCACCGCCACGAAGAAGCCGAGGGGGGCAAGATCGCCAAACGGTAGCACCTGATTACCCGGCACAAAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGCCTTGCCGATCCCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                                           <  1:524881/249‑1 (MQ=255)
                                                                                                               ccGCCACCGCCCCAAAGAAGCCGATGGTGGAAAGAACGCCAAACGGAGGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                      <  1:320052/251‑1 (MQ=255)
                                                                                                                             aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTGCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                          >  1:374425/1‑249 (MQ=255)
                                                                                                                              agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTCCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagaagc                                                                                                       >  2:567921/1‑251 (MQ=255)
                                                                                                                                   ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                  >  1:553996/1‑251 (MQ=255)
                                                                                                                                       tggtggCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGGTACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                              >  1:536743/1‑251 (MQ=255)
                                                                                                                                          gggaaaGATCGCCAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                           <  1:219415/247‑1 (MQ=255)
                                                                                                                                             cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                        >  1:527392/1‑251 (MQ=255)
                                                                                                                                                            gcacaccttgtttccccGGCACACAAACACAAATAAAAAGGGTGATGGGGAAAAAAATCAGGCTGTCGGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTGAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc                                                                         <  1:344005/237‑1 (MQ=255)
                                                                                                                                                                            cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                         >  1:267125/1‑251 (MQ=255)
                                                                                                                                                                            cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTACCGATACCAAC‑‑GAAGTATGTACCAG‑CAGCAACGCCGGATCAAGGCCCATCAGGCACTCCTAACCGCCGAACTTCGCCTGTAAACGACTAAGCTCTTACTAAGCGATGGGCTTTACACCATCCATGATGGGTTTAATCAACCGTGTCATTAGCATCATCTCTGCCGCCGTTTATACCGCCAGCCGCAATACAACtttaaccc                                                         >  2:213239/1‑246 (MQ=255)
                                                                                                                                                                              gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                       >  1:49253/1‑251 (MQ=255)
                                                                                                                                                                                   aaaccaccAATAAAAAGGGTGAGTGGAATAAAAACCAGGCTTGCCAATCCCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                  <  1:566180/244‑1 (MQ=255)
                                                                                                                                                                                    aaaCAGCAATTACAATGGTGAGAGGGATAAAAATCACGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCTTGCTCAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATAAGCAGCATCACTCCCGCCGTTTTTACAGCCAGCTGCCATACACCTTTGACATCGTAAccc                                                 >  2:168179/1‑251 (MQ=255)
                                                                                                                                                                                       cAGCAATTACAATGGTGAGTGGGATCCAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAAAGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGCTGGTCGTTAAACCATCCATGACGGGTTTAATCACCCTGGGCATTAGCCGAATCACTGCCGCCGTATTTACCGCAAACTGCCATACACATTTGACATCGTAACAcgcg                                              >  2:656878/1‑251 (MQ=255)
                                                                                                                                                                                          tacttAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTATACCATCCATGCTGGGTTTAATCACCCGCGGCATTAGCAGCATCCCTGCCGCCGTTTTTCCCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcgcgg                                           >  2:27663/4‑248 (MQ=255)
                                                                                                                                                                                               aaaaGGGTGAGTGGAATAAAAACCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                      <  1:439507/251‑1 (MQ=255)
                                                                                                                                                                                                       gAGTGGTATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGACCCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCa                              <  1:12577/251‑1 (MQ=255)
                                                                                                                                                                                                               aaaaaacaggggtgcgcacaaccaccCC‑CCGTTTGTCCCAGACA‑CACCGCCGGAGTAAGGCCAATAAGGAACCCCGGACCTCCGAACTTCGCCTTTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGTTGTGTTTAATCACCAGGTGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATCAGCCCCAt                       <  2:195678/231‑1 (MQ=255)
                                                                                                                                                                                                                    aTCAGGCTTGCCGCTA‑‑ACACCGCCGAATGTCCCCA‑CAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAAGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCCTCCATGATGGGTTTAATCACCCGTGGCATCAGCAGCATCCCAGCCGCCGTTTTTACCGCCAGCTGCAATAACCCTTTGACATCGTAACCCGCGCAGATGCCGCTAATCAGCCCACTCACaa                  >  2:310437/1‑250 (MQ=255)
                                                                                                                                                                                                                     tccggCTCGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCCGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATAAGCCCCCTCACaaa                 >  2:100421/4‑250 (MQ=255)
                                                                                                                                                                                                                       aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAcc              >  1:367696/1‑251 (MQ=255)
                                                                                                                                                                                                                              gcccagtacagccgccgTATGCCCCAG‑AATCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCGCCTGTAAACGACTACGATCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATCAGCCCCATCTCAAAACCCAc           <  1:374757/238‑1 (MQ=255)
                                                                                                                                                                                                                                   ataccaccgccgTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                      >  2:271328/7‑121 (MQ=255)
                                                                                                                                                                                                                                   ataccaccgccgTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                      <  1:271328/115‑1 (MQ=255)
                                                                                                                                                                                                                                     accaccgccgTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg  >  1:501981/5‑251 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                      
CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGG  >  NC_000913/2173130‑2173609

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