Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,751 T→C 13.4% intergenic (+31/‑208) 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,7510TC13.4% 132.3 / 12.9 60intergenic (+31/‑208)serS/dmsAserine‑‑tRNA ligase/dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base C (5/3);  total (35/25)
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.15e-01

GTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTA  >  NC_000913/940610‑940893
                                                                                                                                             |                                                                                                                                              
gTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTg                                                                                                                                               >  1:221374/1‑143 (MQ=255)
    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                                                                                                       >  2:651938/1‑49 (MQ=25)
                                                                                                                                      aaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                                                                       <  1:651938/49‑1 (MQ=37)
                                                                                                                                       aGCGGCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGGACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa        <  2:309891/143‑1 (MQ=255)
                                                                                                                                          gcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATTTACCTATCTCTTTACTTTTCAAATTATTCGATGTATACAATCCTATATAGCTAACta     >  2:251268/1‑142 (MQ=255)
                                                                                                                                             tGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAATGCta  >  2:78373/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTA  >  NC_000913/940610‑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: