Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 239,412 C→T intergenic (‑676/+7) ykfM ← / ← yafV uncharacterized protein YkfM/2‑oxoglutaramate amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913239,4120CT97.1% 107.8 / NA 34intergenic (‑676/+7)ykfM/yafVuncharacterized protein YkfM/2‑oxoglutaramate amidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/12);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCG  >  NC_000913/239185‑239657
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
tGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACt                                                                                                                                                                                                                                >  1:36621/1‑251 (MQ=255)
  gTAATGAATGACGCCGAATACAAAATAAGTAGAGGAAACGCAGCAGCATTATTTTCAGGGATAAGGCACATTGCTATTAATCTTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCACCACTCACCACAACCTAAtctcc                                                                                                                                                                                                                              >  2:35438/1‑250 (MQ=18)
       gaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAATTCATCCGc                                                                                                                                                                                                                         >  1:224397/1‑251 (MQ=255)
                 gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGCCTTCGGAGTAAATACCCGGCTCTCCCCGTTCTTTAAACACATCCCCACTCCACACCACCTAAAATCCACCCCATCCTGcccat                                                                                                                                                                                                              >  2:101000/1‑248 (MQ=14)
                         aaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAAGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGaa                                                                                                                                                                                                        >  1:230781/1‑250 (MQ=255)
                                     aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCtata                                                                                                                                                                                           <  1:111944/251‑1 (MQ=255)
                                     aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCtat                                                                                                                                                                                            >  2:125476/1‑250 (MQ=255)
                                        cGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTc                                                                                                                                                                                        <  1:73544/251‑1 (MQ=255)
                                                           gggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGAc                                                                                                                                                                     <  1:181531/251‑1 (MQ=255)
                                                                            tATTAATATTTTGACGAATGAGAAGGTATTCACGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACATAAACTCATCCGCATCCTGCCATGCCGGCAACTTTTCTCTATATTTCCGCAATGCTGCCATCGACAGCTCCGCATCAATgcg                                                                                                                                                    >  2:2413/1‑251 (MQ=255)
                                                                                   atTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGcc                                                                                                                                             >  1:25756/1‑251 (MQ=255)
                                                                                      ttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTg                                                                                                                                           <  1:110880/250‑1 (MQ=255)
                                                                                         aCGAATGAGAAGGTATTCAAGTCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATCCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATgc                                                                                                                                       <  2:44756/251‑1 (MQ=255)
                                                                                         aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATg                                                                                                                                        >  2:184053/1‑250 (MQ=255)
                                                                                          cGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGCAAAGCAGCCATGGACAGAAACTACCTGGAGTCAGTCCTTGTGGGGAGCGGGCTTTAGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCAGCATCAATGCGCGCTGCCTGATgcg                                                                                                                                      >  1:241335/1‑251 (MQ=38)
                                                                                              tGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCACTGGACAGAAACTACCTGGCGTCCGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCACTGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTc                                                                                                                                  >  1:251600/1‑251 (MQ=255)
                                                                                                     gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTccc                                                                                                                                                                                      >  1:24561/1‑192 (MQ=255)
                                                                                                     gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTccc                                                                                                                                                                                      <  2:24561/192‑1 (MQ=255)
                                                                                                          cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGa                                                                                                                       <  1:92731/250‑1 (MQ=255)
                                                                                                                      aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTg                                                                                                          <  1:19815/251‑1 (MQ=255)
                                                                                                                           gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTACTTGTGGGGAGCGGGCTTTCGTAGTAATACAACGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTCTATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGAt                                                                                                     >  1:112049/1‑251 (MQ=255)
                                                                                                                                      aaggcagcaATGGACAGAACCTCCCTGGCGGCAGTCCTTGTGGGGACCGGGCGTCCGAAGTAATACCCTGACTCTCCCCGTCCTTAACCACAACCCCCCCTCCCATCAACCTAACCTCACCGGCATCCTGCCATGCCGGAAACTTGTCTCTATATTCCCGCAAGGTTGCCACGGCCAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGGTATATCTCGCCTTGCGGATTAATCTCCCGGc                                                                                         <  1:180937/248‑1 (MQ=255)
                                                                                                                                      aaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCGCCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCgg                                                                                          >  1:34771/1‑251 (MQ=255)
                                                                                                                                                            cccTGGCGTCAGGCCGTGTGGGGAGCGGGCTTTCGTAGGAATCCCCCGACTCTCCCCGTCCTTAAACACAATCCCCCCTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCa                                                                      <  2:216934/248‑1 (MQ=255)
                                                                                                                                                              cTGGCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCg                                                                  >  1:174000/1‑251 (MQ=255)
                                                                                                                                                                ggCGTCAGTCCTTGTGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGtt                                                                >  1:35796/1‑251 (MQ=255)
                                                                                                                                                                               ggggAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGAc                                                 <  2:25756/251‑1 (MQ=255)
                                                                                                                                                                                  gAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGACCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGAAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCCTGCGGATTAATCCCCCGGATGACCCCGCGATAATGGCAGCCGTTGCCCTCGCTGCCGAcgcg                                              >  2:203706/1‑251 (MQ=255)
                                                                                                                                                                                                       aCCCCCACTCTCACCGTCCTTAAACACAATCCCCACTCACCACACCCTAAACTCATCCGCATCCTGCCCTGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCCCCACAtccc                          >  2:133844/1‑248 (MQ=255)
                                                                                                                                                                                                           cGACTCTCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTg                      <  1:125476/251‑1 (MQ=255)
                                                                                                                                                                                                           cGACTCTCCCCGGCCTTAAACACAATCCCCACTCACCACAACCTAAACCCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCGGATGCGCGTCGGCAGTAGCGATAATCTCCCCTGGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCGGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTg                      <  1:35438/251‑1 (MQ=255)
                                                                                                                                                                                                              ctctCCCCGTCCTTAAACACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGAtt                   >  1:36764/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaaCACAATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCggg     >  1:152676/1‑251 (MQ=255)
                                                                                                                                                                                                                                cacaATCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCg  >  1:232539/1‑251 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
TGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCGTTGCCTGATGCGCGTCGGCAGTAGCGATAATCTCTCCTTGCGGATTAATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCG  >  NC_000913/239185‑239657

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