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

CTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAC  >  NC_000913/1850670‑1850935
                                                                                                                                        |                                                                                                                                 
cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                 <  2:139939/139‑1 (MQ=255)
         gcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCAc                                                                                                                       >  1:705280/1‑139 (MQ=255)
             aCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGtat                                                                                                                   >  1:175478/1‑139 (MQ=255)
                                 tGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                               <  2:1078877/139‑1 (MQ=255)
                                        ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGtata                                                                                                                  >  1:1030876/1‑113 (MQ=255)
                                        ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGtata                                                                                                                  <  2:1030876/113‑1 (MQ=255)
                                           aCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                 <  1:935643/96‑1 (MQ=255)
                                           aCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                 >  2:935643/1‑96 (MQ=255)
                                             aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatca                                                                                   >  2:708051/1‑139 (MQ=255)
                                                 cgtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCa                                                                                                                        <  2:516954/98‑1 (MQ=255)
                                                 cgtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCa                                                                                                                        >  1:516954/1‑98 (MQ=255)
                                                      gCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  2:350802/1‑91 (MQ=255)
                                                      gCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  1:350802/91‑1 (MQ=255)
                                                         acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  1:344990/88‑1 (MQ=255)
                                                         acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  2:344990/1‑88 (MQ=255)
                                                           attattTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          >  1:27401/1‑86 (MQ=255)
                                                           attattTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                          <  2:27401/86‑1 (MQ=255)
                                                                    tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAtt                                                            >  1:182104/1‑139 (MQ=255)
                                                                     aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                            >  2:1330237/1‑106 (MQ=255)
                                                                     aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                            <  1:1330237/106‑1 (MQ=255)
                                                                         cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:1213498/1‑104 (MQ=255)
                                                                         cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:1213498/104‑1 (MQ=255)
                                                                                        ccGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggcgg                                        <  2:182104/139‑1 (MQ=255)
                                                                                           cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:776796/1‑86 (MQ=255)
                                                                                           cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:776796/86‑1 (MQ=255)
                                                                                            tGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  1:573306/1‑85 (MQ=255)
                                                                                            tGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  2:573306/85‑1 (MQ=255)
                                                                                            tGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCa                                                                         >  2:1259337/1‑102 (MQ=255)
                                                                                            tGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCa                                                                         <  1:1259337/102‑1 (MQ=255)
                                                                                              tttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc                                                                                                                                  <  1:207784/44‑1 (MQ=255)
                                                                                              tttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc                                                                                                                                  >  2:207784/1‑44 (MQ=255)
                                                                                              tttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCg                                                                                                  >  2:370199/1‑75 (MQ=255)
                                                                                              tttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCg                                                                                                  <  1:370199/75‑1 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  1:14941/82‑1 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  2:14941/1‑82 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                            <  2:131893/80‑1 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                            >  1:131893/1‑80 (MQ=255)
                                                                                                                tATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          >  1:1036002/1‑65 (MQ=255)
                                                                                                                tATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                          <  2:1036002/65‑1 (MQ=255)
                                                                                                                 aTCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAg               <  2:175478/139‑1 (MQ=255)
                                                                                                                  tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                    <  2:518642/53‑1 (MQ=255)
                                                                                                                  tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                    >  1:518642/1‑53 (MQ=255)
                                                                                                                  tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                        <  2:388204/97‑1 (MQ=255)
                                                                                                                  tCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                        >  1:388204/1‑97 (MQ=255)
                                                                                                                    cGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAg            >  2:486236/1‑139 (MQ=255)
                                                                                                                     gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                            >  2:993233/1‑58 (MQ=255)
                                                                                                                     gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtca                                                                                            <  1:993233/58‑1 (MQ=255)
                                                                                                                      cATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                    <  2:218390/49‑1 (MQ=255)
                                                                                                                      cATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                    >  1:218390/1‑49 (MQ=255)
                                                                                                                      cATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAggg          >  1:311127/1‑139 (MQ=255)
                                                                                                                           tATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTtata     >  2:481725/1‑139 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                     <  1:676074/51‑1 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                     >  2:676074/1‑51 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAc  >  1:996349/1‑134 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAc  <  2:996349/134‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
CTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAC  >  NC_000913/1850670‑1850935

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