Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,341 T→C 16.0% intergenic (‑600/+58) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,3410TC16.0% 67.2 / 4.8 25intergenic (‑600/+58)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base C (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCGAAGAGCTACCTTTCACCTTCAGATCTTTAAT  >  NC_000913/4326112‑4326554
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
tAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑Gt                                                                                                                                                                                                                 <  1:2571/238‑1 (MQ=14)
                                             tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGCCTCGCAATCGTCGCATCAGGCATCTTGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GGTGGCCTTATCCGGCTTACCGATCtcaat                                                                                                                                                                  >  2:159839/1‑238 (MQ=14)
                                                      aaCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGt                                                                                                                                                         <  1:100005/238‑1 (MQ=255)
                                                            ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGt                                                                                                                                                         >  1:156144/1‑232 (MQ=255)
                                                                                                        aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATc                                                                                                                                                                       >  1:133941/1‑174 (MQ=255)
                                                                                                        aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATc                                                                                                                                                                       <  2:133941/174‑1 (MQ=255)
                                                                                                            gcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCa                                                                                                  <  1:159839/239‑1 (MQ=255)
                                                                                                                                 cGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCCGGCAACGGTGCCAGATGCCTGATGCGACGCTTGCGCGTCTTATCAGGCTTACCGATCACAATACAATCAGTTCt                                                                                                                                                      <  2:225354/175‑1 (MQ=11)
                                                                                                                                 cGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCCGGCAACGGTGCCAGATGCCTGATGCGACGCTTGCGCGTCTTATCAGGCTTACCGATCACAATACAATCAGTTCt                                                                                                                                                      >  1:225354/1‑175 (MQ=11)
                                                                                                                                              ggATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGc                                                                                                                                                                                >  1:131469/1‑127 (MQ=255)
                                                                                                                                              ggATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGc                                                                                                                                                                                <  2:131469/127‑1 (MQ=255)
                                                                                                                                                   cggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCCGGCGGATGtttt                                                            >  1:125447/1‑238 (MQ=255)
                                                                                                                                                                ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGAt                                                                                                                                                                        >  2:98401/1‑117 (MQ=255)
                                                                                                                                                                ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGAt                                                                                                                                                                        <  1:98401/117‑1 (MQ=255)
                                                                                                                                                                ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATGTAGCAATCGATGTTATGGTCGCCTTCAACCATGCGGACGTTTTTCACTTTGTTGCCa                                              >  2:219550/1‑239 (MQ=255)
                                                                                                                                                                        ccTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCa                                       <  1:219550/238‑1 (MQ=255)
                                                                                                                                                                                 ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTAt                                                                                                                                                                                     <  1:75385/87‑1 (MQ=255)
                                                                                                                                                                                 ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTAt                                                                                                                                                                                     >  2:75385/1‑87 (MQ=255)
                                                                                                                                                                                        gAATTTTGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAAtttt                                         <  1:228232/220‑1 (MQ=255)
                                                                                                                                                                                             tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCg                                 <  1:61634/223‑1 (MQ=255)
                                                                                                                                                                                             tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCg                                 >  2:61634/1‑223 (MQ=255)
                                                                                                                                                                                                  gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGTTGATGGTATTCACCTTGGGGCCAATTTTCGGCATTGaa                               >  2:228232/1‑220 (MQ=38)
                                                                                                                                                                                                    cTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCGAAGAGCTACCtt                     >  2:134337/1‑228 (MQ=255)
                                                                                                                                                                                                    cTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCGAAGAGCTACCtt                     <  1:134337/228‑1 (MQ=255)
                                                                                                                                                                                                             cgcgCAAGCGTCGCATCGGGCATCTGGCC‑‑‑‑‑‑CG‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTAGCCGGCTTACCGATCCCAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCGAAGAGCTACCTTTCACCTTCAGATCTTTAat  <  2:156144/238‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
TAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCA‑‑‑‑‑‑CC‑GTTGCCGGATGCGGCGC‑‑GAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCATCGAAGAGCTACCTTTCACCTTCAGATCTTTAAT  >  NC_000913/4326112‑4326554

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑