Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,357,202 C→T intergenic (‑92/+221) sapA ← / ← ymjA putative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2020CT93.3% 29.8 / ‑5.8 15intergenic (‑92/+221)sapA/ymjAputative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (4/10);  total (4/11)
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

GACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGT  >  NC_000913/1356991‑1357408
                                                                                                                                                                                                                   |                                                                                                                                                                                                              
gACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTTCCGTTTGTTTTCAGGGGTCTTAATGAGGGCCGTGGGTTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                                                                                                           <  1:368210/249‑1 (MQ=255)
                     gTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTTGTTGGTCGTTGAACAGGTGACGTCGCCAtn                                                                                                                                                                             <  1:112859/226‑2 (MQ=255)
                             tatCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACTATAATGCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTTCCGTTTGTGTTCAGGGGTCGTAATGAGGACTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTc                                                                                                                                                                        <  1:74429/223‑1 (MQ=255)
                                                               ttccTGACCACTCCCAAGTCCAACAATCAACAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAgg                                                                                                                                                                                          <  1:261025/169‑1 (MQ=255)
                                                                 ccTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                                                                   <  1:3473/176‑1 (MQ=255)
                                                                    gACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGGCGGCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                                                                                                      <  1:380607/250‑1 (MQ=255)
                                                                         cTCACAAGTCCGGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGAGTTTATGTAAAGAGGCTATCGTACTAATATTTAATGACATTTGCCATTACCGTGTGTGTTCAGGGGTCGTAATGAGGGCTgtgtn                                                                                                                                                                                                           <  1:53446/144‑2 (MQ=255)
                                                                                                    naCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTg                                                                                                                                                                                        >  1:92024/2‑136 (MQ=255)
                                                                                                    gaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGCGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGc                                                                     >  1:413403/1‑251 (MQ=255)
                                                                                                                                                                       gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTGAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGt  <  1:157241/251‑1 (MQ=255)
                                                                                                                                                                          aTTTGNCCTTGNCNNTTGTGTTCAGGGGTCGTAANGAGGGCTGTGTGGTTNGTCGNTGAACNGGNgn                                                                                                                                                                                       >  2:169727/1‑66 (MQ=255)
                                                                                                                                                                                     ccGTTTGTGGTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGATCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCGCTATAACTATAACGCTCAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCgtg                                                             <  1:20282/178‑1 (MQ=255)
                                                                                                                                                                                         ttGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataan                                                                                                                                           <  1:92503/96‑2 (MQ=255)
                                                                                                                                                                                           ntgtTNAGGGGTCGTAATGAGNGCTGTGTGGTTCGTCGTTANNCNGGTGANGTCGc                                                                                                                                                                                 <  2:192891/55‑1 (MQ=255)
                                                                                                                                                                                           gtgtTCAGGGGACGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                                                                                                                 >  1:192891/1‑56 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                              
GACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGT  >  NC_000913/1356991‑1357408

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