Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,719:1 +T 12.5% intergenic (+26/‑141) 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,7191.T12.5% 126.6 / 15.8 40intergenic (+26/‑141)sppA/ansAprotease IV, a signal peptide peptidase/L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base . (23/12);  new base T (3/2);  total (26/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 = 6.33e-01

CCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATC  >  NC_000913/1850577‑1850859
                                                                                                                                               |                                                                                                                                            
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)
                                                                                                                                            aTCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaatatc  <  2:634763/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
CCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATC  >  NC_000913/1850577‑1850859

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