Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,830,112 G→C 7.4% S218S (TCC→TCG)  astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,1120GC 7.4% 161.1 / 3.6 51S218S (TCC→TCG) astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base C (2/2);  total (25/26)
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

GCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCAGAACGGTGAATAGCCGTG  >  NC_000913/1829965‑1830258
                                                                                                                                                   |                                                                                                                                                  
gcgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTgg                                                                                                                                                    >  8:377395/1‑148 (MQ=255)
     cgcgGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaa                                                                                                                                               <  2:161741/148‑1 (MQ=255)
            gATGTCATACTCAATCGTCGGACCTCCGTTAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGAACCGCGCGGGCATGCTCTGTTTGCGGATGTACCTGTCCGATGACGTACTGGGCTTCCTGTGTTATATAGTg                                                                                                                                         <  4:39435/147‑1 (MQ=255)
            gATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTg                                                                                                                                         <  4:121264/147‑1 (MQ=255)
                       cAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGta                                                                                                                                     >  3:201571/1‑140 (MQ=255)
                        aaGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATa                                                                                                                                                <  3:71435/128‑1 (MQ=255)
                        aaGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATa                                                                                                                                                >  4:71435/1‑128 (MQ=255)
                        aaGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGa                                                                                                                            <  1:40396/148‑1 (MQ=255)
                         aGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGta                                                                                                                                     >  4:161015/1‑138 (MQ=255)
                            gTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGtt                                                                                                                        >  5:264297/1‑148 (MQ=255)
                             tCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTTGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGttt                                                                                                                       <  5:134903/148‑1 (MQ=255)
                                   cACCGTCAAAGATGTCGATATAGGTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCAt                                                                                                                 >  4:40782/1‑148 (MQ=255)
                                      cGTCAAAGATGTCGATATAGTTACGGTAGCGAACACCTTCTTTCTCCAGAACAGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGATTCCTGGGATAAAAAGTGGGTATAGATCGGAGGTTTCGGCATCAg                                                                                                              >  2:70476/1‑148 (MQ=255)
                                      cGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAg                                                                                                              >  4:348640/1‑148 (MQ=255)
                                        tCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGtt                                                                                                            >  8:317680/1‑148 (MQ=255)
                                         cAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                                                                           <  6:360773/148‑1 (MQ=255)
                                           aaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTg                                                                                                         >  4:54002/1‑148 (MQ=255)
                                             gATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCa                                                                                                       <  6:264297/148‑1 (MQ=255)
                                             gATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAGGTCCCTGACCGATGACGTCCTGGGCTTCCGGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                                                                                        >  8:119035/1‑147 (MQ=255)
                                                 tCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaa                                                                                                   <  8:42460/148‑1 (MQ=255)
                                                 gcGATATACTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaa                                                                                                   <  1:318758/147‑1 (MQ=255)
                                                   gATATAGTTACGGTATCGAAAACCTTCTTTCTCGAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTTGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaaa                                                                                                  <  6:120666/147‑1 (MQ=255)
                                                   gATATAGTTACGGTAGCGTAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGTTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGTGATAAAAAGTTGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATtaa                                                                                                  <  4:88487/147‑1 (MQ=255)
                                                       taGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCt                                                                                             >  1:14868/1‑148 (MQ=255)
                                                          ttACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGCGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                                                                                      <  1:345168/136‑1 (MQ=255)
                                                          ttACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGCGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                                                                                      >  2:350454/1‑136 (MQ=255)
                                                          ttACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGCGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                                                                                      >  2:345168/1‑136 (MQ=255)
                                                          ttACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGCGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                                                                                      <  1:350454/136‑1 (MQ=255)
                                                               gTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGTGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCg                                                                                     >  6:366374/1‑148 (MQ=255)
                                                               gTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                                                                           <  5:244517/126‑1 (MQ=255)
                                                               gTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                                                                           >  6:244517/1‑126 (MQ=255)
                                                                 aGCGAAAACCTTCTATCTCCAGCACCGCGCGGGCAGGCGCGGGTGGCGGAGGTACCTGACCGATGACCTCCTGGGCTTCCGGGGATAAATAGTGGGTATAGAACGGATGTTTGGGCCTCCGTTGTGCAATAAATGCGTTATGccatgt                                                                                   >  8:342323/1‑144 (MQ=255)
                                                                  gCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTg                                                                                  <  3:54002/148‑1 (MQ=255)
                                                                     aaaaCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                                                                               >  1:135004/1‑148 (MQ=255)
                                                                       aaCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                                                                               >  8:731/1‑146 (MQ=255)
                                                                             ctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAACGATGACGTCCTGGGATTCCTGGGATAAAAAGTGTTTATAGATCGGATTTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTTCCGCAgaga                                                                         <  1:70476/146‑1 (MQ=255)
                                                                                    cAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcg                                                                <  2:33245/148‑1 (MQ=255)
                                                                                     aGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                                                               <  7:377395/148‑1 (MQ=255)
                                                                                                     gcgcGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCg                                                >  4:251600/1‑147 (MQ=255)
                                                                                                        cgGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGaa                                              >  1:31296/1‑146 (MQ=255)
                                                                                                           tttGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAaag                                           >  8:95122/1‑146 (MQ=255)
                                                                                                             tGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAaagaag                                        <  5:366374/147‑1 (MQ=255)
                                                                                                               cGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCg                                      <  7:317680/147‑1 (MQ=255)
                                                                                                                 gATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGtt                                    >  8:97205/1‑147 (MQ=255)
                                                                                                                           aCCGATGACGTCCTGGGCTTCCTGGGATAAAAAGGGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCtt                         <  6:366737/148‑1 (MQ=255)
                                                                                                                            ccGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCttt                        >  1:161820/1‑148 (MQ=255)
                                                                                                                             cGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTg                       <  3:251600/148‑1 (MQ=255)
                                                                                                                                tGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCa                    <  2:161820/148‑1 (MQ=255)
                                                                                                                                       cTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCAGAACGGt             <  2:31296/148‑1 (MQ=255)
                                                                                                                                          ggCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATAAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGAGTTTACCGAGGCTTTGCCAGAACGGTGaa          >  3:134618/1‑148 (MQ=255)
                                                                                                                                               ccTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCAGAACGGTGAATAg       <  3:40782/146‑1 (MQ=255)
                                                                                                                                                 tGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCAGAACGGTGAATAGCCgt   >  4:24427/1‑148 (MQ=255)
                                                                                                                                                  gggATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCAGAACGGTGAATAGCCgtg  >  8:102898/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGAAAAGAAGCGTTTACCGAGGCTTTGCCAGAACGGTGAATAGCCGTG  >  NC_000913/1829965‑1830258

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