Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,754 G→A 13.6% intergenic (+34/‑205) serS → / → dmsA serine‑‑tRNA ligase/dimethyl sulfoxide reductase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913940,7540GA13.6% 130.6 / 13.7 59intergenic (+34/‑205)serS/dmsAserine‑‑tRNA ligase/dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base G (30/21);  new base A (5/3);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

GTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTA  >  NC_000913/940614‑940893
                                                                                                                                            |                                                                                                                                           
gTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACAGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCggg                                                                                                                                           <  2:344819/143‑1 (MQ=255)
     aCGCTACTTGCAGTAATGGAAAACTATCAGCAGGCTGATGCTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtt                                                                                                                                      <  1:240880/143‑1 (MQ=255)
      cGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttt                                                                                                                                     <  1:665858/143‑1 (MQ=255)
       gCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttt                                                                                                                                    <  1:53563/143‑1 (MQ=255)
             tGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCATTTTTTGtc                                                                                                                              >  1:213618/1‑143 (MQ=255)
              gCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                             >  1:52748/1‑143 (MQ=255)
              gCAGTAATGGAAAACTAGCAGCAGGCTGATGGTCGTACTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                             <  1:555481/143‑1 (MQ=255)
              gCAGTAATGAATAACTATCAGCAGGCTGATGGTCGCATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                             <  1:123236/143‑1 (MQ=255)
               cAGTAAAGGAAAACTATCCGCAGGCTGATGGGCGTATTGGAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATCTATTTGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattca                                                                                                                            <  1:367880/143‑8 (MQ=255)
                aGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                                                           >  1:813095/1‑143 (MQ=255)
                  tAATGGAAAACTATCAGCGGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                                                                         <  2:612193/143‑1 (MQ=255)
                  tAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                                                                         >  2:626299/1‑143 (MQ=255)
                     tGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTg                                                                                                                      >  2:106140/1‑143 (MQ=255)
                        aaaaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                                                                   >  1:574811/1‑143 (MQ=255)
                             tgtcggcagGCTGATGGTCGGATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCGAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                              <  1:355278/138‑1 (MQ=255)
                             tactagcagGCTGATGGTCTTATTTAAGTATCAGATGTTTTGCGTACGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                              <  1:26126/139‑1 (MQ=255)
                             tATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                              >  2:613928/1‑143 (MQ=255)
                                 agcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaat                                                                                                          >  1:848996/1‑143 (MQ=255)
                                   cagGCTGATGGTCGTATTGAAGAACCAGAAGTTCTGCGTCCGTATAGGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                        <  1:520830/143‑1 (MQ=255)
                                    agGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTTAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                                                                       <  1:273001/143‑1 (MQ=255)
                                      gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                                                                     <  2:2242/143‑1 (MQ=255)
                                              tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCAc                                                                                             >  1:724958/1‑143 (MQ=255)
                                                gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATAACTACTCACtt                                                                                           >  1:425264/1‑143 (MQ=255)
                                                 tATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTa                                                                                          >  1:204873/1‑143 (MQ=255)
                                                       aGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTa                                                                                    >  1:848651/1‑143 (MQ=255)
                                                        gTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATCCCGAACCAAAATTACTCCTCACTTACACCTaa                                                                                   >  2:38871/1‑143 (MQ=255)
                                                        gTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                                                                   <  1:549528/143‑1 (MQ=255)
                                                           ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGCATATATTGGCTAATACCCAGTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGATTTTTTTGTCTGCCGTTGATAGCGAACAGTAATTACTCGAGACTTACAAGTaggag                                                                                >  1:37615/1‑139 (MQ=255)
                                                           ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtac                                                                                >  1:757749/1‑143 (MQ=255)
                                                           ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAACAATTACTCCTCACTTACACGTAAtac                                                                                >  2:511030/1‑143 (MQ=255)
                                                             cgcggtccTGCGTCCCTATATCACCGGACTGGAAAATATCCGCCAACCCCCAACTTCACCGAATCTCAACACCGCCTGCGGGCGCTTTTTTTGTCTACCAAAGCCGCCAGCCAATCATTACTCCTCACTTACACGTAAtacta                                                                              <  2:151262/136‑1 (MQ=255)
                                                             aGAAGTTCTGCGTCCGTATATGAACGGACTGGACTATATTGGCTAATACCCACTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAACAATTACTGCTCACTTACACGTAAtacta                                                                              >  2:656796/1‑143 (MQ=255)
                                                                     tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                                                                      >  1:134828/1‑143 (MQ=255)
                                                                              atatGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCt                                                             >  1:540435/1‑143 (MQ=255)
                                                                               tatGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTa                                                            <  1:106140/143‑1 (MQ=255)
                                                                                atGAACGGACTGGAATATATTCGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTTTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                           <  2:315700/143‑1 (MQ=255)
                                                                                         cTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTCTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCtt                                                  <  1:132138/143‑1 (MQ=255)
                                                                                            gAATATATTGGCTAATACCCACTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTCATACCGAACAATAATTACTCCTAACTCACACGTAATGCTACTTTCGAGTGAAGAACTACCTATATCttgtg                                               >  2:11200/1‑140 (MQ=255)
                                                                                              acagatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtt                                             <  1:153449/139‑1 (MQ=255)
                                                                                                atatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGCtttt                                           >  1:783349/1‑143 (MQ=255)
                                                                                                    tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                       <  1:613928/143‑1 (MQ=255)
                                                                                                        tAATACCCAATTTTTCTGATTCTTAAAAGCGTCTGTGGGCGTTTTTTTTGTCTGGCTTGGTGAGCGGACAATGATTGCTCCTCACTTACGCGTAATACGACTGTCGAGTGAAAATCTACCTATCGCTTTGATTTTCAAattat                                   <  1:360255/143‑1 (MQ=255)
                                                                                                             cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGAt                              <  1:626299/143‑1 (MQ=255)
                                                                                                             cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGAt                              >  2:354077/1‑143 (MQ=255)
                                                                                                               cAATTTTTCTGGATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGt                            >  2:798565/1‑143 (MQ=255)
                                                                                                                        tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGCACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGcc                   >  2:593385/1‑143 (MQ=255)
                                                                                                                        tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCTATGTATACAAGcc                   >  2:382595/1‑143 (MQ=255)
                                                                                                                        tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATTTATACAAGcc                   >  1:209085/1‑143 (MQ=255)
                                                                                                                            tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata               >  1:349208/1‑143 (MQ=255)
                                                                                                                            tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAAGAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGCAAATCTAACTAACTATTTGATTTTCACATTATTCGATGTATACAAGCCtata               >  2:564140/1‑143 (MQ=255)
                                                                                                                             cTAACAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatat              <  2:687460/143‑1 (MQ=255)
                                                                                                                             aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATCCAAGCCtatat              >  1:174706/3‑143 (MQ=255)
                                                                                                                             aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatat              >  1:810999/3‑143 (MQ=255)
                                                                                                                               aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAg            >  2:285764/1‑143 (MQ=255)
                                                                                                                                aaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGc           >  1:799294/1‑143 (MQ=255)
                                                                                                                                  aaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                                                                       <  1:651938/49‑1 (MQ=37)
                                                                                                                                  aaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                                                                       >  2:651938/1‑49 (MQ=25)
                                                                                                                                   aGCGGCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGGACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa        <  2:309891/143‑1 (MQ=255)
                                                                                                                                      gcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATTTACCTATCTCTTTACTTTTCAAATTATTCGATGTATACAATCCTATATAGCTAACta     >  2:251268/1‑142 (MQ=255)
                                                                                                                                         tGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAATGCta  >  2:78373/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTA  >  NC_000913/940614‑940893

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