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
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cagcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGATCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGt                                                                                                                                                                                                                                                <  1:416872/251‑1 (MQ=255)
     aTACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCg                                                                                                                                                                                                                                           >  1:487356/1‑251 (MQ=255)
          gATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCtt                                                                                                                                                                                                                                      <  1:287102/251‑1 (MQ=255)
              gCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTc                                                                                                                                                                                                                                   <  1:9939/250‑1 (MQ=255)
               ccACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGc                                                                                                                                                                                                                                 >  1:398969/1‑251 (MQ=255)
                     aCGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCggat                                                                                                                                                                                                                           >  1:166485/1‑251 (MQ=255)
                       ccgcccTGACCGTAACCCCAGGTGACGCCGTCCGGGAACTGGAATACCTGAACGCCGGGTTTCAGCACCTGATCGCGGGTTGCGGTAGTGCTGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCAGGCCATTGACGTAGCCGTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCGTTTCGCCCggatgg                                                                                                                                                                                                                         <  2:663083/247‑1 (MQ=255)
                       ccTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTGTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCggatgg                                                                                                                                                                                                                         <  2:3982/250‑1 (MQ=255)
                            cTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTc                                                                                                                                                                                                                    >  1:305136/1‑251 (MQ=255)
                               cccGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCGGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATTCCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTGTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTg                                                                                                                                                                                                                 <  2:349655/250‑1 (MQ=255)
                                 cGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTa                                                                                                                                                                                                               <  1:861164/251‑1 (MQ=255)
                                  gTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTAc                                                                                                                                                                                                              >  1:121800/1‑251 (MQ=255)
                                         cAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAAc                                                                                                                                                                                                       <  1:205536/251‑1 (MQ=255)
                                          aGGTGCCGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCGGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACt                                                                                                                                                                                                      <  2:19955/251‑1 (MQ=255)
                                          aGGTGACGCCTGCCGGGAACTGAAATTCCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGGTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACt                                                                                                                                                                                                      <  2:852868/251‑1 (MQ=255)
                                                 gacgtcagtgcactgtaatcccttgaCGCCGGTTTTCAGCACCGGATCGCGGGTGGCGGTAGTGCCGTGTGGAATAGCAATTCCCGGGCAGAGTAACGTTTAGATTGGTTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTGCTGCTACATTCCCGGCCTGCACCTGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCGCGGATGGATGTCCTGTACGGATTACTGTAACAt                                                                                                                                                                                               <  2:152293/226‑1 (MQ=255)
                                                     tCCGGGAACTGACATACCTGAACGCCGGTTTTAAGCGCCAGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACACTTGCCGAGGAACGTTGAGGTTTACTGTTCGCGCGCCAGCATGCCATTGGCGTAGCCTTCTGCTACATTACCGGCCTGCACCCGCGCCGCAGAGACCCTGCGAATCGCCTCCTTTTGTATCCGGCCTTTTCGCCCGGGTGAAGGTCCTGTACGGATAACTGGAACATAGtt                                                                                                                                                                                           >  2:613854/1‑251 (MQ=255)
                                                      ccGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTtc                                                                                                                                                                                          >  1:349702/1‑251 (MQ=255)
                                                      cagggAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGACGCAGCGACCTGGCGCAGCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGCTGTCCTGTACGGATAACTGGAACATAGttt                                                                                                                                                                                          >  2:285002/3‑250 (MQ=255)
                                                           aattGAAATACGTGAACGCCGGTGTTCAGCAATGGATCGCGGGGGTCGTTAGGGTAGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGACAGCAGGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCCCCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGACTTGTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTctcctc                                                                                                                                                                                     <  2:191/248‑1 (MQ=255)
                                                                aaaTACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTg                                                                                                                                                                                <  1:345392/251‑1 (MQ=255)
                                                                     ccTGAACGCCGGTTTTCAGCACCTGATCTCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGa                                                                                                                                                                            >  1:534791/1‑250 (MQ=255)
                                                                        gcaCGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGACCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGCTGGATGTCCTGTACGGCTAAATGGACCATAGTTCTCCTCTCTTGATGAATTg                                                                                                                                                                        >  2:692692/3‑251 (MQ=255)
                                                                        gAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTg                                                                                                                                                                        >  1:870389/1‑251 (MQ=255)
                                                                        gAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTg                                                                                                                                                                        <  1:47271/251‑1 (MQ=255)
                                                                                ggttgcAGCTCCTGATCGCGGGTGTCGGTAGTGTCGTGTGGAATAGCAATGCCATTGCCGAGGAACGTTGAGGTTTGCTGGTCGCGCGCCAGCATGCCATTGCCGTAGCCTTCTGCTACATTGCCGGCCTGCACCAGCGCCGCATCGTCCTGGCGATTCGCCTCTTTTTGTCTCCGGCCTTGTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTTAAACGAtt                                                                                                                                                                <  2:5930/246‑1 (MQ=255)
                                                                                  tttcagacccTGATCGCGGTGGTCGGGAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTc                                                                                                                                                               <  2:487356/242‑1 (MQ=255)
                                                                                  tttCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTc                                                                                                                                                               <  1:692692/250‑1 (MQ=255)
                                                                                             gATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGATTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCctctctg                                                                                                                                                     >  1:524866/1‑246 (MQ=255)
                                                                                               tCGCGGGTGTAGGTAGTGACGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGGTAACTGGAACATCGTTCTCCTCTCTTGCTGAATTGACAAGATTCAGCCTCTATgtaa                                                                                                                                                  >  2:408998/1‑247 (MQ=255)
                                                                                                cgcgGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGaaa                                                                                                                                                <  1:842856/251‑1 (MQ=255)
                                                                                                   gggTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGaaaaa                                                                                                                                              >  1:485942/1‑250 (MQ=255)
                                                                                                              gTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCt                                                                                                                                  >  1:29489/1‑251 (MQ=255)
                                                                                                               ggCCGTGTGGAATAGCCATACCAGTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTg                                                                                                                                 <  2:236682/250‑1 (MQ=255)
                                                                                                                 ccGTGTGGAATAGCAATACCATTGCCGAGGACCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTCCATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGc                                                                                                                               <  2:485942/251‑1 (MQ=255)
                                                                                                                      tgGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGAAGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAgg                                                                                                                          >  2:323674/1‑251 (MQ=255)
                                                                                                                         aaTAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGCAGCCGTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGtt                                                                                                                       <  2:563794/251‑1 (MQ=255)
                                                                                                                                 tACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACa                                                                                                                >  1:90867/1‑250 (MQ=255)
                                                                                                                                     aTTGCAGAGAAACGTTGAGGTTTGCTGTTCTCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTAACGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTTTCTCCGGCCTTTTCGCCCGGATGGATGTCATGTACGGATAAATGGCACATAGTTCTCCTCTCTAGCTGAATTGAAACGATTCAGCCTATATGAGAAAAAAAGAGCCAACCTGGCTTAGGGTTaacgaataaatc                                                                                                           >  2:581567/1‑241 (MQ=255)
                                                                                                                                          cGAGGAACTTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCt                                                                                                      <  2:653293/251‑1 (MQ=255)
                                                                                                                                               aaCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTGCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAAc                                                                                                 <  2:201026/251‑1 (MQ=255)
                                                                                                                                                aCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTATACTTCTGCTACATTACAGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCATGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTCTGAGAAAAAAAGCGCCAACCTGGCTTAGGGATAAAGACAAGATCTCGCTGAAACg                                                                                                >  2:296878/1‑251 (MQ=255)
                                                                                                                                                      tggTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGCTGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGTTTCAGCCTCTATGCGAAAAAAAGAGCCAACCTGGCTTAGGGTTAAAGACAAAATAGCGATGAAACGTTTCa                                                                                           >  2:218892/2‑250 (MQ=255)
                                                                                                                                                      aGGTTTGCTGTTCGCGCTCCAGCATGCCATTGTCGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCATCGACCTGTCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACGTGAACATAGTTCTCCTCTCTTGCTGAATTGAATCGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCaa                                                                                          <  1:840918/251‑1 (MQ=255)
                                                                                                                                                      aGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGttcca                                                                                           <  1:448472/251‑4 (MQ=255)
                                                                                                                                                         tttGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGa                                                                                        >  2:94610/1‑250 (MQ=255)
                                                                                                                                                               gTTCTCGCGCCAGACTGCCACTTACGTTGCCTTCTGCTACTTTACCGGCCATCACCAGCGCCTCAGCGACCTGGCGAATTGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATTGATGTCCTTTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCGCTATGAGAAAAAAAGCGCCAACCTGGCTTAGTGTTAAAGACAAGATCGCGCTGAACCTTTTCAAGAAAGCata                                                                                 <  2:39612/251‑1 (MQ=255)
                                                                                                                                                                  cgcgcgCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGCTTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCataat                                                                               >  1:693226/1‑250 (MQ=255)
                                                                                                                                                                   ccgcgcCAGCATGCCATTGACGTACCCGTCTGCTACATTACCGGCCTGCCCCAGCGCAGCAGCGACCTGGCTAATCTCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGGCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCACTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATAc                                                                             <  2:249150/250‑1 (MQ=255)
                                                                                                                                                                     gcgcCAGCATGCCCTTGCCGTTGCCTTCTGCTACATTCCCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCAGAAACGTTTCAAGAAAGCATAATACtt                                                                           <  2:398969/251‑1 (MQ=255)
                                                                                                                                                                       cccAGCAGGCCATTGACGTAGCCTTCTGCTACCTTAACGGCCTGCACCAGCGCCGCAGCGACCTTGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGCAACGTTTCAAGAAAGCATAATACTTCt                                                                         <  2:693226/250‑1 (MQ=255)
                                                                                                                                                                        ccAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTg                                                                        <  1:85943/251‑1 (MQ=255)
                                                                                                                                                                            cATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCACTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTc                                                                    >  1:488487/1‑251 (MQ=255)
                                                                                                                                                                                  aTTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACg                                                              <  1:371921/251‑1 (MQ=255)
                                                                                                                                                                                  aTTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACg                                                              <  1:115774/251‑1 (MQ=255)
                                                                                                                                                                                        gTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCAGCAGCGACATGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCATGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGTATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGGTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCacgg                                                         >  1:60678/1‑247 (MQ=255)
                                                                                                                                                                                          aGCCTTCTGCTCCATTACCGGCCTGCACCAGCGCCGCAGCTACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGaa                                                       <  2:349702/250‑1 (MQ=255)
                                                                                                                                                                                               tCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTATTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGcc                                                   <  1:123032/249‑1 (MQ=255)
                                                                                                                                                                                                  gCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTaa                                               >  1:198755/1‑250 (MQ=255)
                                                                                                                                                                                                       attcccGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTAt                                         <  1:649971/247‑1 (MQ=255)
                                                                                                                                                                                                             ggCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCa                                    >  2:48038/1‑250 (MQ=255)
                                                                                                                                                                                                                   cacatgcgcCGCAGCGACTTTGCGAATCGCCTCTGTTTGTCGTCGGCCTTGTCGCCCGGATGGATGTCATGTACGGATAACTTGCACATAGTTCTCCTCTCTTTCTGAATTGAAACGATTCTGCCTCTATGAGAAAAAAAGCGCCAACCGGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTACAAGAAAGCATAATACTACTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGtta                             <  2:10107/246‑2 (MQ=255)
                                                                                                                                                                                                                   cACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGttt                             <  1:8794/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGGCGAATCGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAtttt           >  2:409145/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aTCGCCTCTTTTTGTATCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAt    >  1:741160/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cGCCTCTTTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACATAGTTCTCCTCTCTTGCTGAATTGAAACGATTCAGCCTCTATGAGAAAAAAAGCGCCAACCTGGCTTAGGGTTAAAGACAAGATCGCGCTGAAACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttt  <  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: