Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,259,422 A→C intergenic (‑126/+297) psuK ← / ← fruA putative pseudouridine kinase/fructose‑specific PTS multiphosphoryl transfer protein FruA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,4220AC100.0% 113.3 / NA 31intergenic (‑126/+297)psuK/fruAputative pseudouridine kinase/fructose‑specific PTS multiphosphoryl transfer protein FruA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/16);  total (15/16)

GCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTC  >  NC_000913/2259182‑2259659
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
gCGTAAATTTTATTTTACCTGGCTTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGcc                                                                                                                                                                                                                                     >  2:148428/1‑251 (MQ=255)
 cGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGcc                                                                                                                                                                                                                                     >  1:366388/1‑250 (MQ=255)
              ttaCCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAAt                                                                                                                                                                                                                        <  2:66272/250‑1 (MQ=255)
                aCCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATggg                                                                                                                                                                                                                     >  1:210452/1‑251 (MQ=255)
                                gCATAATTGAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATCACCAATCCTGTCAATCCACAAGAATATAGATTCTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGATAGGCAAAATGGGCAGACTTGAGCTGGCt                                                                                                                                                                                                     <  2:330694/251‑1 (MQ=255)
                                                           ccGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAAt                                                                                                                                                                           <  2:366388/250‑1 (MQ=255)
                                                            cGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAAtt                                                                                                                                                                          >  2:412113/1‑250 (MQ=255)
                                                             ggCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACCTGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCTCAAAATGCGATTAAATGCCTGAATTAc                                                                                                                                                                        <  2:254303/251‑1 (MQ=255)
                                                                          ttCGCCGAACCTATAATTAAGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAACATGCGATTAAATGCCTGACTTACAAAAGCGAGGAAt                                                                                                                                                           >  2:298279/1‑251 (MQ=255)
                                                                            cGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCAGGTAAAAAATTACCAACCCTGTAAATCCACAAGAATATAGATACTGCAAAGAAGGTTTAAACATGTGTTAATATTTAAACAAGTGCTCATAGTATGCATGTTTTAAGTAAAGTCAATTACCGAAAGAGAGGAAAAATGGGCAGATTTGAGCAGGCACACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAAttg                                                                                                                                                         >  1:29498/1‑250 (MQ=255)
                                                                                                      ccTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATGACAAAAGCGAGGAATTTGTGACGGGATGCACGGAACAGAACaa                                                                                                                               <  2:406731/251‑1 (MQ=255)
                                                                                                         tttCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAAt                                                                                                                              <  2:331680/249‑1 (MQ=255)
                                                                                                          ttCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCACTCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATaa                                                                                                                            >  1:167487/1‑250 (MQ=255)
                                                                                                                            tctGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATc                                                                                                         >  2:201041/1‑251 (MQ=255)
                                                                                                                                                      aaTCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAgtgt                                                                               <  1:372166/251‑1 (MQ=255)
                                                                                                                                                                 atatAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGt                                                                    >  2:172839/1‑251 (MQ=255)
                                                                                                                                                                     aGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGa                                                                <  1:201041/251‑1 (MQ=255)
                                                                                                                                                                                aGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGttt                                                                  <  1:328520/238‑1 (MQ=255)
                                                                                                                                                                                aGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGttt                                                                  >  2:328520/1‑238 (MQ=255)
                                                                                                                                                                                 gACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGa                                                    <  2:64690/251‑1 (MQ=255)
                                                                                                                                                                                     gTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGaa                                                   >  1:158715/1‑248 (MQ=255)
                                                                                                                                                                                        tAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATa                                              <  2:26700/250‑1 (MQ=255)
                                                                                                                                                                                         aaaNATGTTTGCATACTTAAACCTGTGCACATAATATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATTCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTTCCCTTGAGGAAATAACCTTAGCGGTTGATGAACGTGTGAAGTTTTTACGGAAATATaaa                                            <  2:374248/251‑1 (MQ=255)
                                                                                                                                                                                          aaCATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaaa                                           <  1:172839/251‑1 (MQ=255)
                                                                                                                                                                                                     atatTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAAtt                                >  2:154897/1‑251 (MQ=255)
                                                                                                                                                                                                         ttAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATTGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAt                             <  1:148428/250‑1 (MQ=255)
                                                                                                                                                                                                                 gtgCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAAGCTGCCCTTGATGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAGAAAATTACa                    <  2:62032/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTcac               >  1:258009/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAAc           >  1:278232/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ttttACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtc  >  2:352855/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cggaaagtcaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAt                                       >  2:131594/4‑209 (MQ=255)
                                                                                                                                                                                                                                        cGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAt                                       <  1:131594/209‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
GCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTC  >  NC_000913/2259182‑2259659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: