Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,204,395 T→A intergenic (+106/‑201) yehK → / → yehL uncharacterized protein YehK/putative AAA(+) MoxR family ATPase YehL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,204,3950TA95.5% 62.6 / NA 22intergenic (+106/‑201)yehK/yehLuncharacterized protein YehK/putative AAA(+) MoxR family ATPase YehL
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (8/13);  total (8/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 = 1.00e+00

TAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGAT  >  NC_000913/2204156‑2204589
                                                                                                                                                                                                                                               |                                                                                                                                                                                                  
taataaATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTgg                                                                                                                                                                                         >  1:483269/1‑251 (MQ=255)
 aataaATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCGTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTgg                                                                                                                                                                                         <  1:383577/250‑1 (MQ=255)
      atgatgCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACg                                                                                                                                                                                   >  1:436723/1‑251 (MQ=255)
             gTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGc                                                                                                                                                                            <  1:256519/251‑1 (MQ=255)
                  tttATATAAGAGATGCGGGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCAAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCAAAgg                                                                                                                                                                       >  2:435023/1‑251 (MQ=255)
                         aaGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTTTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGTTgtgg                                                                                                                                                                 <  1:1878/250‑1 (MQ=255)
                               tGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAgg                                                                                                                                                                       <  1:418519/238‑1 (MQ=255)
                               tGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAgg                                                                                                                                                                       >  2:418519/1‑238 (MQ=255)
                                gAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGc                                                                                                                                                         <  1:483824/251‑1 (MQ=255)
                                      aTTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAAGAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGGCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGaa                                                                                                                                                    <  2:345182/250‑1 (MQ=255)
                                         gatgaGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAAt                                                                                                                                                 >  1:253174/1‑250 (MQ=255)
                                                                                   gCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttattt                                                                                                       >  1:138258/1‑250 (MQ=255)
                                                                                                gTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGACCCGGCGCAACTTCCCAATGGCCTGGATTATAAAGCTCATTAGCTTAATGGCAATGGTGTCCATCCGTGGTTATAAATCGTGCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTTCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCtt                                                                                         <  1:304306/251‑1 (MQ=255)
                                                                                                          aaaTATATTCAACCAAGAAACAGAATAATATACTCTAAATAAATCAAATTGTTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTTAAAATAGTAAATAAAATAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGatatg                                                                               <  1:377054/250‑2 (MQ=255)
                                                                                                                                 aataatATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgcgc                                                        >  1:342998/1‑251 (MQ=255)
                                                                                                                                       ataCTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTAGCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAAc                                                   <  1:475807/250‑1 (MQ=255)
                                                                                                                                               aataatTCAAATTGGACCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGATGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTACAGATATTTTTCGCTCTAACAACTAATGCGCCAAACATTTATTgt                                          >  1:331328/1‑250 (MQ=255)
                                                                                                                                                           tggtactaTCCGTCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAAGTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTTCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAatat                               <  1:93332/243‑1 (MQ=255)
                                                                                                                                                               ccGAACCGGCGCAACGTCCCAATGGCCTGGATTATAAAACTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTAAAAAACGTTCGTTGGAACACTATGTTGCACAGGCTGGGGTGATAGGCCTGAAAATAGTAAATAAAAGAACCTGCCACACCAGCAGGTTTTTTTAGTTATTGTGATCTGCTTTCCAGATATTTGTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGttt                          <  1:285195/251‑1 (MQ=255)
                                                                                                                                                                     cggcgaatcgccccAATGGCCTGGATTAAAAATCTCATTATCTGACTTGCAACGTGGTCCATCCGTGGTTAGAATTCGTTCGTTTGAACACGATGATGCACAAGCTGTTGTGTTAGGCCTGAAAATAGTAAATAAAAGTACCTGCCTCACCAGCAGGTTTTTTCATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGTCAAACATTTATTGCGCGTAAAATATCGTGTATTTc                     <  1:286526/239‑1 (MQ=255)
                                                                                                                                                                             cGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACa            <  2:223130/251‑1 (MQ=255)
                                                                                                                                                                                       gCCTGGTTTATAAAGCCTATTAACTGAATTGAAACGGGGTTCAGCCGTGGTTATAAATCGTGCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGACAAACATTTATTGCGTGTAAAATATCGTTTATTTCATTAATACATTGCAGGGAt  <  2:76944/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                  
TAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGAT  >  NC_000913/2204156‑2204589

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