Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 788,084 C→T A252T (GCC→ACC)  galM ← galactose‑1‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913788,0840CT100.0% 127.4 / NA 37A252T (GCC→ACC) galMgalactose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/21);  total (16/21)

GTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTC  >  NC_000913/787948‑788202
                                                                                                                                        |                                                                                                                      
gttcCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTg                                                                                                                     <  2:739937/140‑1 (MQ=255)
 ttcCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGt                                                                                                                    <  2:885827/140‑1 (MQ=255)
    cccGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAAc                                                                                                                 >  2:356961/1‑140 (MQ=255)
     ccGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACa                                                                                                                >  2:464318/1‑140 (MQ=255)
       gcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAg                                                                                                              <  1:223869/140‑1 (MQ=255)
       gcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAg                                                                                                              <  2:91431/140‑1 (MQ=255)
        cgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGa                                                                                                             <  2:782074/140‑1 (MQ=255)
             gtgtGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCg                                                                                                        <  1:681446/140‑1 (MQ=255)
              tgtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGt                                                                                                       >  1:774012/1‑140 (MQ=255)
              tgtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGt                                                                                                       <  2:103971/140‑1 (MQ=255)
              tgtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGt                                                                                                       >  1:967446/1‑140 (MQ=255)
                   cgccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCg                                                                                                  <  1:482963/140‑1 (MQ=255)
                    gccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGt                                                                                                 >  2:40700/1‑140 (MQ=255)
                    gccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGt                                                                                                 <  1:356961/140‑1 (MQ=255)
                     ccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTa                                                                                                <  1:337225/140‑1 (MQ=255)
                        aGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCg                                                                                                        <  1:97168/129‑1 (MQ=255)
                        aGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCg                                                                                                        >  2:97168/1‑129 (MQ=255)
                        aGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAAcc                                                                                             >  2:1148504/1‑140 (MQ=255)
                        aGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAAcc                                                                                             >  1:1152415/1‑140 (MQ=255)
                               tGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACt                                                                                      <  2:774012/140‑1 (MQ=255)
                               tGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACt                                                                                      <  2:967446/140‑1 (MQ=255)
                                gCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACtt                                                                                     <  1:464318/140‑1 (MQ=255)
                                     aGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACt                                                                                      <  1:259036/134‑1 (MQ=255)
                                     aGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACt                                                                                      >  2:259036/1‑134 (MQ=255)
                                                   gCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGc                                                                  <  1:1148504/140‑1 (MQ=255)
                                                    ccGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCa                                                                 <  2:1171036/140‑1 (MQ=255)
                                                                    gACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCga                                                 <  2:1152415/140‑1 (MQ=255)
                                                                         tCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGAtt                                            >  2:407/1‑140 (MQ=255)
                                                                         tCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGAtt                                            >  1:484761/1‑140 (MQ=255)
                                                                                     tttCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAGACTCACTGGCGATGATTTTGGCGCTGCgg                                <  2:770058/140‑1 (MQ=255)
                                                                                                        aTGCGCCGCCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCgg             >  1:1220815/1‑140 (MQ=255)
                                                                                                           cgccgcCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGa          <  2:1023496/140‑1 (MQ=255)
                                                                                                           cgccgcCACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGa          >  1:1023496/1‑140 (MQ=255)
                                                                                                                ccACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCgg             >  1:95781/1‑132 (MQ=255)
                                                                                                                ccACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCgg             <  2:95781/132‑1 (MQ=255)
                                                                                                                ccACTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGAt     >  2:1007311/1‑140 (MQ=255)
                                                                                                                   cTTTCTTGCCATCGCCTTTGGTCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTc  >  1:341258/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                      
GTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTC  >  NC_000913/787948‑788202

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