Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,684,971 G→A G179G (GGC→GGT glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,9710GA100.0% 136.9 / NA 44G179G (GGC→GGTglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/19);  total (25/19)

GTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACC  >  NC_000913/2684837‑2685099
                                                                                                                                      |                                                                                                                                
gtaACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCgg                                                                                                                              <  1:940221/139‑1 (MQ=255)
 taACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGa                                                                                                                             <  1:861413/139‑1 (MQ=255)
          gAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCaga                                                                                                                    >  2:445194/1‑139 (MQ=255)
          gAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATACGCaga                                                                                                                    >  1:1024592/1‑139 (MQ=255)
                 aGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacca                                                                                                             >  2:482002/1‑139 (MQ=255)
                   gAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAaccacc                                                                                                           <  2:521432/139‑1 (MQ=255)
                       cGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATa                                                                                                       >  1:640888/1‑139 (MQ=255)
                        gggTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATACGCAGAGAAACCACCGATaa                                                                                                      <  2:1024592/139‑1 (MQ=255)
                                   gACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                           <  1:445194/139‑1 (MQ=255)
                                      gCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTtg                                                                                        <  1:78318/139‑1 (MQ=255)
                                          gcagcaACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGtt                                                                                    >  1:428331/1‑139 (MQ=255)
                                              caACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCttt                                                                                >  2:281754/1‑139 (MQ=255)
                                              caACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCttt                                                                                >  1:879851/1‑139 (MQ=255)
                                                     gCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGt                                                                         >  1:518885/1‑139 (MQ=255)
                                                          cAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                                                                    <  1:281754/139‑1 (MQ=255)
                                                                         cAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTc                                                     >  1:161199/1‑139 (MQ=255)
                                                                           aCGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGa                                                   <  2:640888/139‑1 (MQ=255)
                                                                                cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                           >  1:632193/1‑94 (MQ=255)
                                                                                cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGc                                                                                           <  2:632193/94‑1 (MQ=255)
                                                                                cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGa                                              >  1:882577/1‑139 (MQ=255)
                                                                                cAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGa                                              >  1:515836/1‑139 (MQ=255)
                                                                                 aGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGAc                                             >  2:74815/1‑139 (MQ=255)
                                                                                  ggTAAGCACCGATGCTGTCTGCGATTTCACGCATTTTCGCCCAGTCCACAACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGCCTTGTGTTCTTTGGCTTGGTTCACAAGATCGGCCTAGCCGATATGAcc                                            >  2:99836/1‑139 (MQ=255)
                                                                                     aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                         <  2:428331/139‑1 (MQ=255)
                                                                                     aaGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGt                                         <  2:879851/139‑1 (MQ=255)
                                                                                        cACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacca                                                                                                             >  2:965090/1‑68 (MQ=255)
                                                                                        cACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAacca                                                                                                             <  1:965090/68‑1 (MQ=255)
                                                                                          ccGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCAt                                    >  2:254069/1‑139 (MQ=255)
                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAgaga                                                                                                                  <  1:73287/60‑1 (MQ=255)
                                                                                           cGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAgaga                                                                                                                  >  2:73287/1‑60 (MQ=255)
                                                                                            gATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCg                                  <  2:518885/139‑1 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAgaga                                                                                                                  <  1:519378/58‑1 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAgaga                                                                                                                  >  2:519378/1‑58 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGc                                                           <  2:625252/113‑1 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGc                                                           >  1:625252/1‑113 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                          <  2:896985/114‑1 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCg                                                          >  1:896985/1‑114 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                 >  2:189620/1‑139 (MQ=255)
                                                                                             aTGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGa                                 >  1:599397/1‑139 (MQ=255)
                                                                                                 tGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATAcc                             <  1:189620/139‑1 (MQ=255)
                                                                                                     aGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTgg                                                                              <  2:589868/86‑1 (MQ=255)
                                                                                                     aGCGATTTCACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTgg                                                                              >  1:589868/1‑86 (MQ=255)
                                                                                                             cACGCATTTTCGCCCAGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGAt                 >  2:926009/1‑139 (MQ=255)
                                                                                                                            aGTCCACCACACCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTAcc  >  2:13700/1‑139 (MQ=255)
                                                                                                                                   cacaCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCgg                                                            <  2:599002/74‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
GTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGATATGACCGGTAGCATCGATACCGTAAGGAACGATGTTGTACAGTTTACC  >  NC_000913/2684837‑2685099

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