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% 189.2 / NA 46intergenic (+113/‑54)sppA/ansAprotease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/25);  total (21/25)

GACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCG  >  NC_000913/1850672‑1850937
                                                                                                                                      |                                                                                                                                   
gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccct                                                                                                                                 <  2:88230/139‑2 (MQ=255)
 aCCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccctg                                                                                                                                >  1:1440251/1‑137 (MQ=255)
            cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGGTCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGtata                                                                                                                    <  2:602800/139‑1 (MQ=255)
                cGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                   >  2:47544/1‑121 (MQ=255)
                cGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                   <  1:47544/121‑1 (MQ=255)
                    aaGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTg                                                                                                            <  2:1440251/139‑1 (MQ=255)
                     aGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGc                                                                                                           >  1:1267276/1‑139 (MQ=255)
                          tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                      <  1:929161/139‑1 (MQ=255)
                          tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                      >  2:1023711/1‑139 (MQ=255)
                          tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                      <  2:1267276/139‑1 (MQ=255)
                             aCTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGa                                                                                                   <  1:597950/139‑1 (MQ=255)
                                 aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                               >  2:209083/1‑139 (MQ=255)
                                    ggCCGGCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCATCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTCTCAGGTGCGCTCCCCTGCGTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  1:1078155/139‑1 (MQ=255)
                                    ggCCGACAGATCGTTGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  1:209083/139‑1 (MQ=255)
                                    ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  1:1023711/139‑1 (MQ=255)
                                        gACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            <  2:747129/103‑1 (MQ=255)
                                        gACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            >  1:747129/1‑103 (MQ=255)
                                                  cgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  2:117231/1‑125 (MQ=255)
                                                  cgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  1:117231/125‑1 (MQ=255)
                                                                   aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  2:1024014/1‑108 (MQ=255)
                                                                   aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  1:1024014/108‑1 (MQ=255)
                                                                       cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTTCCCTGCCt                                                                                                                            <  2:1192926/72‑1 (MQ=255)
                                                                       cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGt                                                         >  1:1433258/1‑139 (MQ=255)
                                                                       cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTACCCTGCCt                                                                                                                            >  1:1192926/1‑72 (MQ=255)
                                                                                       cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccctgc                                                                                                                               >  2:636445/1‑51 (MQ=255)
                                                                                       cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccctgc                                                                                                                               <  1:636445/54‑4 (MQ=255)
                                                                                       cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggcggg                                         <  2:1433258/139‑1 (MQ=255)
                                                                                             ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGAc                                       <  2:896255/135‑1 (MQ=255)
                                                                                             ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGAc                                       >  1:896255/1‑135 (MQ=255)
                                                                                                 ccccATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCggg                                >  2:89179/1‑138 (MQ=255)
                                                                                                 ccccATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCggg                                <  1:89179/138‑1 (MQ=255)
                                                                                                           ccTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGAc                                       >  1:714501/1‑121 (MQ=255)
                                                                                                           ccTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGAc                                       <  2:714501/121‑1 (MQ=255)
                                                                                                            cTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  1:531004/1‑67 (MQ=255)
                                                                                                            cTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  2:531004/67‑1 (MQ=255)
                                                                                                               aTCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg                 <  2:416628/139‑1 (MQ=255)
                                                                                                                tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg                 >  2:1625896/1‑138 (MQ=255)
                                                                                                                tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg                 <  1:1625896/138‑1 (MQ=255)
                                                                                                                   gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                              >  1:1534885/1‑58 (MQ=255)
                                                                                                                   gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                              <  2:1534885/58‑1 (MQ=255)
                                                                                                                   gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa               <  2:1307672/137‑1 (MQ=255)
                                                                                                                   gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAGCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa               >  1:1307672/1‑137 (MQ=255)
                                                                                                                          aTCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg                 >  2:334675/1‑128 (MQ=255)
                                                                                                                          aTCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg                 <  1:334675/128‑1 (MQ=255)
                                                                                                                           tCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  2:889869/52‑1 (MQ=255)
                                                                                                                           tCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  1:889869/1‑52 (MQ=255)
                                                                                                                           tCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTtatata     >  1:272459/1‑139 (MQ=255)
                                                                                                                              ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAg                                                                       <  2:1773423/70‑1 (MQ=255)
                                                                                                                              ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAg                                                                       >  1:1773423/1‑70 (MQ=255)
                                                                                                                              ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                               >  1:824887/1‑78 (MQ=255)
                                                                                                                              ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                               <  2:824887/78‑1 (MQ=255)
                                                                                                                              ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCg  >  2:160775/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCG  >  NC_000913/1850672‑1850937

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