Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,400 +GAAA coding (483/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3981.A100.0% 209.5 / NA 59E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (25/34);  total (25/34)
*NC_0009134,008,3982.A100.0% 211.7 / NA 59E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (25/34);  total (25/34)
*NC_0009134,008,3983.G100.0% 181.3 / NA 59E161? (GAA→NAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (25/34);  total (25/34)
*NC_0009134,008,3984.A100.0% 209.4 / NA 59E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (25/34);  total (25/34)

GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458
                                                                                             ||||                                                                 
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGCAGCGGTGGCG‑‑‑‑a                                                                  >  2:2584046/1‑93 (MQ=255)
gCGAGCGCGTCAACGTCAGTTCCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                  <  1:3422101/93‑1 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 <  1:3182967/93‑1 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                 >  2:1849334/1‑93 (MQ=255)
  gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  2:2756240/1‑92 (MQ=255)
  gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                <  1:3581944/93‑2 (MQ=255)
  gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                <  1:2674900/93‑2 (MQ=255)
  gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                <  1:575495/93‑2 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:2329246/1‑91 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:2283480/1‑91 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  1:1052050/1‑91 (MQ=255)
     cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaaa                                                             <  2:2369141/93‑5 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:2592961/93‑1 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:4432993/93‑1 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 >  2:3254248/1‑93 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:3793677/93‑3 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:3632859/93‑3 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  1:1240529/1‑90 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:1051163/93‑5 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             >  1:3702178/1‑89 (MQ=255)
                                               aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  2:1077404/1‑93 (MQ=255)
                                               aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:760103/93‑1 (MQ=255)
                                               aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:229684/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:3374702/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:1894017/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:1566229/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:401856/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCCGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:3408020/93‑1 (MQ=255)
                                                  gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:3117785/93‑1 (MQ=255)
                                                  gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:1467147/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:4425676/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:1119358/93‑1 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                >  1:4457103/1‑93 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:3526130/93‑1 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               >  1:4576177/1‑93 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:4322943/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  2:4523682/1‑93 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:2669860/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:2958163/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  2:516542/1‑93 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgggg             >  1:2074904/1‑90 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:2058391/1‑93 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  2:4150932/1‑93 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  2:1382101/93‑1 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  1:2192640/1‑93 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:1265185/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:1671891/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:1470272/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCAGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:4571105/1‑93 (MQ=255)
                                                               canaaGNNCATGCGCTGGAAGNNGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:4559497/90‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:2234580/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:3847674/91‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:1013108/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:2030356/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:2669759/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      <  2:4362651/92‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:4466518/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:2475525/2‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:1050044/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:470735/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:4590471/1‑93 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  2:1078745/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:4450069/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:2907437/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:3864513/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:4631955/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:1507239/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4654945/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:3115899/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:2559408/93‑1 (MQ=255)
                                                                    ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:3964472/93‑1 (MQ=255)
                                                                                             ||||                                                                 
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458

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