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.89.7% 140.9 / 8.7 39intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (27/8);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 6.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01
*NC_0009132,173,3640C.89.7% 141.6 / 8.0 39intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (27/8);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 6.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC  >  NC_000913/2173124‑2173600
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                            
cATTAGCCGCCAGTTGGGTGTGCAAACAAATAGTTTGCGTCGCGATCCACAGGGGGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGCACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGACGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAACTCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                  >  2:499688/1‑251 (MQ=255)
   tAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTGCCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGT‑TACCACC‑‑GCCGTATGTCCcag                                                                                                                                                                                                                                <  1:611782/250‑1 (MQ=255)
     gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGCGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCcatc                                                                                                                                                                                                                               >  2:495782/1‑247 (MQ=255)
                      cAAACCAATATTTTCCGTCGCGATCCACAGGGTGAGCCCCAAAATGAGGCCACCCGAGATTAAGGTGCGGAACAGATTCCCCCGATGCACGGCGTCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTCCCCGGCCCACAACCGGCAATAAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCAA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGTATCAAGGc                                                                                                                                                                                                            <  1:598370/251‑1 (MQ=255)
                          ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCTGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                        >  2:503812/1‑251 (MQ=255)
                               agttggggTCGCGATCCACAGGGGGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCTACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCATCACGTGATTACCCTGCGCACAAACAGCAAGTAAAA‑TTGTGAGTGGGATAAAATTCAGGCTTGC‑GAGTTCGACC‑‑TCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                                                                                                                                   <  1:439269/244‑1 (MQ=255)
                                      gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                            >  1:58262/1‑251 (MQ=255)
                                                                    atgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                              >  1:304538/1‑251 (MQ=255)
                                                                       acacCCGAGATTAAGGAGCGGAACAGATTTCCACGATGCACGGCGACCGCCATAGCCACGAAGAAGCCGATGGTGGCAAGATCGCCCAAC‑GGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATAAGGAACTCATGACCGCCGAACTTCGCCTGTAAACGCCTAcacac                                                                                                                                                           >  2:642438/1‑247 (MQ=255)
                                                                         ccctgagaTTAAGGTGCGAAACAGATTCCCTGGAGGCCCGGGCACCGCCAACGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAAG‑GGCGGCACCTGATTACCGGGCACAAAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGTCTGTAAACGACTACGAGCCt                                                                                                                                                         <  1:592894/247‑1 (MQ=255)
                                                                           ccGAGATTAAGGTGCGAAACAGATTTCCACGAGGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGAAAGACCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGc                                                                                                                                                       <  1:134155/251‑1 (MQ=255)
                                                                           acgcgaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACAGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAGCTTCGCCTGTAAACGACAACGAGACTGc                                                                                                                                                       >  2:519481/5‑251 (MQ=255)
                                                                            cGAGATTAAGGTGCGGACCAGCTTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGACGCCGATGGTGGCACGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCCTGCt                                                                                                                                                      >  2:274648/1‑251 (MQ=255)
                                                                                  tAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                 >  1:485704/1‑250 (MQ=255)
                                                                                  tAAGGAGCGGAAAAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGAGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                 >  2:659517/1‑250 (MQ=255)
                                                                                         cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                         >  1:366802/1‑251 (MQ=255)
                                                                                                     ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTACAccatccat                                                                                                                             >  2:391262/1‑251 (MQ=255)
                                                                                                           tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                       >  1:8449/1‑251 (MQ=255)
                                                                                                                ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                                                                                                                             >  2:277661/1‑208 (MQ=255)
                                                                                                                ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                                                                                                                             <  1:277661/208‑1 (MQ=255)
                                                                                                                 gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                >  1:335308/1‑156 (MQ=255)
                                                                                                                  agACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTCAGCGATGGGCGGTAACACATCAATGATGGGTTTactc                                                                                                                >  2:393795/2‑248 (MQ=255)
                                                                                                                    aCCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAAC‑GGCAGGACCTGATTACCCGGCACACAAACAGCAATTACAC‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                    >  2:550397/1‑149 (MQ=255)
                                                                                                                       gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                           >  1:609855/1‑251 (MQ=255)
                                                                                                                            cacaccGAAGAAGCCGA‑GGGGGCAAAAGCCCCAAACCGGAAGCACCTGATACCCCGGCACAAAAACAGCAAGTAAAA‑GGGTGAGGGGGAAAAAAATCAGGCTGCCCGA‑‑ACC‑CACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCa                                                                                                      <  2:154338/246‑1 (MQ=255)
                                                                                                                                  gaagaaGCCGATGGTGGCAAGACCGCCAAAC‑GGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAACCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagca                                                                                                <  1:311363/251‑1 (MQ=255)
                                                                                                                                   aagaagCCGATGGGGGCAAGACCGCCAAAC‑GGCAGCACCGGATTACCCGGCAAACAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcag                                                                                               <  1:68176/251‑1 (MQ=255)
                                                                                                                                     gaagCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                              >  1:65616/1‑250 (MQ=255)
                                                                                                                                          ggagggtgcaaaGACCGCAAAAC‑GGCAGCCCCGAATAACCCGGCACACAAACAGCAATAAAAA‑GGGGGGGTGGGAAAAAAACCGGGCTGGCCGA‑‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTCTAATCACCCGTGGCATTAGCAGCATCACt                                                                                        <  2:8449/241‑1 (MQ=255)
                                                                                                                                                     aGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                             >  1:184500/1‑251 (MQ=255)
                                                                                                                                                                        aCCTGATTACCCGGCACCCCAACAGCAATTAAAA‑TGGTGAGTGGGCTAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTCGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                           >  2:270634/1‑251 (MQ=255)
                                                                                                                                                                         ccTGATTAACAGGCACACAAACAGCAATTAAAA‑TGGAGCGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGAACCGCCGGATCAAGGCCAATCAGGAACTCCAGACAGCCGAACTTCGCCTGTAAAAGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATTGGGTTAATCACCCGAGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                           >  1:438972/1‑250 (MQ=255)
                                                                                                                                                                                          aaaaaaagtaTTAAAAGTGGTGAGGGGGAAAAAAATCAGGCTGCCCGA‑T‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCATCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACGCCTTTGACATCGTAAc                                          <  1:557890/242‑1 (MQ=255)
                                                                                                                                                                                            aaCAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcg                                       >  2:442846/1‑251 (MQ=255)
                                                                                                                                                                                                   aTTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                                >  1:538199/1‑251 (MQ=255)
                                                                                                                                                                                                     ggttgctcgcccggtctaaaaaaaCCAGGCTGGCCGA‑T‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAGCGCTTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTGTACCTCGAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg                               <  1:279793/232‑1 (MQ=255)
                                                                                                                                                                                                             gTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAAt                       >  1:62662/1‑251 (MQ=255)
                                                                                                                                                                                                              tGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACC‑‑GCCGAATGCACCGGCAGCAACGCCGGATCAAGCACAAACAGGACATCCTGACCGCCGGACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGACGCCGATAATc                      >  1:439940/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaTCAGGCTTGCCGA‑TACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca          >  1:156794/1‑250 (MQ=255)
                                                                                                                                                                                                                            aaaacaaggcggcccgatccccccGCCGTCGGTCCCAGCAGCGACCCCGGACAAAGGCAAACAAGGAACCCCTGCCCGCCGAACTTCGCCTGTAAACGACGACGAGCCTGCTTAGCGAGGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGGGGCATGAGCAGCATCACTGCCGCCGTTTTGACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac           <  2:62662/231‑1 (MQ=255)
                                                                                                                                                                                                                                       ggccgatccccccGCCGTAGGCCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACCCCTTGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCTTGGCATTGGCAGCATCACTTCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc  <  1:107641/240‑1 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                            
CATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑TACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC  >  NC_000913/2173124‑2173600

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: ‑