Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 244,775 G→A G150D (GGC→GAC)  yafJ → putative glutamine amidotransferase YafJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913244,7750GA96.7% 95.6 / ‑3.6 31G150D (GGC→GAC) yafJputative glutamine amidotransferase YafJ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (10/19);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

AATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAATCGACTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACA  >  NC_000913/244535‑245010
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           
aaTCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGCCAACATggcg                                                                                                                                                                                                                                   >  2:357588/1‑251 (MQ=255)
                                                       ggAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGttt                                                                                                                                                                            <  2:74936/251‑1 (MQ=255)
                                                          aaaTACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCCGCAGACggggggtttc                                                                                                                                                                          >  2:206532/1‑243 (MQ=255)
                                                              aCTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACAt                                                                                                                                                                     <  1:344430/251‑1 (MQ=255)
                                                                tCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACAtgc                                                                                                                                                                   >  1:66970/1‑251 (MQ=255)
                                                                 cACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACAtgc                                                                                                                                                                   <  1:174643/250‑1 (MQ=255)
                                                                 cACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATGTAAATATATCGCCTCACTGGCGGATTAACTGCGGCAGAAGGGCGTTTTCAACAtg                                                                                                                                                                    <  2:125041/249‑1 (MQ=255)
                                                                           cccGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGAc                                                                                                                                                        <  1:162319/251‑1 (MQ=255)
                                                                              gcgAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCAAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACggg                                                                                                                                                     <  1:194625/251‑1 (MQ=255)
                                                                                gAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGTCGAAGTGCGTTTTTCACATGCTGCTTTCTGACggggg                                                                                                                                                   >  2:69406/1‑249 (MQ=255)
                                                                                     aTGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCTCTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACAGGCTGCTTGCGGACGGGCGCTAt                                                                                                                                               <  2:41900/250‑1 (MQ=255)
                                                                                             gTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCg                                                                                                                                      <  1:39045/251‑1 (MQ=255)
                                                                                                aCTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTAt                                                                                                                                   <  1:360479/251‑1 (MQ=255)
                                                                                                  tGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATGTATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTg                                                                                                                                 >  1:143638/1‑251 (MQ=255)
                                                                                                  tGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGGGTTTTCAAAATGCTGCTTTCGGACGGGCGCTATGTAATGGCggaatg                                                                                                                                 >  2:300919/1‑246 (MQ=255)
                                                                                                          tGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATTAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCggct                                                                                                                          >  2:93120/1‑247 (MQ=255)
                                                                                                                         aCTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCGTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTGCGGACGGGCGCTATGTAATGGCGTAGGGCTCGACTAATTTACGCTGGATCa                                                                                                          <  2:208820/251‑1 (MQ=255)
                                                                                                                           ggACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCGCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATAGATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCTGACGGGCGCTATGTAATGGCGTGTTGCTCGACTAATTTACACTTGATCAcc                                                                                                        <  2:307083/250‑1 (MQ=255)
                                                                                                                             aCGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACccg                                                                                                      <  2:301671/251‑1 (MQ=255)
                                                                                                                                    aCAAATCACTGGAAACCGGCAACGTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATATAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACTCTGGATCACCCGCcgcgcg                                                                                               <  2:249339/251‑1 (MQ=255)
                                                                                                                                          cACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGtt                                                                                          <  2:52683/250‑1 (MQ=255)
                                                                                                                                                               ccccGGTAGGCGACACCGACAGCGAAAAAGCCTTCGGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACCCCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGAt                                                                    <  2:196623/250‑1 (MQ=255)
                                                                                                                                                                   ggTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAAAGGCGTTTTTCTCGACTAATTTACTCCGGGTTACCCGCCCCGCGCCGTTTGGGGTGGCCACGTTGCTGGGTCa                                                                  >  2:9102/1‑249 (MQ=255)
                                                                                                                                                                         cGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTgg                                                          <  2:129418/251‑1 (MQ=255)
                                                                                                                                                                           aaaCCGACAGCGAAAAAGCCTTCTGCGGGCTCCTGCAGAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGAGGAACTGCGGCAGAAGGGCGTTGTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCCCGACTAATTTAGACTGGATCACCCGCCGCGCGCTGTTTGGCTGGGCAACGTTGCTGGATTAGGATGTGGaa                                                        <  2:360761/251‑1 (MQ=255)
                                                                                                                                                                                           aGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAATCGACTTCAGCTCGc                                        >  1:280401/1‑251 (MQ=255)
                                                                                                                                                                                                tCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAATCGACTTCAGCTCGCAGAc                                    <  1:68578/250‑1 (MQ=255)
                                                                                                                                                                                                       gCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAATCGACTTCAGCCCGCAGACCACACCg                             >  2:334799/1‑250 (MQ=255)
                                                                                                                                                                                                          ccTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGGGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGGGGAAATCGACTTCAGCTCGCAGACCACACCGAatg                         >  2:139498/1‑251 (MQ=255)
                                                                                                                                                                                                                     tAACGCAGCGTTACCCGCGCACACCGGACAACATGGCGGCGGTATTTAAATATATCGCCGCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAAGCGACTGCAGCTCGCAGACCACACCGAATGATGTGGTCAc               <  2:222048/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aCCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAATCGACTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAaca  <  1:315477/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           
AATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCGGTAGGCGAAACCGACAGCGAAAAAGCCTTCTGCTGGCTCCTGCATAAATTAACGCAGCGTTACCCGCGCACACCGGGCAACATGGCGGCGGTATTTAAATATATCGCCTCACTGGCGGATGAACTGCGGCAGAAGGGCGTTTTCAACATGCTGCTTTCGGACGGGCGCTATGTAATGGCGTATTGCTCGACTAATTTACACTGGATCACCCGCCGCGCGCCGTTTGGCGTGGCAACGTTGCTGGATCAGGATGTGGAAATCGACTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACA  >  NC_000913/244535‑245010

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