Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,850,810 (C)5→4 intergenic (+117/‑50) sppA → / → ansA protease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,850,8060C.100.0% 208.1 / NA 51intergenic (+113/‑54)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (25/26);  total (25/26)

TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATA  >  NC_000913/1850676‑1850934
                                                                                                                                  |                                                                                                                                
tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  2:997929/139‑1 (MQ=255)
  cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  2:624262/137‑1 (MQ=255)
  cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  1:624262/1‑137 (MQ=255)
        cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGtata                                                                                                                 <  2:1161222/139‑1 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTc                                                                                                                        <  2:1048409/111‑1 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTc                                                                                                                        >  1:1048409/1‑111 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                            >  1:214583/1‑139 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                            >  2:58982/1‑139 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  1:104931/139‑1 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  2:214583/139‑1 (MQ=255)
                                  ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  1:1555004/105‑1 (MQ=255)
                                  ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  2:1555004/1‑105 (MQ=255)
                                           cgtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  2:1575851/1‑96 (MQ=255)
                                           cgtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  2:1435111/128‑1 (MQ=255)
                                           cgtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  1:1435111/1‑128 (MQ=255)
                                           cgtcgGCCACATTAGTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  1:1575851/96‑1 (MQ=255)
                                               ggCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATg                                                                          <  1:1901453/139‑1 (MQ=255)
                                                   acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  1:1031102/120‑1 (MQ=255)
                                                   acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  2:1031102/1‑120 (MQ=255)
                                                             ttACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                            >  2:1962817/1‑139 (MQ=255)
                                                              tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  1:219346/77‑1 (MQ=255)
                                                              tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  2:219346/1‑77 (MQ=255)
                                                              tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  1:789509/109‑1 (MQ=255)
                                                              tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  2:789509/1‑109 (MQ=255)
                                                               aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  1:1334686/1‑76 (MQ=255)
                                                               aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  2:1334686/76‑1 (MQ=255)
                                                                 gTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  1:2040459/106‑1 (MQ=255)
                                                                 gTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  2:2040459/1‑106 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                <  1:971028/66‑1 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                >  2:971028/1‑66 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGat                                                                                               >  2:425616/1‑98 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGat                                                                                               <  1:425616/98‑1 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  1:1754468/1‑104 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  2:1754468/104‑1 (MQ=255)
                                                                           ccTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTaca                                              <  2:663529/139‑1 (MQ=255)
                                                                                 cccGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  1:1146503/1‑58 (MQ=255)
                                                                                 cccGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  2:1146503/58‑1 (MQ=255)
                                                                                   cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                           >  2:90542/1‑86 (MQ=255)
                                                                                   cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                           <  1:90542/86‑1 (MQ=255)
                                                                                         ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         >  2:1947497/1‑50 (MQ=255)
                                                                                         ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                         <  1:1947497/50‑1 (MQ=255)
                                                                                                       ccTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTc                  >  2:2042952/1‑139 (MQ=255)
                                                                                                               gCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  1:1265096/1‑60 (MQ=255)
                                                                                                               gCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  2:1265096/60‑1 (MQ=255)
                                                                                                               gCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTaca                                              <  1:111038/103‑1 (MQ=255)
                                                                                                               gCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTaca                                              >  2:111038/1‑103 (MQ=255)
                                                                                                                 aTAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                       <  2:547999/92‑1 (MQ=255)
                                                                                                                 aTAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                       >  1:547999/1‑92 (MQ=255)
                                                                                                                     tATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa            >  2:1972801/1‑131 (MQ=255)
                                                                                                                     tATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa            <  1:1972801/131‑1 (MQ=255)
                                                                                                                       tCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTtatata  >  1:849091/1‑139 (MQ=255)
                                                                                                                          ggTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         <  2:513725/49‑1 (MQ=255)
                                                                                                                          ggTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                         >  1:513725/1‑49 (MQ=255)
                                                                                                                                  |                                                                                                                                
TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATA  >  NC_000913/1850676‑1850934

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