Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,751 T→C 10.1% 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,7510TC10.1% 129.1 / 4.4 56intergenic (+31/‑208)serS/dmsAserine‑‑tRNA ligase/dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (4/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

GTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCT  >  NC_000913/940610‑940892
                                                                                                                                             |                                                                                                                                             
gTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTg                                                                                                                                              >  2:218745/1‑143 (MQ=255)
 ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGc                                                                                                                                             >  1:72918/1‑143 (MQ=255)
 gtGGGCGTACGCTGGGTGCAGTGATGGAATACTGTCAGCGGGCTGATGGTCGTAGTGAAGTACCAGAAGTTCTGCGTCCGTATGTGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGc                                                                                                                                             <  2:250734/142‑1 (MQ=255)
   ggTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCGGATGGTCGTATTGAAGTACTAGATGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGTCTTTTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCgg                                                                                                                                           <  2:167793/143‑1 (MQ=255)
       gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATCCCCACTTTTTCTGAATCTAAAAAGCGCCTGCGAgcgc                                                                                                                                       >  1:2365/1‑143 (MQ=255)
            cTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATCGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATATCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttt                                                                                                                                  <  1:652100/143‑1 (MQ=255)
             tGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttttt                                                                                                                                 >  2:301449/1‑143 (MQ=255)
                      tAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                                                                        >  1:462817/1‑143 (MQ=255)
                         aagcgtaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTg                                                                                                                     <  2:442237/137‑1 (MQ=255)
                          ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAACAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGa                                                                                                                    >  1:25264/1‑143 (MQ=255)
                          ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGAATGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGa                                                                                                                    >  2:291469/1‑143 (MQ=255)
                           gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGGGCCTGAGGGCGCTTTTTTTGTCTCCCGTTGAt                                                                                                                   >  2:281809/1‑143 (MQ=255)
                             aaaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAc                                                                                                                 >  1:510453/1‑143 (MQ=255)
                                       cagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                       >  2:526310/1‑143 (MQ=255)
                                       cagGCTGATCGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                       <  1:211576/143‑1 (MQ=255)
                                            tGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                                                                  <  2:183950/143‑1 (MQ=255)
                                             gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATTTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                 <  2:504438/143‑1 (MQ=255)
                                               tGGTCGTACTGAAGTACCAGAAGCTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                               <  2:659421/143‑1 (MQ=255)
                                                 gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                             <  1:194734/143‑1 (MQ=255)
                                                  tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAGTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGGCTCCCTTTGATACCGAACAATAATTACGCCTCAc                                                                                            >  2:116255/1‑143 (MQ=255)
                                                            ggaCCAGGAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTTCTCCTCACTTACACGTaa                                                                                  <  2:247450/141‑1 (MQ=255)
                                                              aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAta                                                                                <  1:301449/143‑1 (MQ=255)
                                                                cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtact                                                                              >  2:26034/1‑143 (MQ=255)
                                                                 aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtacta                                                                             <  1:2156/143‑1 (MQ=255)
                                                                   aaGTTCTGCGTCCGTATATGAACGTACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTTAATAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATactact                                                                           <  1:588303/143‑1 (MQ=255)
                                                                               cGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACAACTCTCGAGTGaaac                                                               >  2:257448/1‑142 (MQ=255)
                                                                                   tatGAACGGACTGGAATATATTGGCGAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTGTGTGTCTCCCCCTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTa                                                           <  1:484109/143‑1 (MQ=255)
                                                                                    aggaccGGACTGGACAATATTGGCTAATACCCAATTTTTCTGAATCTCAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                          <  1:281809/138‑1 (MQ=255)
                                                                                           gACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGAGCGCGCTTTTTTTGTCTCCTTTTGATACCCAACAATAATTACTCCTCACTTACACGGAATACTACTTTCGAGTGAAAATCTGCCTAtcta                                                   >  1:659379/1‑142 (MQ=255)
                                                                                              tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt                                                >  1:229552/1‑143 (MQ=255)
                                                                                              tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTATCTCCCTTTGATACCAAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAACTATCTCttt                                                >  2:7404/1‑143 (MQ=255)
                                                                                                gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAAAATAATTACTCCTCACTTACACGTAATACTAATTTCGAGTGAAAGTCTGGCTATCTCTTTGa                                              >  1:267070/1‑143 (MQ=255)
                                                                                                       ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                       >  1:434887/1‑143 (MQ=255)
                                                                                                               tACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGGGAAAATCTACCGATCTCTTTGATTTTCAAATTAGTCg                               >  1:218284/1‑143 (MQ=255)
                                                                                                                 cacAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAACTTATTCgtt                             >  1:307245/3‑141 (MQ=255)
                                                                                                                  ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTACTACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATg                            >  2:111287/1‑143 (MQ=255)
                                                                                                                    aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGCACAATAATTACTCCTCACTTACACGTAATACTACTTACGAGTGAAAATCTACCTATCTCTTTGAATTCAAACTTATTCGATGta                          >  2:257510/1‑143 (MQ=255)
                                                                                                                    aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGta                          <  2:20494/143‑1 (MQ=255)
                                                                                                                       ttttCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACATAACAATAATTACTCCTCACTTACACTTAATACTACTTTCGAGTGAAAATCTACCTATATCTTTGATTTTCAAATTATTCGATGTATAc                       >  2:534619/1‑143 (MQ=255)
                                                                                                                           cTGACTCTAACAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCTCTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGc                   <  2:647615/143‑1 (MQ=255)
                                                                                                                             gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                 <  2:480311/143‑1 (MQ=255)
                                                                                                                             gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                 <  2:462817/143‑1 (MQ=255)
                                                                                                                                caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGCTGTATGCAAGCCtata              >  1:408513/4‑143 (MQ=255)
                                                                                                                                 cTAACAAGCGGCTGCGGTCGCTTTTTTTGTCTCCCTGTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatat             <  1:617919/143‑1 (MQ=255)
                                                                                                                                   aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAg           >  2:299562/1‑143 (MQ=255)
                                                                                                                                   aaaaaGCGCCCGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATTTATACAAGCCTATATAg           >  2:86927/1‑143 (MQ=255)
                                                                                                                                   aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAGTCTACCTATCTCTTTGATGTTCAAATTATTCGATGTATACAAGCGgatatgg           >  1:1017/1‑141 (MQ=255)
                                                                                                                                    aaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGc          <  1:218745/143‑1 (MQ=255)
                                                                                                                                       aGCTCCTGCGGGCGCTTTTGTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATCGCGaa       <  2:550057/143‑1 (MQ=255)
                                                                                                                                       aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa       <  2:434887/143‑1 (MQ=255)
                                                                                                                                       aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa       <  2:563699/143‑1 (MQ=255)
                                                                                                                                       aGCGCCTGCCGGCGCTTTTGTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa       <  1:302225/143‑1 (MQ=255)
                                                                                                                                       aGCCCCTGCTGGTGCTCTTTCTGTCTCCCTTTGATACCGAACAATACCTACTCCTCACTTACACGTAATACTACTGTCGAGTGAAAGTCTACCTATCTCATTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa       <  1:295588/143‑1 (MQ=255)
                                                                                                                                        gcgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc      >  1:383278/1‑143 (MQ=255)
                                                                                                                                           ccTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACATATCTCTTTGATTTTCAAAGTATTCGATGTATACAAGGCTATATAGCGAACTGc   >  1:44176/1‑143 (MQ=255)
                                                                                                                                            cTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAGCAATAATTAATCCTCCCTTACACGGAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCt  <  1:656149/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCT  >  NC_000913/940610‑940892

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