Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,640,924 A→G intergenic (+153/‑18) ytjE → / → yjtD protein YtjE/putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,640,9240AG97.1% 109.3 / NA 35intergenic (+153/‑18)ytjE/yjtDprotein YtjE/putative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/14);  minor base T (0/1);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCG  >  NC_000913/4640679‑4641155
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
cTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCa                                                                                                                                                                                                                                      >  1:77121/1‑249 (MQ=255)
        atcatACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTgg                                                                                                                                                                                                                            >  1:371853/1‑251 (MQ=255)
         tcatACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCCACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTACGCTAAGCAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGc                                                                                                                                                                                                                           >  2:453603/1‑251 (MQ=255)
                aCTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGc                                                                                                                                                                                                                    <  1:262407/251‑1 (MQ=255)
                      aaaCGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATaa                                                                                                                                                                                                              <  2:77121/251‑1 (MQ=255)
                       aaCGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAAc                                                                                                                                                                                                             <  1:31616/251‑1 (MQ=255)
                           tggtgAATCCAATTGTCGAGATTTATTTTTTATAAAATTAACCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTCCAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTATAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGatt                                                                                                                                                                                                         >  1:422689/1‑251 (MQ=255)
                                           cGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCg                                                                                                                                                                                         >  1:437644/1‑251 (MQ=255)
                                                            aaaaTTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATAtt                                                                                                                                                                        >  2:99606/1‑251 (MQ=255)
                                                                 tATCCTAAGGAAACAGAAGGATATATGGCATTATTTAACAATACAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTGAGGATTGTGAAATTCGCACCGCGAAATAAGTAGACAGTGATCACCCTTTTCGCGGTTATTTGTTCAAGAAGAGTGGCAAAATGCGTAGAAGGATTATTCTGGGCTCACCCGCCAGAGCAGAAAATATTgggg                                                                                                                                                                    <  2:246658/250‑1 (MQ=255)
                                                                     cTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCACCCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcg                                                                                                                                                               >  2:397936/1‑251 (MQ=255)
                                                                             aCAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGc                                                                                                                                                       <  1:466974/251‑1 (MQ=255)
                                                                                     atttGTAGCATTTTTGAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTCCTGCGACTTACGTTAAGGATTTGTAAATTCGCCCCGCGTAATAAGTTGACATTGATCACCCTGGTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTTTAACGATTATTCTGGTCGCGCCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGaaaa                                                                                                                                               <  2:39995/248‑1 (MQ=255)
                                                                                                     aacaatTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTgggg                                                                                                                                                                    >  1:56576/1‑214 (MQ=255)
                                                                                                     aacaatTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTGGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTgggg                                                                                                                                                                    <  2:56576/214‑1 (MQ=255)
                                                                                                            cAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCgg                                                                                                                        <  1:303594/251‑1 (MQ=255)
                                                                                                              aCCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGAt                                                                                                                      >  1:19796/1‑251 (MQ=255)
                                                                                                                           tCAGGAACTAATTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTCAGGATTTGTAAATTCGCACCGCGTAATAAGTTGCCAGTGATCACCCTGTTCGCGGATATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGACACTATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAgac                                                                                                           >  2:438371/1‑247 (MQ=255)
                                                                                                                                       ttAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTAAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGa                                                                                              <  2:371853/250‑1 (MQ=255)
                                                                                                                                          gCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGagcc                                                                                           >  2:150182/1‑250 (MQ=255)
                                                                                                                                           ccTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCCCCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGagccag                                                                                         <  1:150182/251‑1 (MQ=255)
                                                                                                                                               ttttAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGTGGCAGCGGCGCGGGCACTGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAggaccccc                                                                                                  >  2:399484/1‑232 (MQ=255)
                                                                                                                                               ttttAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACAcc                                                                                                  <  1:399484/238‑1 (MQ=255)
                                                                                                                                                          gTAAAGGGCTTTTTCTGAGCCTTCCGATAAGGATTTGCGAATTCGCACCGCGTAATAAGTCGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAACATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGGAAATATTGGGGCAGCGGCGCGGGCGATGAAAACGATGGGTTATCGCGTACTGCGGATTGTCGTTAGCCAGGCACACCTGGAGCAAGCCACCCGCTggggcg                                                                          >  2:82576/1‑248 (MQ=255)
                                                                                                                                                                         tGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGcgcg                                                                                                                                                          <  2:260945/156‑1 (MQ=255)
                                                                                                                                                                         tGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGcgcg                                                                                                                                                          >  1:260945/1‑156 (MQ=255)
                                                                                                                                                                           cGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCACTATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGCTTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGGGat                                                         >  2:289724/1‑251 (MQ=255)
                                                                                                                                                                           cGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCCGCGgcgc                                                                                                                                                           >  2:328874/1‑153 (MQ=255)
                                                                                                                                                                           cGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGgcgc                                                                                                                                                           <  1:328874/153‑1 (MQ=255)
                                                                                                                                                                                         gatgtGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAatat                                           <  1:99606/247‑1 (MQ=255)
                                                                                                                                                                                          aTTTGTAAATTCGCTCCGCGTAATAAGCTGATTGGGATCTCCCTGTTCGCGGTTTTTTGGTCAAGAAGAGTGGAAATTGGCGTTTAACGATTATGTTGGGCGCGCTCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAAGGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGAGAATATTGATAATAtt                                          <  2:324651/251‑1 (MQ=255)
                                                                                                                                                                                              gTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGGGATATTATTGATAATAttca                                        >  2:198596/1‑247 (MQ=255)
                                                                                                                                                                                                                        gACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCACTATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGCACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCACCCCCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGtt             >  2:214135/1‑251 (MQ=255)
                                                                                                                                                                                                                              gATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACg        >  2:403231/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cctgTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATgacg  >  2:339142/4‑247 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       
CTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCG  >  NC_000913/4640679‑4641155

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: