Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 860,154 A→T 7.5% intergenic (‑126/+20) ybiV ← / ← ybiW sugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,1540AT 7.5% 256.7 / 10.8 80intergenic (‑126/+20)ybiV/ybiWsugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base T (3/3);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-03

TGATTACCGCTGGCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATT  >  NC_000913/859889‑860425
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                               
tGATTACCGCTGGCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACNTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGCAGCAGGNCACAAATTATATTTCGAATGAAAGAAAGGATGAAAATGATGAAGGNANTGATGTAGCCCGGCGACGACGCCACCgg                                                                                                                                                                                                                                                                        >  3:33548/1‑275 (MQ=255)
 gATTACCGCTGGCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCggg                                                                                                                                                                                                                                                                       >  2:286097/1‑275 (MQ=255)
   ttACCGCTGGCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTGCTGGATCGAAAATGAAACCGGAACAGTGTAATAACAATGTGACGCGGAGCAAAAATTATATTTCGAATGAGAGTAAGGCTGAAATTGATGAGGTGAAGGATTTAGCCCGGCAACGACGCCGCCgggac                                                                                                                                                                                                                                                                     >  2:351575/1‑273 (MQ=255)
   ttACCGCTGGCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCACCGAGcc                                                                                                                                                                                                                                                                     >  2:245260/1‑275 (MQ=255)
           ggCAACAACGAACTTAAGGCCGTGCTTTTTAAGTTCCTGTTATTTCGCCATAAAACGTGGTTGGTTGGACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATAAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAt                                                                                                                                                                                                                                                             <  2:506272/275‑1 (MQ=255)
            gCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                                                                                                                                                                                                            >  2:420326/1‑275 (MQ=255)
                   ggAAGTTAATGCCGGGCTTGTTCAGTTCCTGATATTGAGCCAGAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTGAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTAGATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCat                                                                                                                                                                                                                                                     <  4:97655/274‑1 (MQ=255)
                   caaaCTGAATGCCGCGCGTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGGCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCat                                                                                                                                                                                                                                                     <  4:107621/273‑1 (MQ=255)
                             gCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTTGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATATAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTa                                                                                                                                                                                                                                           <  2:342374/275‑1 (MQ=255)
                               cgcgCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACg                                                                                                                                                                                                                                         <  4:401095/275‑1 (MQ=255)
                                        tCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatgat                                                                                                                                                                                                                                <  2:317543/275‑1 (MQ=255)
                                                 gATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTAc                                                                                                                                                                                                                                                          <  2:2083/240‑1 (MQ=255)
                                                 gATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTg                                                                                                                                                                                                                       <  4:420325/275‑1 (MQ=255)
                                                        cgcCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGATGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGACATTGATGATGTGACTGATTTAGCACGGCGACGAAGCCGCCGGGCCGAGGAGATTACAGCATATGTGCAATACGGGCGATGATATAGTCCTGAGCGTCt                                                                                                                                                                                                                >  1:6641/1‑275 (MQ=255)
                                                          caaTAAAACGTTGTTGTTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTgg                                                                                                                                                                                                              <  4:163576/273‑1 (MQ=255)
                                                             tAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgaag                                                                                                                                                                                                           >  2:532342/1‑273 (MQ=255)
                                                                                    ttGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAACATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGCCCGAGGAGATTACGGCATAGATTCAGTACGGGCGATGAAATCGGTCTGAGCGTCGGGAAGGAGGGTGGTGAAAGACGCGGAAt                                                                                                                                                                                    >  4:380803/1‑275 (MQ=255)
                                                                                     tGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTg                                                                                                                                                                                                                       <  4:460203/239‑1 (MQ=255)
                                                                                     tGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTg                                                                                                                                                                                                                       >  2:460204/1‑239 (MQ=255)
                                                                                        cgtcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGATGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAAAGAGCGCGGTAAAGAACGCGGAAAAGc                                                                                                                                                                                 >  4:56390/1‑274 (MQ=255)
                                                                                        cgtcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATAGCGTCCTGAGCGTCTGGAGAGGGCGCGGTGAAGGACGGGGaaacgc                                                                                                                                                                                 >  4:127597/1‑270 (MQ=255)
                                                                                        cgtcgtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGGGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCTAATGAAAGTAAGGATGAAATTGATGATGTGAATGACTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAAc                                                                                                                                                                                           <  2:380804/264‑1 (MQ=255)
                                                                                            gtTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc                                                                                                                                                                                                                                                                      >  4:178385/1‑185 (MQ=255)
                                                                                             ttAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTGACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTAACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGa                                                                                                                                                                           <  1:33548/275‑1 (MQ=255)
                                                                                                aGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgag                                                                                                                                                                                                                                                                  <  2:178386/185‑1 (MQ=255)
                                                                                                aGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGaca                                                                                                                                                                         <  4:519010/274‑1 (MQ=255)
                                                                                                 gAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACg                                                                                                                                                                       <  2:160978/275‑1 (MQ=255)
                                                                                                       tACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGgagaga                                                                                                                                                                                                         >  2:89115/1‑235 (MQ=255)
                                                                                                         ccGTCCATGTCTGTGACGATAACTTTTACGCTAAGAGGTGTGCTCCTGGCTAGAAAATGAACCCTTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCa                                                                                                                                                               <  1:133761/275‑1 (MQ=255)
                                                                                                           gTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAgcgc                                                                                                                                                                                                     <  4:89114/235‑1 (MQ=255)
                                                                                                                            tAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCgg                                                                                                                                             >  4:456730/1‑274 (MQ=255)
                                                                                                                                          aTAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGa                                                                                                                                                                           >  2:284350/1‑230 (MQ=255)
                                                                                                                                              gtgtGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATGGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACg                                                                                                                                                                       <  4:284349/230‑1 (MQ=255)
                                                                                                                                                 tgCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGTCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGACCGCGGAATAGCCCGCGACCCGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTTGACTCcagca                                                                                                                        >  1:106583/1‑274 (MQ=255)
                                                                                                                                                  gCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAg                                                                                                                                                              >  2:420691/1‑235 (MQ=255)
                                                                                                                                                      cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATcg                                                                                                                                                          <  4:420690/235‑1 (MQ=255)
                                                                                                                                                        ggCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTcg                                                                                                                >  2:66822/1‑275 (MQ=255)
                                                                                                                                                                aaTGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGaaa                                                                                                         <  2:183822/274‑1 (MQ=255)
                                                                                                                                                                    aaaCCGTAACAGTGTAATAACAATGTGACGCAGAGAACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAtg                                                                                                    <  4:286096/275‑1 (MQ=255)
                                                                                                                                                                         gTAACAGTGTAATAACAATGTGACGCAGAGCAAAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTTAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTa                                                                                               <  2:224551/275‑1 (MQ=255)
                                                                                                                                                                                tgtAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGGAGAACGCGGACTAGCCCGCGACACGCACTACCAGATCGCGGTACTGAGCGGGATGCTTTTTCGCATCCAGCAGCGTTTAGAGGGAAACGAAGTTGTACTGAat                                                                                         >  2:383339/1‑274 (MQ=255)
                                                                                                                                                                                  taataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGGGGGAAACGATGTTTTACTGAattt                                                                                       >  4:81341/1‑272 (MQ=255)
                                                                                                                                                                                   aataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCata                                                                                                                                                                                                                                                    <  5:985141/116‑1 (MQ=255)
                                                                                                                                                                                   aataaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCata                                                                                                                                                                                                                                                    >  6:985141/1‑116 (MQ=255)
                                                                                                                                                                                        cAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGa                                                                                                                                                                           >  2:61081/1‑184 (MQ=255)
                                                                                                                                                                                          aTGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGAtt                                                                                                                                                                                                                                                            >  1:28645/1‑101 (MQ=255)
                                                                                                                                                                                            gtgACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACg                                                                                                                                                                       <  4:61080/184‑1 (MQ=255)
                                                                                                                                                                                              gACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAg                                                                                                                                                                                                                                                        <  3:28645/101‑1 (MQ=255)
                                                                                                                                                                                               aCGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATg                                                                          <  2:456731/274‑1 (MQ=255)
                                                                                                                                                                                                cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAGATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgaggag                                                                                                                                                                                                                                                               <  8:672197/92‑1 (MQ=255)
                                                                                                                                                                                                cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCgaggag                                                                                                                                                                                                                                                               >  7:672197/1‑92 (MQ=255)
                                                                                                                                                                                                cGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGCAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACCGGTCGGGATGTTTTTGCTCATCCACGAGCGTTTCGCGGGAACAGATTTTGTTATGACTATGCCAGCCCTTATGcc                                                                        >  4:445492/1‑274 (MQ=255)
                                                                                                                                                                                                            aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGtcgccg                                                                                                                                                                                                                                                                         >  7:583763/1‑70 (MQ=39)
                                                                                                                                                                                                            aTTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGtcgccg                                                                                                                                                                                                                                                                         <  8:583763/70‑1 (MQ=38)
                                                                                                                                                                                                                 atTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACCGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGCGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACga                                                       >  2:73490/1‑274 (MQ=255)
                                                                                                                                                                                                                      gatagacagtaagtaTGGAATTGATGCTGTGAATGATGTACCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGAGANCGTCCTTAGCGTCTGGAGAGAGCGCGGTGAAGCACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGcc                                                                                    <  3:106583/234‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                                                                                                                                                                                                           >  6:563521/1‑105 (MQ=255)
                                                                                                                                                                                                                       aaTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGt                                                                                                                                                                                                                           <  5:563521/105‑1 (MQ=255)
                                                                                                                                                                                                                            aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAaagaa                                                              >  2:423067/1‑257 (MQ=255)
                                                                                                                                                                                                                            aaGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGatcat                                            <  2:56391/275‑1 (MQ=255)
                                                                                                                                                                                                                                aaGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTa                                                          <  4:423066/257‑1 (MQ=255)
                                                                                                                                                                                                                                  ggATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGtt                                       >  4:224157/1‑274 (MQ=255)
                                                                                                                                                                                                                                    aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGt                                                                                                                                                                                                                                                 <  8:508323/70‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGt                                                                                                                                                                                                                                                 >  7:508323/1‑70 (MQ=255)
                                                                                                                                                                                                                                       aaaTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCCGTTCGCGGGAACCGATGTTGTACTTAATATGCCAGCCTTTATGCATTTCAAAGAAGGTACGCCGCAGGATCATCAGTTTcct                                   >  4:457510/1‑271 (MQ=255)
                                                                                                                                                                                                                                         aTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCgg                                                                                                             >  4:13190/1‑197 (MQ=255)
                                                                                                                                                                                                                                         aTTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCCGTTTCTGCtt                                >  2:409485/1‑274 (MQ=255)
                                                                                                                                                                                                                                          ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                                                                                                                                                                                                     <  7:132663/76‑1 (MQ=255)
                                                                                                                                                                                                                                          ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCg                                                                                                                                                                                                                                     >  8:132663/1‑76 (MQ=255)
                                                                                                                                                                                                                                          ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGGTATCGTCCTGAg                                                                                                                                                                                                                     >  7:369560/1‑92 (MQ=255)
                                                                                                                                                                                                                                          ttGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGGTATCGTCCTGAg                                                                                                                                                                                                                     <  8:369560/92‑1 (MQ=255)
                                                                                                                                                                                                                                             atgatgTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGaaa                                                                                                         <  2:13190/197‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgtgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGAtg                                                                                                    <  1:9672/198‑1 (MQ=255)
                                                                                                                                                                                                                                                   tgAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCg                      <  2:224158/274‑1 (MQ=255)
                                                                                                                                                                                                                                                        gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatgat                                                                                                                                                                                                                                >  7:516338/1‑67 (MQ=255)
                                                                                                                                                                                                                                                        gATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCgatgat                                                                                                                                                                                                                                <  8:516338/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                                                                                                                                                               <  6:528882/42‑1 (MQ=25)
                                                                                                                                                                                                                                                                  ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTc                                                                                                                                                                                                                                               >  5:528882/1‑42 (MQ=25)
                                                                                                                                                                                                                                                                  ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGa                                                                                                                                                             <  7:194859/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggcggcGTCGTCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGa                                                                                                                                                             >  8:194859/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                    cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATa                                                                                                                                                     >  7:1020296/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                    cgacgaCGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATa                                                                                                                                                     <  8:1020296/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cgacgCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGAtt  <  1:14665/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                               
TGATTACCGCTGGCAACAACGAACTTAATGCCGCGCTTTTTCAGTTCCTGATATTGCGCCATAAAACGTGGTTGGTTGTACGTTTTGGCGTCGTTAAGAAAAGTACCGTCCATGTCTGTGACGATAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCATCCAGCAGCGTTTCGCGGGAAACGATGTTGTACTGAATATGCCAGCCTTTATGCACTTCAAAGAAGGTACGCAGCAGGATCATCAGTTTCTGCTTGTCAGATTCGTTCTCCAGCGTTGCCGGATT  >  NC_000913/859889‑860425

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: