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% 165.7 / NA 40intergenic (+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/19);  total (21/19)

CTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCG  >  NC_000913/1850670‑1850937
                                                                                                                                        |                                                                                                                                   
cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                   <  2:1321004/139‑1 (MQ=255)
        cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCa                                                                                                                          <  2:1550067/139‑1 (MQ=255)
                               aCTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGa                                                                                                   >  1:512072/1‑139 (MQ=255)
                                   aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccctgc                                                                                                                               >  1:198839/1‑105 (MQ=255)
                                   aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTcccctgc                                                                                                                               <  2:198839/108‑4 (MQ=255)
                                      ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  2:512072/139‑1 (MQ=255)
                                                  gtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            <  1:136411/95‑1 (MQ=255)
                                                  gtcgGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            >  2:136411/1‑95 (MQ=255)
                                                         acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                 >  1:1579970/1‑115 (MQ=255)
                                                         acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                 <  2:1579970/115‑1 (MQ=255)
                                                             tattTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  1:1435820/1‑116 (MQ=255)
                                                             tattTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  2:1435820/116‑1 (MQ=255)
                                                                     aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            >  1:992939/1‑76 (MQ=255)
                                                                     aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            <  2:992939/76‑1 (MQ=255)
                                                                         cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       <  2:497142/109‑1 (MQ=255)
                                                                         cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       >  1:497142/1‑109 (MQ=255)
                                                                                     ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                   >  1:780197/1‑54 (MQ=255)
                                                                                     ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTccc                                                                                                                                   <  2:780197/54‑1 (MQ=255)
                                                                                         cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       >  2:1521129/1‑93 (MQ=255)
                                                                                         cGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       <  1:1521129/93‑1 (MQ=255)
                                                                                          gCTGTTTCGCCCCATATTTCTTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            <  2:508583/55‑1 (MQ=255)
                                                                                          gCTGTTTCGCCCCATATTTCTTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCt                                                                                                                            >  1:508583/1‑55 (MQ=255)
                                                                                           cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCa                                                                                                                          <  1:1011100/56‑1 (MQ=255)
                                                                                           cTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCa                                                                                                                          >  2:1011100/1‑56 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  1:837577/82‑1 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  2:837577/1‑82 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGaaa                                                                    <  1:224546/106‑1 (MQ=255)
                                                                                               ttCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGaaa                                                                    >  2:224546/1‑106 (MQ=255)
                                                                                                   ccccATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgg                                              >  1:521161/1‑124 (MQ=255)
                                                                                                   ccccATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcgg                                              <  2:521161/124‑1 (MQ=255)
                                                                                                                tATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            >  1:999941/1‑65 (MQ=255)
                                                                                                                tATCCGCATAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt                                                                                            <  2:999941/65‑1 (MQ=255)
                                                                                                                       aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                 <  1:255813/53‑1 (MQ=255)
                                                                                                                       aTAGTATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGata                                                                                                 >  2:255813/1‑53 (MQ=255)
                                                                                                                           tATCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTtata       >  1:908565/1‑139 (MQ=255)
                                                                                                                             tCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       <  2:1635467/57‑1 (MQ=255)
                                                                                                                             tCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       >  1:1635467/1‑57 (MQ=255)
                                                                                                                             tCAGGTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTtatata     >  2:104069/1‑139 (MQ=255)
                                                                                                                                ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                          >  2:71211/1‑83 (MQ=255)
                                                                                                                                ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg                                                          <  1:71211/83‑1 (MQ=255)
                                                                                                                                ggTGCGCTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCg  >  2:371039/1‑139 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       <  2:1314971/51‑1 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat                                                                                       >  1:1314971/1‑51 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCggg                                >  2:1615357/1‑106 (MQ=255)
                                                                                                                                   gcgcTCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCggg                                <  1:1615357/106‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGACCATCGGGATGCAGCGTTCCGAGCAGGGTTATATACCG  >  NC_000913/1850670‑1850937

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: