Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,263,441 Δ1 bp coding (55/1131 nt) fruB ← fructose‑specific PTS multiphosphoryl transfer protein FruB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,263,4410C.98.5% 293.2 / ‑0.3 66coding (55/1131 nt)fruBfructose‑specific PTS multiphosphoryl transfer protein FruB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (28/36);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTT  >  NC_000913/2263195‑2263684
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cagcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGATCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGt                                                                                                                                                                                                                                                <  1:416872/251‑1 (MQ=255)
     aTACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCg                                                                                                                                                                                                                                           >  1:487356/1‑251 (MQ=255)
          gATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCtt                                                                                                                                                                                                                                      <  1:287102/251‑1 (MQ=255)
              gCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTc                                                                                                                                                                                                                                   <  1:9939/250‑1 (MQ=255)
               ccACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGc                                                                                                                                                                                                                                 >  1:398969/1‑251 (MQ=255)
                     aCGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCggat                                                                                                                                                                                                                           >  1:166485/1‑251 (MQ=255)
                       ccgcccTGACCGTAACCCCAGGTGACGCCGTCCGGGAACTGGAATACCTGAACGCCGGGTTTCAGCACCTGATCGCGGGTTGCGGTAGTGCTGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCAGGCCATTGACGTAGCCGTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCGTTTCGCCCggatgg                                                                                                                                                                                                                         <  2:663083/247‑1 (MQ=255)
                       ccTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTGTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCggatgg                                                                                                                                                                                                                         <  2:3982/250‑1 (MQ=255)
                            cTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTc                                                                                                                                                                                                                    >  1:305136/1‑251 (MQ=255)
                               cccGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCGGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATTCCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑GTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTg                                                                                                                                                                                                                 <  2:349655/250‑1 (MQ=255)
                                 cGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTa                                                                                                                                                                                                               <  1:861164/251‑1 (MQ=255)
                                  gTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTAc                                                                                                                                                                                                              >  1:121800/1‑251 (MQ=255)
                                         cAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAAc                                                                                                                                                                                                       <  1:205536/251‑1 (MQ=255)
                                          aGGTGCCGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCGGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACt                                                                                                                                                                                                      <  2:19955/251‑1 (MQ=255)
                                          aGGTGACGCCTGCCGGGAACTGAAATTCCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGGTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACt                                                                                                                                                                                                      <  2:852868/251‑1 (MQ=255)
                                                 gacgtcagtgcactgtaatcccttgaCGCCGGTTTTCAGCACCGGATCGCGGGTGGCGGTAGTGCCGTGTGGAATAGCAATTCCCGGGCAGAGTAACGTTTAGATTGGTTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTGCTGCTACATTCCCGGCCTGCACCTGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCGCGGATGGATGTCCTGTACGGATTACTGTAACAt                                                                                                                                                                                               <  2:152293/226‑1 (MQ=255)
                                                     tCCGGGAACTGACATACCTGAACGCCGGTTTTAAGCGCCAGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACACTTGCCGAGGAACGTTGAGGTTTACTGTTCGCGCGCCAGCATGCCATTGGCGTAGCCTTCTGCTACATTACCGGCCTGCACCCGCGCCGCAGAGACCCTGCGAATCGCCTCCT‑TTTGTATCCGGCCTTTTCGCCCGGGTGAAGGTCCTGTACGGATAACTGGAACATAGtt                                                                                                                                                                                           >  2:613854/1‑251 (MQ=255)
                                                      ccGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTtc                                                                                                                                                                                          >  1:349702/1‑251 (MQ=255)
                                                      cagggAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGACGCAGCGACCTGGCGCAGCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGCTGTCCTGTACGGATAACTGGAACATAGttt                                                                                                                                                                                          >  2:285002/3‑250 (MQ=255)
                                                           aattGAAATACGTGAACGCCGGTGTTCAGCAATGGATCGCGGGGGTCGTTAGGGTAGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGACAGCAGGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCCCCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGACTTGTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTctcctc                                                                                                                                                                                     <  2:191/248‑1 (MQ=255)
                                                                aaaTACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTg                                                                                                                                                                                <  1:345392/251‑1 (MQ=255)
                                                                     ccTGAACGCCGGTTTTCAGCACCTGATCTCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGa                                                                                                                                                                            >  1:534791/1‑250 (MQ=255)
                                                                        gcaCGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGACCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGCTGGATGTCCTGTACGGCTAAATGGACCATAGTTCTCCTCTCTTGATGAATTg                                                                                                                                                                        >  2:692692/3‑251 (MQ=255)
                                                                        gAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTg                                                                                                                                                                        >  1:870389/1‑251 (MQ=255)
                                                                        gAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTg                                                                                                                                                                        <  1:47271/251‑1 (MQ=255)
                                                                                ggttgcAGCTCCTGATCGCGGGTGTCGGTAGTGTCGTGTGGAATAGCAATGCCATTGCCGAGGAACGTTGAGGTTTGCTGGTCGCGCGCCAGCATGCCATTGCCGTAGCCTTCTGCTACATTGCCGGCCTGCACCAGCGCCGCATCGTCCTGGCGATTCGCCTCTT‑TTTGTCTCCGGCCTTGTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTTAAACGAtt                                                                                                                                                                <  2:5930/246‑1 (MQ=255)
                                                                                  tttcagacccTGATCGCGGTGGTCGGGAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTc                                                                                                                                                               <  2:487356/242‑1 (MQ=255)
                                                                                  tttCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTc                                                                                                                                                               <  1:692692/250‑1 (MQ=255)
                                                                                             gATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGATTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCctctctg                                                                                                                                                     >  1:524866/1‑246 (MQ=255)
                                                                                               tCGCGGGTGTAGGTAGTGACGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGGTAACTGGAACATCGTTCTCCTCTCTTGCTGAATTGACAAGATTCAGCCTCTATgtaa                                                                                                                                                  >  2:408998/1‑247 (MQ=255)
                                                                                                cgcgGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGaaa                                                                                                                                                <  1:842856/251‑1 (MQ=255)
                                                                                                   gggTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGaaaaa                                                                                                                                              >  1:485942/1‑250 (MQ=255)
                                                                                                              gTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCt                                                                                                                                  >  1:29489/1‑251 (MQ=255)
                                                                                                               ggCCGTGTGGAATAGCCATACCAGTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTg                                                                                                                                 <  2:236682/250‑1 (MQ=255)
                                                                                                                 ccGTGTGGAATAGCAATACCATTGCCGAGGACCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTCCATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGc                                                                                                                               <  2:485942/251‑1 (MQ=255)
                                                                                                                      tgGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGAAGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAgg                                                                                                                          >  2:323674/1‑251 (MQ=255)
                                                                                                                         aaTAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGCAGCCGTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGtt                                                                                                                       <  2:563794/251‑1 (MQ=255)
                                                                                                                                 tACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACa                                                                                                                >  1:90867/1‑250 (MQ=255)
                                                                                                                                     aTTGCAGAGAAACGTTGAGGTTTGCTGTTCTCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTAACGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTTTCTCCGGCCTTTTCGCCCGGATGGATGTCATGTACGGATAAATGGCACATAGTTCTCCTCTCTAGCTGAATTGAAACGATTCAGCCTATATGAGAAAAAAAGAGCCAACCTGGCTTAGGGTTaacgaataaatc                                                                                                           >  2:581567/1‑241 (MQ=255)
                                                                                                                                          cGAGGAACTTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCt                                                                                                      <  2:653293/251‑1 (MQ=255)
                                                                                                                                               aaCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTGCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAAc                                                                                                 <  2:201026/251‑1 (MQ=255)
                                                                                                                                                aCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTATACTTCTGCTACATTACAGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCATGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTCTGAGAAAAAAAGCGCCAACCTGGCTTAGGGATAAAGACAAGATCTCGCTGAAACg                                                                                                >  2:296878/1‑251 (MQ=255)
                                                                                                                                                      tggTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGCTGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGTTTCAGCCTCTATGCGAAAAAAAGAGCCAACCTGGCTTAGGGTTAAAGACAAAATAGCGATGAAACGTTTCa                                                                                           >  2:218892/2‑250 (MQ=255)
                                                                                                                                                      aGGTTTGCTGTTCGCGCTCCAGCATGCCATTGTCGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCATCGACCTGTCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACGTGAACATAGTTCTCCTCTCTTGCTGAATTGAATCGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCaa                                                                                          <  1:840918/251‑1 (MQ=255)
                                                                                                                                                      aGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGttcca                                                                                           <  1:448472/251‑4 (MQ=255)
                                                                                                                                                         tttGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGa                                                                                        >  2:94610/1‑250 (MQ=255)
                                                                                                                                                               gTTCTCGCGCCAGACTGCCACTTACGTTGCCTTCTGCTACTTTACCGGCCATCACCAGCGCCTCAGCGACCTGGCGAATTGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATTGATGTCCTTTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCGCTATGAGAAAAAAAGCGCCAACCTGGCTTAGTGTTAAAGACAAGATCGCGCTGAACCTTTTCAAGAAAGCata                                                                                 <  2:39612/251‑1 (MQ=255)
                                                                                                                                                                  cgcgcgCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGCTTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCataat                                                                               >  1:693226/1‑250 (MQ=255)
                                                                                                                                                                   ccgcgcCAGCATGCCATTGACGTACCCGTCTGCTACATTACCGGCCTGCCCCAGCGCAGCAGCGACCTGGCTAATCTCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGGCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCACTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATAc                                                                             <  2:249150/250‑1 (MQ=255)
                                                                                                                                                                     gcgcCAGCATGCCCTTGCCGTTGCCTTCTGCTACATTCCCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCAGAAACGTTTCAAGAAAGCATAATACtt                                                                           <  2:398969/251‑1 (MQ=255)
                                                                                                                                                                       cccAGCAGGCCATTGACGTAGCCTTCTGCTACCTTAACGGCCTGCACCAGCGCCGCAGCGACCTTGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGCAACGTTTCAAGAAAGCATAATACTTCt                                                                         <  2:693226/250‑1 (MQ=255)
                                                                                                                                                                        ccAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTg                                                                        <  1:85943/251‑1 (MQ=255)
                                                                                                                                                                            cATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCACTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTc                                                                    >  1:488487/1‑251 (MQ=255)
                                                                                                                                                                                  aTTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACg                                                              <  1:371921/251‑1 (MQ=255)
                                                                                                                                                                                  aTTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACg                                                              <  1:115774/251‑1 (MQ=255)
                                                                                                                                                                                        gTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCAGCAGCGACATGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCATGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGTATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGGTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCacgg                                                         >  1:60678/1‑247 (MQ=255)
                                                                                                                                                                                          aGCCTTCTGCTCCATTACCGGCCTGCACCAGCGCCGCAGCTACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGaa                                                       <  2:349702/250‑1 (MQ=255)
                                                                                                                                                                                               tCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTATT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGcc                                                   <  1:123032/249‑1 (MQ=255)
                                                                                                                                                                                                  gCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTaa                                               >  1:198755/1‑250 (MQ=255)
                                                                                                                                                                                                       attcccGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTAt                                         <  1:649971/247‑1 (MQ=255)
                                                                                                                                                                                                             ggCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCa                                    >  2:48038/1‑250 (MQ=255)
                                                                                                                                                                                                                   cacatgcgcCGCAGCGACTTTGCGAATCGCCTCT‑GTTTGTCGTCGGCCTTGTCGCCCGGATGGATGTCATGTACGGATAACTTGCACATAGTTCTCCTCTCTTTCTGAATTGAAACGATTCTGCCTCTATGAGAAAAAAAGCGCCAACCGGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTACAAGAAAGCATAATACTACTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGtta                             <  2:10107/246‑2 (MQ=255)
                                                                                                                                                                                                                   cACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGttt                             <  1:8794/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGGCGAATCGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAtttt           >  2:409145/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aTCGCCTCTT‑TTTGTATCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAt    >  1:741160/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cGCCTCTT‑TTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttt  <  1:48038/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTT  >  NC_000913/2263195‑2263684

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