Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,464,629 G→A S27S (AGC→AGT treC ← trehalose‑6‑P hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,464,6290GA100.0% 114.6 / NA 36S27S (AGC→AGTtreCtrehalose‑6‑P hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/17);  total (19/17)

CGCTACGTCGTAACCGTTATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCGCTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGC  >  NC_000913/4464506‑4464745
                                                                                                                           |                                                                                                                    
cGCTACGTCGTAACCGTTATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTcc                                                                                                       >  2:865487/1‑139 (MQ=255)
   tACGTCGTAACCGTTATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGa                                                                                                    >  1:783933/1‑139 (MQ=255)
   tACGTCGTAACCGTTATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGa                                                                                                    >  1:660103/1‑139 (MQ=255)
           aaCCGTTATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCtt                                                                                            <  1:814020/139‑1 (MQ=255)
                ttATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGAt                                                                                       <  1:305155/139‑1 (MQ=255)
                        cTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGa                                                                                        >  2:518127/1‑130 (MQ=255)
                        cTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGa                                                                                        <  1:518127/130‑1 (MQ=255)
                                     tAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGt                                                                  >  1:383610/1‑139 (MQ=255)
                                              gTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGccacc                                                         <  2:816293/139‑1 (MQ=255)
                                                     aGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTgggg                                                  <  2:383610/139‑1 (MQ=255)
                                                     aGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTgggg                                                  <  1:865487/139‑1 (MQ=255)
                                                            aTCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGa                                           >  1:35748/1‑139 (MQ=255)
                                                                                 aTAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGccac                                                          <  2:390261/103‑1 (MQ=255)
                                                                                 aTAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGccac                                                          >  1:390261/1‑103 (MQ=255)
                                                                                  tAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATaa                                                                                     >  1:818292/1‑75 (MQ=255)
                                                                                  tAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATaa                                                                                     <  2:818292/75‑1 (MQ=255)
                                                                                     tCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAAc                                                                       <  2:440346/86‑1 (MQ=255)
                                                                                     tCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAAc                                                                       >  1:440346/1‑86 (MQ=255)
                                                                                           tgtTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCt                          >  1:425732/1‑125 (MQ=255)
                                                                                           tgtTGGATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCt                          <  2:425732/125‑1 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCgtg                                                                                                             <  1:309/38‑1 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCgtg                                                                                                             >  2:309/1‑38 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGccacc                                                         >  1:376743/1‑90 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGccacc                                                         <  2:376743/90‑1 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCAt                                      <  2:694490/109‑1 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCAt                                      >  1:694490/1‑109 (MQ=255)
                                                                                               ggATAACGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGgcgc        >  1:856059/1‑139 (MQ=255)
                                                                                                   aaCGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCAt                                      >  2:865041/1‑105 (MQ=255)
                                                                                                   aaCGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCAt                                      <  1:865041/105‑1 (MQ=255)
                                                                                                    aCGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGaaa                                                                                                  >  2:85308/1‑44 (MQ=255)
                                                                                                    aCGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGaaa                                                                                                  <  1:85308/44‑1 (MQ=255)
                                                                                                     cGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGa                                                                                        <  1:126739/53‑1 (MQ=255)
                                                                                                     cGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGa                                                                                        >  2:126739/1‑53 (MQ=255)
                                                                                                     cGCCACGTAAATCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTgc  >  2:602764/1‑139 (MQ=255)
                                                                                                                tCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGAt                                             <  2:314540/85‑1 (MQ=255)
                                                                                                                tCGCCGGTACCACTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGAt                                             >  1:314540/1‑85 (MQ=255)
                                                                                                                           |                                                                                                                    
CGCTACGTCGTAACCGTTATCGACCTGGGGAGAGACATAAAAGGGGGTTAGCCAGATGGCATCAACGCCCAGTTTATGCAGATAGTCCAGGTGTTGGATAACGCCACGTAAATCGCCGGTACCGCTACCCGTGGTGTCCTGAAAACTCTTTGGATAAATCTGGTAGATAACGCCGTTTTGCCACCAGTGGGGAAGATGAGTCATTACGTTATTCCTGCGAATGCGAGGGGGCGCAATTGC  >  NC_000913/4464506‑4464745

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