Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,944 G→A A122T (GCC→ACC)  iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,9440GA100.0% 180.8 / NA 56A122T (GCC→ACC) iapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/29);  total (27/29)

ATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTA  >  NC_000913/2876808‑2877067
                                                                                                                                        |                                                                                                                           
aTTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACAcc                                                                                                                           >  1:1094152/1‑139 (MQ=255)
 tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACAccc                                                                                                                          >  1:795138/1‑139 (MQ=255)
         tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGc                                                                                                                  <  2:548344/139‑1 (MQ=255)
          cGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                 >  2:909013/1‑139 (MQ=255)
                     taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTATGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                      <  2:424566/139‑1 (MQ=255)
                       ccgcccgcGATAATCGTAAGAGTTGGCATAACGTGACGGGAAGTACGGTGTTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCCTTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                    >  2:259678/1‑139 (MQ=255)
                       ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTg                                                                                                           >  1:868966/1‑132 (MQ=255)
                       ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTg                                                                                                           <  2:868966/132‑1 (MQ=255)
                       ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                    >  2:363018/1‑139 (MQ=255)
                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                  >  2:625535/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                 >  1:1164384/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                 >  2:408634/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                 >  1:528519/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                 >  2:1057324/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                 >  1:885751/1‑139 (MQ=255)
                              cgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                             >  1:250666/1‑139 (MQ=255)
                                  aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                      <  1:166458/126‑1 (MQ=255)
                                  aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                      >  2:166458/1‑126 (MQ=255)
                                          gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                 <  1:1057324/139‑1 (MQ=255)
                                           agCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAg                                                                                                                   >  1:253204/1‑104 (MQ=255)
                                           agCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAg                                                                                                                   <  2:253204/104‑1 (MQ=255)
                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCTCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                              <  2:795138/139‑1 (MQ=255)
                                               ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                            >  2:993085/1‑139 (MQ=255)
                                                gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                             >  2:441573/1‑121 (MQ=255)
                                                gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                             <  1:441573/121‑1 (MQ=255)
                                                gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                           >  2:171159/1‑139 (MQ=255)
                                                 cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                          <  1:243287/139‑1 (MQ=255)
                                                 cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                          <  1:1041142/139‑1 (MQ=255)
                                                    aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                       <  1:909013/139‑1 (MQ=255)
                                                    aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                       <  2:401792/139‑1 (MQ=255)
                                                       gTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                    >  2:462021/1‑139 (MQ=255)
                                                          aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                       >  1:1058836/1‑101 (MQ=255)
                                                          aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                       <  2:1058836/101‑1 (MQ=255)
                                                          aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                                 >  1:248525/1‑139 (MQ=255)
                                                                     ggTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                      <  1:462021/139‑1 (MQ=255)
                                                                       tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                    <  2:885751/139‑1 (MQ=255)
                                                                                 tCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                         <  1:460584/108‑1 (MQ=255)
                                                                                 tCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                         >  2:460584/1‑108 (MQ=255)
                                                                                  cATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                         <  2:1094152/139‑1 (MQ=255)
                                                                                   aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtgt                                        <  1:171159/139‑1 (MQ=255)
                                                                                       aGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                    <  2:1164384/139‑1 (MQ=255)
                                                                                       aGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                    <  1:408634/139‑1 (MQ=255)
                                                                                         gCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                           >  2:391972/1‑98 (MQ=255)
                                                                                         gCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                           <  1:391972/98‑1 (MQ=255)
                                                                                         gCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCt                                  <  2:528519/139‑1 (MQ=255)
                                                                                                         gATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGc                  <  1:625535/139‑1 (MQ=255)
                                                                                                                 ttATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaaa          <  2:250666/139‑1 (MQ=255)
                                                                                                                    tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                    >  1:217338/1‑110 (MQ=255)
                                                                                                                    tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                    <  2:217338/110‑1 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                           >  2:513595/1‑66 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                           <  1:513595/66‑1 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa  <  1:1208926/139‑1 (MQ=255)
                                                                                                                          aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                  <  2:1056783/42‑1 (MQ=255)
                                                                                                                          aTCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                  >  1:1056783/1‑42 (MQ=255)
                                                                                                                            cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                            >  1:950785/1‑78 (MQ=255)
                                                                                                                            cTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                            <  2:950785/78‑1 (MQ=255)
                                                                                                                                        |                                                                                                                           
ATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTA  >  NC_000913/2876808‑2877067

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