Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,844,805 T→C L281P (CTC→CCC)  adeD → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,8050TC100.0% 157.6 / NA 43L281P (CTC→CCC) adeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/22);  total (21/22)

AGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCTCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCGACGGCGCG  >  NC_000913/3844679‑3844937
                                                                                                                              |                                                                                                                                    
aGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGa                                                                                                                          >  2:1068161/1‑139 (MQ=255)
     gaagaaGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTa                                                                                                                     <  1:462077/139‑1 (MQ=255)
       agaagGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                                                   <  2:651906/139‑1 (MQ=255)
           ggACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGt                                                                                                               >  2:871602/1‑139 (MQ=255)
            gACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTg                                                                                                              >  1:757896/1‑139 (MQ=255)
               gCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGga                                                                                                           >  1:774258/1‑139 (MQ=255)
                  ggAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGAt                                                                                                        <  2:451051/139‑1 (MQ=255)
                      aTTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCGACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGcc                                                                                                    >  1:625016/1‑139 (MQ=255)
                          cAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCGACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATg                                                                                                <  2:625016/139‑1 (MQ=255)
                                  cATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGacac                                                                                        <  1:621511/139‑1 (MQ=255)
                                       cGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGgag                                                                                                          <  1:255928/116‑1 (MQ=255)
                                       cGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGgag                                                                                                          >  2:255928/1‑116 (MQ=255)
                                               atCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                                                   >  1:26471/1‑99 (MQ=255)
                                               atCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAAc                                                                                                                   <  2:26471/99‑1 (MQ=255)
                                                tCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAAt                                                                          >  1:613700/1‑139 (MQ=255)
                                                    cgAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTcgc                                                                      >  1:728780/1‑139 (MQ=255)
                                                      aaGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTcgccg                                                                    <  2:774258/139‑1 (MQ=255)
                                                      aaGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATGTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGTTGCCTTAATTcgccg                                                                    <  1:787241/139‑1 (MQ=255)
                                                          gTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTg                                                                >  1:1023511/1‑139 (MQ=255)
                                                           tCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGa                                                               <  1:1169776/139‑1 (MQ=255)
                                                             cGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATc                                                             <  1:1068161/139‑1 (MQ=255)
                                                               cTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGa                                                           <  2:390095/139‑1 (MQ=255)
                                                                    cGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGaacaac                                                      >  1:784402/1‑139 (MQ=255)
                                                                            cAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGaa                                                                                              >  2:31417/1‑91 (MQ=255)
                                                                            cAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGaa                                                                                              <  1:31417/91‑1 (MQ=255)
                                                                            cAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCg                                                            >  1:910118/1‑125 (MQ=255)
                                                                            cAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCg                                                            <  2:910118/125‑1 (MQ=255)
                                                                            cAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAAtgtg                                              >  2:751395/1‑139 (MQ=255)
                                                                               cgcgCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCg                                           <  1:961691/139‑1 (MQ=255)
                                                                               cgcgCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCg                                           <  1:220420/139‑1 (MQ=255)
                                                                                 cgcTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCt                                         <  1:751395/139‑1 (MQ=255)
                                                                                    tGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCa                                      >  1:371629/1‑139 (MQ=255)
                                                                                                    cGAATTTAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTc                      >  1:891537/1‑139 (MQ=255)
                                                                                                         ttAACAGTCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAAGACACATCg                                                                                    >  1:772027/1‑72 (MQ=255)
                                                                                                         ttAACAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCg                                                                                    <  2:772027/72‑1 (MQ=255)
                                                                                                           aaCAGCCCGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCt               >  2:319865/1‑139 (MQ=255)
                                                                                                                cccGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCg          >  2:1476/1‑139 (MQ=255)
                                                                                                                  cGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGc                                                                                                     >  2:1136373/1‑46 (MQ=255)
                                                                                                                  cGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGc                                                                                                     <  1:1136373/46‑1 (MQ=255)
                                                                                                                  cGCAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCGAc        >  1:143036/1‑139 (MQ=255)
                                                                                                                    cAATGCATGCCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCGACgg      <  1:196289/139‑1 (MQ=255)
                                                                                                                      atgcatgcCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCGACGgcg    <  2:613700/139‑1 (MQ=255)
                                                                                                                        gcatgcCCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCGACGgcgcg  <  2:379417/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
AGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGCTCTGTACCGATGACCGTAACCCGTGGGAGATCGCCCATGAAGGACACATCGATGCCTTAATTCGCCGCCTGATCGAACAACACAATGTGCCGCTGCATGTGGCATATCGCGTCGCCAGCTGGTCGACGGCGCG  >  NC_000913/3844679‑3844937

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: