Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,778,119 C→T intergenic (+337/+27) ypjF → / ← ypjA toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter YpjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,778,1190CT87.5% 55.9 / ‑1.6 26intergenic (+337/+27)ypjF/ypjAtoxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter YpjA
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (13/8);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACAAATTCTTTGCGCACAG  >  NC_000913/2777882‑2778360
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
tttCTTTTAGGACAAACAATGGGCCCAAAGTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGATACACAACCAATTCAGAATCCGCCAGCCCACCAATCATGATTGGAAGGAGAAAGAACAACACCAACAAAACCAGGATGTTAGAAGTCTCAGCAAAACACCGACCAGACGGAGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGATCATGGTAACATGCAAGCCAACAGACGTggact                                                                                                                                                                                                                                      >  2:3963/1‑248 (MQ=255)
               acaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTc                                                                                                                                                                                                                       >  1:16437/1‑251 (MQ=255)
                         cAAAATGCACATTGAACTGATGATTACAATCAATTTAACAACAAATTACACAACCCATTCAGACTCCGCCAGCACACCAATCATGATTGGACGGAGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGACTCAGCAAAACACCGACCAGACGATGAGGAGACATAAAAGGATACGCAAAGTAGCACCGGCTCCTGGTAACATGAAAACCACCAGATGTGGTCTTTTCCGTTTATGGTCAAAACgccc                                                                                                                                                                                                              >  2:321262/1‑247 (MQ=255)
                                                       caatTTAACAACAAGTTACACAACCAAATGAAACTCCGCCAGCACAACAATAATGATTGGACGATGTAAGGACAACACCAATAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGATGGTGAGGAGACATAAAAGGATACGCAAAGTAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGTGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTCCCGCGTTCCACgggg                                                                                                                                                                                <  2:281710/250‑1 (MQ=255)
                                                         atTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGAGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGAtt                                                                                                                                                                             >  1:316532/1‑251 (MQ=255)
                                                             aacaacAAGTTACACAACCAATTCATACTCCTCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCcac                                                                                                                                                                         <  1:162917/251‑1 (MQ=255)
                                                                    aGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGGAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGc                                                                                                                                                                   >  1:210930/1‑250 (MQ=255)
                                                                               cAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCtat                                                                                                                                                        >  1:241788/1‑250 (MQ=255)
                                                                                                                 ggACGGTATAAGGACAACACCAACAAAAACAGGAAGATATAAGTCTCAGCAAATCCCCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGCGCCGCGGCTCCTGGTAACATGAATGCCTCCAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCTCGTTCCACGGGTATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACTTActgctg                                                                                                                     <  2:66752/251‑1 (MQ=37)
                                                                                                                   aCGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACGCCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGa                                                                                                                    >  1:194452/1‑250 (MQ=255)
                                                                                                                       tgtAAGGACAACACCAACAAAAACAGGAAGTTAGACGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATgcg                                                                                                               >  2:77451/1‑249 (MQ=255)
                                                                                                                       agaaaggtccactccAACAAAAACAGGAATATAGAAGGCCAAGTAAAACTCCGACCAGACTGTGAGGAGACATAAAAGGATTCGTAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGTTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCTCGGGGATTCCACACCGGCACTATGGCTATCCCCCACCCCAAGATTCCCGCTTAACGTACTGCTGAATGAg                                                                                                               <  2:316532/238‑1 (MQ=255)
                                                                                                                                     caacaaAAACAGGAAGTAAGAAGTCTCAGCAAAACACGGACCAGACGGTGAGGAGACATAAAAGGAGACGCAAAGGAGCCGCGGCACCTGGTAACATGAAAGCCTACAGAGGTGGGCTTTTTCGTGGATTGTCTGAACGACCAGTTCACACCAGCTACCGCGTTCTACGGTGATTCGACACCGGCACTATGGCTATACCCCTCCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgc                                                                                                 <  1:369905/251‑1 (MQ=255)
                                                                                                                                        caaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcctg                                                                                              <  1:224171/251‑1 (MQ=255)
                                                                                                                                                     ttAGAAGTCTCAGCAAAACAAGGACAAGAAGGTGTGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCATACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACAAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGAACTATGGCTATAACACACCCCAAGATGCCCGCTTAAAGACCTGCTGAATGAGGATTTAATACCTGCATGGTATAttaccagt                                                                                 >  2:78461/1‑245 (MQ=255)
                                                                                                                                                         aaGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGTTACGCAAAGGAGCCGCGGCTCCTGGTCACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGTTGGTCAGAACGACCAGTTCACACCAGCTATCGCGTTCCACGGGGAATCCACACCGGCACAAGGGACAAACCCCACCCCAAGATGCCCGCTTAACGTTCTACTGAATGAGGCTTTAATACCTACCTGGTATTTTCCatgacagc                                                                             >  2:195811/1‑244 (MQ=255)
                                                                                                                                                                        acCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTCCCCGACAAATAATTg                                                               >  1:70686/1‑250 (MQ=255)
                                                                                                                                                                                 aCGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTAc                                                     >  1:194696/1‑251 (MQ=255)
                                                                                                                                                                                        ggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTCCATGGTATATTCCACGTCTGCCCGCCAAATAATTGACGAAATTACCGTCACt                                              >  2:119740/1‑251 (MQ=255)
                                                                                                                                                                                                aaaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCtata                                                                                                                                                       >  1:21965/1‑138 (MQ=255)
                                                                                                                                                                                                aaaaGGATACGCAAAGGAGCCGCGGCGCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTGCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTGCCACGGGGATTCCACACCGGCACTATGGCtata                                                                                                                                                       <  2:21965/138‑1 (MQ=25)
                                                                                                                                                                                                 aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGCTGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCCGCCCGCCAAATAATTGACGGAATTACCGTCACTATTCAAttt                                     >  1:326112/1‑251 (MQ=255)
                                                                                                                                                                                                    ggATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACTGTCACTATTCACTTTCAc                                  <  2:194452/251‑1 (MQ=255)
                                                                                                                                                                                                              aGGAGCCGCGGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTatc                         >  1:320894/1‑250 (MQ=255)
                                                                                                                                                                                                                     gcgGCTCCTGGTAACATGAAAGCCTACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACaa                 <  1:326367/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGAAAGCCCACAGAGGGGGTCTTTTTCGGTGATGGGTAGATCGACCCGTTCACACCAGATACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCGCCCCAAGATGCCCGCTTAACGTACTGCTGAAAGAGGCTTTAATACCTGCCTGGTATATTCCCCGTCTGCCCGACAATTAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGTCAAATTCTTTGCGCACAg  <  2:16437/251‑1 (MQ=17)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
TTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACAAATTCTTTGCGCACAG  >  NC_000913/2777882‑2778360

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