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% 157.1 / NA 38intergenic (+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/20);  total (18/20)

GCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGGTG  >  NC_000913/1850677‑1850940
                                                                                                                                 |                                                                                                                                      
gcgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTGTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTc                                                                                                                              >  1:973778/1‑139 (MQ=255)
                     tGTACTGAGTGGCCGACAGCTCGTCGGCCACATTATTTTTTACGTCGACGAATCGTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         <  1:1295415/139‑1 (MQ=255)
                            aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                                  >  1:862501/1‑139 (MQ=255)
                            aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc                                                                                                  >  1:495100/1‑139 (MQ=255)
                               ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  2:973778/139‑1 (MQ=255)
                               ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  2:333262/139‑1 (MQ=255)
                                      aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc                                                                                                                                       <  2:633228/93‑1 (MQ=255)
                                      aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcc                                                                                                                                       >  1:633228/1‑93 (MQ=255)
                                      aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               >  1:1414115/1‑100 (MQ=255)
                                      aGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               <  2:1414115/100‑1 (MQ=255)
                                                             tACGTCGACGAATCCTCTTCCCGCTGTTTTGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               >  1:1026289/1‑77 (MQ=255)
                                                             tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               <  2:1026289/77‑1 (MQ=255)
                                                                  cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  2:82591/104‑1 (MQ=255)
                                                                  cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               >  1:82591/1‑104 (MQ=255)
                                                                               tCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgggc                                               <  2:862501/139‑1 (MQ=255)
                                                                                cccGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      >  1:94232/1‑52 (MQ=255)
                                                                                cccGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                      <  2:94232/52‑1 (MQ=255)
                                                                                   gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGa                                           >  1:1409025/1‑139 (MQ=255)
                                                                                    cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               <  1:267967/86‑1 (MQ=255)
                                                                                    cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                               >  2:267967/1‑86 (MQ=255)
                                                                                               cATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTc                                                                                                                              <  1:1466444/44‑1 (MQ=255)
                                                                                               cATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTc                                                                                                                              >  2:1466444/1‑44 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               <  1:750423/41‑1 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                               >  2:750423/1‑41 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGtata                                                                                                                       <  1:896210/49‑1 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGtata                                                                                                                       >  2:896210/1‑49 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         >  2:176953/1‑63 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         <  1:176953/63‑1 (MQ=255)
                                                                                                 tatTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAg                             >  1:155507/1‑139 (MQ=255)
                                                                                                    ttCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGt                          <  2:495100/139‑1 (MQ=255)
                                                                                                            ccGCATAGTATCAGGTGTGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa                  <  2:165610/139‑1 (MQ=255)
                                                                                                            ccGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCa                  <  2:1409025/139‑1 (MQ=255)
                                                                                                              gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTg                                                                                                               >  1:440313/1‑44 (MQ=255)
                                                                                                              gCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTg                                                                                                               <  2:440313/44‑1 (MQ=255)
                                                                                                                aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         <  2:197972/48‑1 (MQ=255)
                                                                                                                aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCttt                                                                                                         >  1:197972/1‑48 (MQ=255)
                                                                                                                aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    <  1:1327722/53‑1 (MQ=255)
                                                                                                                aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                    >  2:1327722/1‑53 (MQ=255)
                                                                                                                    tATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGAc                                          >  1:1003123/1‑107 (MQ=255)
                                                                                                                    tATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGAc                                          <  2:1003123/107‑1 (MQ=255)
                                                                                                                       cAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATAc       <  2:155507/139‑1 (MQ=255)
                                                                                                                            gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGgtg  >  1:45592/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                      
GCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCGGTG  >  NC_000913/1850677‑1850940

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