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% 80.0 / NA 22intergenic (‑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 (9/13);  total (9/13)

ACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCT  >  NC_000913/2259177‑2259662
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
accaGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAAtt                                                                                                                                                                                                                                              >  1:317148/1‑250 (MQ=255)
       gTAAATTTTATTTTACCTGGATTTGAAACCGCATAATTTAATGATTCATGTGAATAGCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAACCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCag                                                                                                                                                                                                                                      <  2:107286/251‑1 (MQ=255)
                                aaTCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGc                                                                                                                                                                                                             >  2:336564/1‑251 (MQ=255)
                                                    tCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGAtt                                                                                                                                                                                          >  2:91004/1‑250 (MQ=255)
                                                    tCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGAtt                                                                                                                                                                                          <  1:65538/250‑1 (MQ=255)
                                                                    cGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCagagag                                                                                                                                                                                                                                  <  1:13876/194‑1 (MQ=255)
                                                                    cGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCagagag                                                                                                                                                                                                                                  >  2:13876/1‑194 (MQ=255)
                                                                     gACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACaaa                                                                                                                                                                        <  1:336564/251‑1 (MQ=255)
                                                                           aaTATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGGCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATGACCGCCAGAGATGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAg                                                                                                                                                                  <  2:355646/251‑1 (MQ=255)
                                                                                   ccGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTtgt                                                                                                                                                          <  1:360016/251‑1 (MQ=255)
                                                                                           ataATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGAt                                                                                                                                                  <  1:38086/251‑1 (MQ=255)
                                                                                                           ccTTTTCTCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACaa                                                                                                                                  <  2:186805/251‑1 (MQ=255)
                                                                                                              gttCGCGCATAGCACGTCCTCTGGGAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCAGATTTAAACAGGTGCTCATAGTATGCATGTTTTACGTAAAGGCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATTCGATGAAATGCCTGAATTACAAAAGCGAGGAGGTTGTGACGGGATGCACAGAACAGAACAAt                                                                                                                                 <  2:56909/248‑1 (MQ=255)
                                                                                                                                   tGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAACATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTtgtg                                                                                                                                                         >  2:137538/1‑204 (MQ=255)
                                                                                                                                   tGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTtgtg                                                                                                                                                         <  1:137538/204‑1 (MQ=255)
                                                                                                                                       aaaaaaTTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCt                                                                                                      <  1:91004/251‑1 (MQ=255)
                                                                                                                                                 cAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGGGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATa                                                                                            <  2:51232/251‑1 (MQ=255)
                                                                                                                                                                         taGATATTGCCAAGCCGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAATCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAAAAAATAATCTGCCCTAGAGGAAATAACCTGAGTGTTTGATGAACGtatg                                                                    >  2:116050/1‑248 (MQ=255)
                                                                                                                                                                                      gACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTg                                                        >  1:81798/1‑250 (MQ=255)
                                                                                                                                                                                                      gtTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAttt                                        <  1:91305/250‑1 (MQ=255)
                                                                                                                                                                                                             tttAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACAGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAt                                >  2:252994/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACAGGTCATTctct  >  2:312342/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ACCAGGCGTAAATTTTATTTTACCTGGATTTGAATCCGCATAATTTAATGATTCATGTGAATATCCGGCGACATCAATATTCGCCGAACCTATAATTACGACATAATCCTTTTCGCGCATAGCACGTCCTCTGGTAAAAAATTACCAACCCTGTCAATCCACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCT  >  NC_000913/2259177‑2259662

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