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% 133.8 / 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
                                                                                                                                  |                                                                                                                                      
tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCt                                                                                                                               <  1:324634/139‑1 (MQ=255)
   gcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      >  2:491747/1‑130 (MQ=255)
   gcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      <  1:491747/130‑1 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      <  2:185676/104‑1 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      >  1:185676/1‑104 (MQ=255)
                             aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                                  >  1:817109/1‑139 (MQ=255)
                                ggCCGACATATCGTCGGCCACATTTTTTTTTACGTCCACGAATCCTCTTCCCGCTGGTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  1:927955/139‑1 (MQ=255)
                                ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  1:674646/139‑1 (MQ=255)
                                     aCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               >  2:594255/1‑134 (MQ=255)
                                     aCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  1:594255/134‑1 (MQ=255)
                                                                 gTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCCCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTAc                                                              <  2:665502/139‑1 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCAc                                                                                                                            >  1:401174/1‑75 (MQ=255)
                                                                   cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCAc                                                                                                                            <  2:401174/75‑1 (MQ=255)
                                                                               ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                             <  1:249949/126‑1 (MQ=255)
                                                                               ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                             >  2:249949/1‑126 (MQ=255)
                                                                                    gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    <  1:587886/82‑1 (MQ=255)
                                                                                    gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggcggg                                            <  1:433089/138‑1 (MQ=255)
                                                                                    gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggcggg                                            >  2:433089/1‑138 (MQ=255)
                                                                                    gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTAGata                                                                                                    >  2:587886/1‑82 (MQ=255)
                                                                                     cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc                                                                                                                                     <  2:858949/49‑1 (MQ=255)
                                                                                     cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccc                                                                                                                                     >  1:858949/1‑49 (MQ=255)
                                                                                       gTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                                  <  1:49357/113‑1 (MQ=255)
                                                                                       gTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAt                                                                  >  2:49357/1‑113 (MQ=255)
                                                                                             ccccATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCata                                                                                           <  2:196972/82‑1 (MQ=255)
                                                                                             ccccATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCata                                                                                           >  1:196972/1‑82 (MQ=255)
                                                                                                cATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca                                                                                             >  2:374146/1‑77 (MQ=255)
                                                                                                cATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca                                                                                             <  1:374146/77‑1 (MQ=255)
                                                                                                 atatTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCa                              <  2:817109/139‑1 (MQ=255)
                                                                                                            tCCGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGc                   >  1:530260/1‑139 (MQ=255)
                                                                                                              cGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    <  1:693995/56‑1 (MQ=255)
                                                                                                              cGCATAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    >  2:693995/1‑56 (MQ=255)
                                                                                                                  tAGTATCAGGTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGtt             <  1:761320/139‑1 (MQ=255)
                                                                                                                          ggTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    <  2:305873/44‑1 (MQ=255)
                                                                                                                          ggTGCGCT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    >  1:305873/1‑44 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                          <  2:325889/51‑1 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                          >  1:325889/1‑51 (MQ=255)
                                                                                                                             gcgcT‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGgtg  >  2:159775/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGGTG  >  NC_000913/1850676‑1850940

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