Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,715 Δ1 bp 11.4% intergenic (+22/‑145) sppA → / → ansA protease IV, a signal peptide peptidase/L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,850,7150A.11.4% 108.1 / 11.1 44intergenic (+22/‑145)sppA/ansAprotease IV, a signal peptide peptidase/L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base A (26/12);  new base . (3/2);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCA  >  NC_000913/1850577‑1850849
                                                                                                                                          |                                                                                                                                       
ccTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACGTAGTGGCCGACAGAtc                                                                                                                                     >  1:182622/1‑143 (MQ=255)
 cTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcg                                                                                                                                   >  1:634763/1‑143 (MQ=255)
   gCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTACTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGGCGTGTACTGAGTGGCCGACCGAtcgtc                                                                                                                                 >  2:580207/1‑143 (MQ=255)
    cACCGCTGGCCTCCGTAGCCTCTACCGCTAAAAGTGGAAGTGACAAGCTGGCCGCGTTTAATGACCCCCAACCCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtcg                                                                                                                                <  1:559718/143‑1 (MQ=255)
      tcGCTGGGCGCGGTAGCCTGTAAGGTTAAAAGTGAAAGTGTCAAGCCGGCCGCGTTTAATGAGCCACAATACCGGTATGCGTTTTGCCTGACCTGCGCCAACAGGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGc                                                                                                                              <  2:23210/142‑1 (MQ=255)
      ccGCTGGCCTCGGTAGCCTCTACCGTTACAAGTGAAAGTGACAGGCTGGCCGCGTTGAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGc                                                                                                                              <  2:614290/143‑1 (MQ=255)
      ccGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGGGACAAGCTAGCCGCGTTGAATGACCAACAAAACCGTTATGCGGTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTTGCCGACAGCtagtcgcg                                                                                                                              >  1:515511/1‑141 (MQ=255)
        gCTGGCCTCGGGAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTCAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCca                                                                                                                            <  2:369877/143‑1 (MQ=255)
                   tAGCTTCTACGGTTAAAAGCGCAAGGGACACTCCGGCCGCGTCTAATCCCCCACCACACCGTTCTGCGTTGCGCCTGACCCGCGCCAACAGGCGGTAACTCTTGTACTGAGTGGCCGACATAGCGTCGGCCACattattcgtt                                                                                                                 <  2:634074/143‑5 (MQ=255)
                       ctctACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCCTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGt                                                                                                             >  1:607582/1‑143 (MQ=255)
                              gTTAAAAGTGAAAGTGACAAGCTGGCAGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGCCTTGTACTGAGTGGCAGACAGATCGTCGGCCACATTATTATTTACGTCGACGaa                                                                                                      >  2:589948/1‑143 (MQ=255)
                                 aaaaGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATcc                                                                                                   >  2:917710/1‑143 (MQ=255)
                                                aaGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACCTTATTTTTTACGTCGACGAATCCGCGTCCCGCCGTTTc                                                                                    >  2:366068/1‑143 (MQ=255)
                                                aaGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTc                                                                                    >  1:780860/1‑143 (MQ=255)
                                                  gCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                                                                  <  2:481136/143‑1 (MQ=255)
                                                    tGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTGTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGccc                                                                                <  1:480908/143‑1 (MQ=255)
                                                       ccGCGTTTAATGCCCCCCATATCCGTTATGCGCTTTGCCGTCCCTGCGCCCACATGCGTTAAGTCTTGTTCTGAGTCGCCGACAGATCTTCGGCCACATTATTTCTTACGTCGACGCATCCTCTTCCCGCTGTTTCGCCCCat                                                                             <  2:240887/143‑1 (MQ=255)
                                                       ccGCGCTTAATGCCCCGCAAGCCCGTGATGCGTTTTGCCTGACCTGCGCCCACATGCGTTCGGCCCTGCACTGCCGGGCCCACAGTTCGTCGGCTCCATTAGTTTTTAGGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCat                                                                             <  2:873656/143‑1 (MQ=255)
                                                           gTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGTCAAGATGCGTTAAGTCTTGTATTGAGTGGCTGACAGATTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAttt                                                                         >  1:623545/1‑143 (MQ=255)
                                                           gTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAttt                                                                         >  2:836492/1‑143 (MQ=255)
                                                                 tGACCCACAAAGCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTAt                                                                   <  1:341874/143‑1 (MQ=255)
                                                                     ccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGTCGACGCCTCCACCCGCTGTTTCGCCCCATATTTCCTTATCCGc                                                               >  2:412634/1‑143 (MQ=255)
                                                                     ccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAAGCCCCTTCCCGCCGATTCGCCACATATTTCCTCCTCCGc                                                               >  2:736976/1‑143 (MQ=255)
                                                                       acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGcgt                                                             >  1:73027/1‑141 (MQ=255)
                                                                          aaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATCATCTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATCTCCTTATCCGCAGAGt                                                          >  1:789407/1‑143 (MQ=255)
                                                                           aaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATAGTTCCTTGTCTGCATGGTa                                                         >  1:491101/1‑143 (MQ=255)
                                                                               gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAg                                                     <  1:917710/143‑1 (MQ=255)
                                                                                                     cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGt                               >  1:89455/1‑143 (MQ=255)
                                                                                                     cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGt                               >  1:818229/1‑143 (MQ=255)
                                                                                                      gcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGta                              <  1:621024/143‑1 (MQ=255)
                                                                                                          aCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGtataga                          >  1:441371/1‑141 (MQ=255)
                                                                                                                 ttAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGc                   >  1:871957/1‑143 (MQ=255)
                                                                                                                 ttAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGCGCGCTCCCCCTGCCTCACGTGTATACTTTTGc                   >  1:880669/1‑143 (MQ=255)
                                                                                                                 ttAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGCTGCGCTCCCCCTGCCTCACGTACCTACTTTTGc                   >  1:887468/1‑143 (MQ=255)
                                                                                                                      tCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTCCGCTCCCCCTCCCCCACGTACATACTTTTGCTCttt              >  2:410634/1‑143 (MQ=255)
                                                                                                                        ttGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCCCACGTACGTCCTTTTGCTCTTTCg            >  2:452184/1‑143 (MQ=255)
                                                                                                                              ggcgTGGCCGACAGATCGGCGGCCACATTATGTTTTACGTGGACGCATCCTCTTCCCGCTGTTTCGCCCCACATTACCTTGTCCGCGTAGTATCAGGTGCGCTCCCCCTGCCTCCCGTATATACTTTTGCTCTTTCGATAtca      <  1:353092/140‑1 (MQ=255)
                                                                                                                               atgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGAAGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCCCTCCCCCCGCCTCACGTATACACTTTTGCTCTTTCGATAtcat     >  2:874419/3‑143 (MQ=255)
                                                                                                                                aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCCTATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATattatt    >  1:921540/1‑143 (MQ=255)
                                                                                                                                 gTGGCCGACGATCTGTCGGCCACATTATGTTTTACGTCGACGAATCCTCTTGCCCCTGTTTCGCCCCATGTTTCCTTATCCGCATAGTATCAGGTGCGCCCCCCCTACCGCACGCATATACTTTTGCTCTTTCGACAtcctcc   >  2:182448/1‑139 (MQ=255)
                                                                                                                                  tGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAGGGCTCTTCCCGCTGTTTCGCCCTATATTGCCTTATCCGAATAGTATCAGGTGCGCTCCCACTGCCTCACGTATATCCTTTTGCTCTTTCGCTAtcattca  >  1:190817/1‑143 (MQ=255)
                                                                                                                                  tGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATGtcattgc  >  1:190503/1‑141 (MQ=255)
                                                                                                                                  tGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca  <  2:341225/143‑1 (MQ=255)
                                                                                                                                  tGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGGTTCGCCCCATATTTCCTTATCCGTATAGCATCAGGTGCGCTCCGCCTGCCTCCGCCATCAACCTTTGCTCTCGCGATGtcattct  >  2:682382/1‑142 (MQ=255)
                                                                                                                                          |                                                                                                                                       
CCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCA  >  NC_000913/1850577‑1850849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: