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% 163.8 / NA 39intergenic (+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/18);  total (21/18)

TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGG  >  NC_000913/1850676‑1850927
                                                                                                                                  |                                                                                                                         
tGCTCCAACATGCCTTAAGGTTTGTACTGAGTAGCCGACAGATCGTGGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCTCTGTTTCGCCCCATATTTCCTTATCCGCGTAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                  <  2:323084/139‑1 (MQ=255)
        cATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGtata                                                                                                          >  2:640877/1‑139 (MQ=255)
                            gAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGatat                                                                                      >  1:998255/1‑139 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                  <  2:272922/110‑1 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                  >  1:272922/1‑110 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                     >  1:1117632/1‑139 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  <  2:998255/139‑1 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  <  2:1117632/139‑1 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  <  1:640877/139‑1 (MQ=255)
                                                              tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                            <  2:982270/99‑1 (MQ=255)
                                                              tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                            >  1:982270/1‑99 (MQ=255)
                                                               aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccct                                                                                                                       <  2:882942/72‑2 (MQ=255)
                                                               aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccct                                                                                                                       >  1:882942/1‑71 (MQ=255)
                                                               aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                         >  2:350269/1‑70 (MQ=255)
                                                               aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                         <  1:350269/70‑1 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                  <  1:1059621/72‑1 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                  >  2:1059621/1‑72 (MQ=255)
                                                                         aTCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  >  1:1127182/1‑98 (MQ=255)
                                                                         aTCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  <  2:1127182/98‑1 (MQ=255)
                                                                              ctTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgg                                    >  2:428180/1‑139 (MQ=255)
                                                                                   cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                    <  2:46663/86‑1 (MQ=255)
                                                                                   cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                    >  1:46663/1‑86 (MQ=255)
                                                                                       gTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTTCCTCACGTATATACTTATGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACAGGCGGGACCa                           >  2:341486/1‑139 (MQ=255)
                                                                                          tCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCg                        >  2:411007/1‑139 (MQ=255)
                                                                                                  tatTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                       >  1:169860/1‑68 (MQ=255)
                                                                                                  tatTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                       <  2:169860/68‑1 (MQ=255)
                                                                                                           aTCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGtt            <  1:429238/134‑1 (MQ=255)
                                                                                                           aTCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGtt            >  2:429238/1‑134 (MQ=255)
                                                                                                           aTCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg       <  1:24602/139‑1 (MQ=255)
                                                                                                            tCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGtata                                                                                                          >  1:516641/1‑39 (MQ=255)
                                                                                                            tCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGtata                                                                                                          <  2:516641/39‑1 (MQ=255)
                                                                                                             ccGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa     >  2:658593/1‑139 (MQ=255)
                                                                                                                cATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAggg  >  1:707199/1‑139 (MQ=255)
                                                                                                                     tATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg       >  2:365520/1‑129 (MQ=255)
                                                                                                                     tATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg       <  1:365520/129‑1 (MQ=255)
                                                                                                                          ggTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                >  1:853718/1‑83 (MQ=255)
                                                                                                                          ggTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                <  2:853718/83‑1 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATTTACTTTTGCTCTTTCGATAtcatt                                                                                  <  2:878267/46‑1 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  <  1:702322/46‑1 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  >  2:702322/1‑46 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                  >  1:878267/1‑46 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                     <  2:631230/75‑1 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                     >  1:631230/1‑75 (MQ=255)
                                                                                                                                  |                                                                                                                         
TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGG  >  NC_000913/1850676‑1850927

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: ‑