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

TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGGTG  >  NC_000913/1850676‑1850940
                                                                                                                                  |                                                                                                                                      
tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               <  1:411691/139‑1 (MQ=255)
tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               <  2:215514/139‑1 (MQ=255)
       aCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      <  1:105385/126‑1 (MQ=255)
       aCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      >  2:105385/1‑126 (MQ=255)
                aaGTCTTGTACTGAGTGGCCTACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTg                                                                                                               <  2:451683/139‑1 (MQ=255)
                 aGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGc                                                                                                              >  2:708656/1‑139 (MQ=255)
                       gTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTc                                                                                                        >  1:257606/1‑139 (MQ=255)
                        tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCg                                                                                                       <  2:399143/139‑1 (MQ=255)
                           tGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGGGCTCCCCTGCCTCACGTATGTACTTTTGCTCTTTCGata                                                                                                    <  1:397790/139‑1 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                                  >  1:593052/1‑139 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  2:57759/139‑1 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  1:708656/139‑1 (MQ=255)
                                  ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca                                                                                             >  2:694165/1‑139 (MQ=255)
                                     aCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                          <  2:257606/139‑1 (MQ=255)
                                       aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatca                                                                                        <  2:713633/139‑1 (MQ=255)
                                       aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatca                                                                                        >  1:554104/1‑139 (MQ=255)
                                                  cacaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCaa                                                                             <  1:639609/139‑1 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGt                                                            >  1:261080/1‑139 (MQ=255)
                                                                               ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaata                                                                                     <  2:215879/102‑1 (MQ=255)
                                                                               ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaata                                                                                     >  1:215879/1‑102 (MQ=255)
                                                                               ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTATGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcggg                                                >  1:473004/1‑139 (MQ=255)
                                                                                         ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATc                                      <  1:694165/139‑1 (MQ=255)
                                                                                                           aTCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg                    <  2:593052/139‑1 (MQ=255)
                                                                                                                 aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATg                                >  2:379174/1‑121 (MQ=255)
                                                                                                                 aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATg                                <  1:379174/121‑1 (MQ=255)
                                                                                                                     tATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         <  1:591807/44‑1 (MQ=255)
                                                                                                                     tATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         >  2:591807/1‑44 (MQ=255)
                                                                                                                        cAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAc       <  2:366306/139‑1 (MQ=255)
                                                                                                                          ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGa                                                                                                      >  2:435385/1‑42 (MQ=255)
                                                                                                                          ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGa                                                                                                      <  1:435385/42‑1 (MQ=255)
                                                                                                                          ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCg     >  1:20198/1‑139 (MQ=255)
                                                                                                                             gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGgtg  >  2:401198/1‑139 (MQ=255)
                                                                                                                             gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGgtg  >  1:526816/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGGTG  >  NC_000913/1850676‑1850940

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: