Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,538,491 C→T G92D (GGC→GAC)  nanM ← N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,4910CT100.0% 144.9 / NA 42G92D (GGC→GAC) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/23);  total (19/23)

TTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538356‑4538625
                                                                                                                                       |                                                                                                                                      
ttGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcgc                                                                                                                                     <  2:647999/139‑1 (MQ=255)
  gtgtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  1:95621/139‑1 (MQ=255)
    gtACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  2:143356/139‑1 (MQ=255)
               acaTGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGAt                                                                                                                      <  1:223699/139‑1 (MQ=255)
                      ccGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATc                                                                                                               >  1:785891/1‑139 (MQ=255)
                         ccATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAAt                                                                                                            <  2:561431/139‑1 (MQ=255)
                              cccATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATg                                                                                                       >  2:774603/1‑139 (MQ=255)
                                     gcgcGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                <  1:774603/139‑1 (MQ=255)
                                                     tttAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGt                                                                                <  2:772061/139‑1 (MQ=255)
                                                         aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCagag                                                                                                  <  2:798724/117‑1 (MQ=255)
                                                         aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCagag                                                                                                  >  1:798724/1‑117 (MQ=255)
                                                                         tGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTaa                                                            >  1:273742/1‑139 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 <  1:366823/52‑1 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCaaa                                                                                                                                 >  2:366823/1‑52 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    <  2:629024/113‑1 (MQ=255)
                                                                                           cGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAAt                                                                    >  1:629024/1‑113 (MQ=255)
                                                                                                 tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATgtcg                                                                                                                                      <  1:361264/41‑1 (MQ=255)
                                                                                                 tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  <  1:20267/61‑1 (MQ=255)
                                                                                                 tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGcc                                                                                                                  >  2:20267/1‑61 (MQ=255)
                                                                                                 tAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGc                                    <  2:273981/139‑1 (MQ=255)
                                                                                                  aaaTACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGcc                                   <  2:273742/139‑1 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGt                                                                                >  1:149837/1‑88 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGt                                                                                <  2:149837/88‑1 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     >  2:93104/1‑131 (MQ=255)
                                                                                                        cTGAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTgg                                     <  1:93104/131‑1 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   <  2:514969/35‑1 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCa                                                                                                                                   >  1:514969/1‑35 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                >  2:118630/1‑70 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                <  1:118630/70‑1 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATc                           >  2:782955/1‑139 (MQ=255)
                                                                                                          gAGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCAGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATc                           >  2:737365/1‑139 (MQ=255)
                                                                                                           aGTCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATcc                          >  1:64847/1‑139 (MQ=255)
                                                                                                             tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCt                                         >  2:231406/1‑122 (MQ=255)
                                                                                                             tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCt                                         <  1:231406/122‑1 (MQ=255)
                                                                                                             tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAg                        >  2:790815/1‑139 (MQ=255)
                                                                                                             tCAAGCCCTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAg                        >  2:260162/1‑139 (MQ=255)
                                                                                                                  cccTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                >  1:425804/1‑62 (MQ=255)
                                                                                                                  cccTCGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGt                                                                                                <  2:425804/62‑1 (MQ=255)
                                                                                                                    cTCGCGGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTAc                 <  2:641787/139‑1 (MQ=255)
                                                                                                                      cGCTGTTTTTGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCa               >  2:274390/1‑139 (MQ=255)
                                                                                                                             tttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTa        <  1:782955/139‑1 (MQ=255)
                                                                                                                              ttGCCAATGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTAc       <  2:779858/139‑1 (MQ=255)
                                                                                                                                   aaTGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTgcg  <  2:64847/139‑1 (MQ=255)
                                                                                                                                   aaTGTCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTgcg  <  1:790815/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TTGTGTACAAAAGTCACATGGCCCGCCATGCCCATCGGCGCGTGCGACATCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCCAGGGAATGCAGCTAACGCTGTCCATTTTTTATCTTTGGCCTGTGTATCCAGCTTGTACCATGCCGTACCTGCG  >  NC_000913/4538356‑4538625

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