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

GCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAC  >  NC_000913/1850668‑1850935
                                                                                                                                          |                                                                                                                                 
gCCTGACCTGCGCCAACATGCGTTAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc                                                                                                                                  <  1:874151/139‑1 (MQ=255)
        tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGATTCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  2:757749/139‑1 (MQ=255)
        tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  2:429213/139‑1 (MQ=255)
         gcgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTc                                                                                                                         <  2:811072/139‑1 (MQ=255)
                              tGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                    >  2:755121/1‑139 (MQ=255)
                                  cTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGat                                                                                                <  2:1002275/139‑1 (MQ=255)
                                     aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  2:149120/1‑110 (MQ=255)
                                     aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  1:149120/110‑1 (MQ=255)
                                     aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  2:637551/1‑110 (MQ=255)
                                     aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  1:637551/110‑1 (MQ=255)
                                     aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                             >  1:155775/1‑139 (MQ=255)
                                        ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:1020960/139‑1 (MQ=255)
                                          ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca                                                                                        >  2:331706/1‑139 (MQ=255)
                                                                           cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGt                                                       >  2:377691/1‑139 (MQ=255)
                                                                                  tCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:4878/97‑1 (MQ=255)
                                                                                  tCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:4878/1‑97 (MQ=255)
                                                                                            gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:218254/1‑87 (MQ=255)
                                                                                            gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:218254/87‑1 (MQ=255)
                                                                                             cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  1:457735/54‑1 (MQ=255)
                                                                                             cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  2:457735/1‑54 (MQ=255)
                                                                                                tttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCAt                                  >  2:92652/1‑139 (MQ=255)
                                                                                                 ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGct                                                                                                        >  2:139943/1‑68 (MQ=255)
                                                                                                 ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGct                                                                                                        <  1:139943/68‑1 (MQ=255)
                                                                                                          tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:427488/73‑1 (MQ=255)
                                                                                                          tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:427488/1‑73 (MQ=255)
                                                                                                            tttCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  2:593143/1‑39 (MQ=255)
                                                                                                            tttCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  1:593143/39‑1 (MQ=255)
                                                                                                              tCCTTATCCGTATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGtt                    <  2:516518/139‑1 (MQ=255)
                                                                                                              tCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGtt                    <  2:581938/139‑1 (MQ=255)
                                                                                                                  tATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGa                >  2:1112910/1‑139 (MQ=255)
                                                                                                                    tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGc              >  1:1064239/1‑139 (MQ=255)
                                                                                                                     ccGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa             <  2:1064239/139‑1 (MQ=255)
                                                                                                                     ccGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa             <  1:755121/139‑1 (MQ=255)
                                                                                                                      cGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAg            <  1:331706/139‑1 (MQ=255)
                                                                                                                       gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGc                                                    <  2:1013407/98‑1 (MQ=255)
                                                                                                                       gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGc                                                    >  1:1013407/1‑98 (MQ=255)
                                                                                                                        cATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAggg          <  1:377691/139‑1 (MQ=255)
                                                                                                                                cAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAc  <  1:966096/139‑1 (MQ=255)
                                                                                                                                cAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAc  <  2:155775/139‑1 (MQ=255)
                                                                                                                                     gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:104155/1‑46 (MQ=255)
                                                                                                                                     gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:104155/46‑1 (MQ=255)
                                                                                                                                     gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggcgg                                        >  1:762195/1‑96 (MQ=255)
                                                                                                                                     gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggcgg                                        <  2:762195/96‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                 
GCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAC  >  NC_000913/1850668‑1850935

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