Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,754 G→A 10.5% 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,7540GA10.5% 98.4 / 4.5 47intergenic (+34/‑205)serS/dmsAserine‑‑tRNA ligase/dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base A (3/2);  total (26/21)
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.64e-01

TGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCG  >  NC_000913/940612‑940885
                                                                                                                                              |                                                                                                                                   
tGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCg                                                                                                                                     >  2:401437/1‑143 (MQ=255)
  gTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCATATATGAACGGACTGGAATATATTGGCTAATACCCACTTTTTCTGAATCTAAAAAGGCCCTGCgcg                                                                                                                                   >  2:136209/1‑141 (MQ=255)
   tCGGACGCTGGGTGCAGGAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgc                                                                                                                                  <  1:59060/143‑1 (MQ=255)
   gcGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgc                                                                                                                                  <  1:339280/142‑1 (MQ=255)
       aCGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTAACAAAAGTTCTGCGTCCGTATATGAACGGACAGGAACAAAATGGCTAATAACCAAATTTTGCGAATCTAAAGAGCGCATGCGAGCGCtt                                                                                                                              >  2:252312/1‑143 (MQ=255)
          cTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                                                           >  1:137249/1‑143 (MQ=255)
           tGGTTGCAGCGATCGAAAACGATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                                                          <  1:475937/143‑2 (MQ=255)
            ggTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                         <  1:140244/143‑1 (MQ=255)
             ggtGCAGTAATGGAAAACTATCAGCAGGCTGATAGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                                                        <  1:399667/141‑4 (MQ=255)
                  aGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                                                   >  1:355147/1‑143 (MQ=255)
                        ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCAGGCGCAATTTTTGTCGCCCTTTGa                                                                                                             >  2:54350/1‑143 (MQ=255)
                           aaaCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAGAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAc                                                                                                          <  2:103569/143‑1 (MQ=255)
                             aCTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCg                                                                                                        <  2:504938/143‑1 (MQ=255)
                                 tCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGATCa                                                                                                    >  1:489762/1‑143 (MQ=255)
                                 tCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGAGCCTACGGGAGATTTTTTTGTCGCCCTTTCATACCAAACa                                                                                                    >  2:401314/1‑143 (MQ=255)
                                    gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTCAAAGGCGCCTGCGGGCGCGTTTCTTGTCTCCCTTTGGTTCCGTAcccca                                                                                                 >  1:168339/1‑139 (MQ=255)
                                      agGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATCTTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                                                               >  1:75631/1‑143 (MQ=255)
                                           gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                          <  1:401437/143‑1 (MQ=255)
                                           caTGGTCGTATTGCAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                          <  1:186537/142‑1 (MQ=255)
                                               gttggaTTGACGTACCAGAATTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGTGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCa                                                                                      <  2:415627/138‑1 (MQ=255)
                                                tCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCAc                                                                                     >  1:346195/1‑143 (MQ=255)
                                                        aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAAGAAGTCCCCCTCACTTACACGt                                                                             >  2:81800/1‑143 (MQ=255)
                                                        aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACAGAACAATAATTAGTCCTCAATCacacct                                                                             >  2:534367/1‑141 (MQ=255)
                                                                      cggCGTCCGTATATGAACGGACTGGAATATATTGGCCAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                               <  1:125648/141‑1 (MQ=255)
                                                                      cTGCGTCCGTACACGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                               <  1:219495/143‑1 (MQ=255)
                                                                        gCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACCCGTAATACTCCTTTCGAGt                                                             >  2:394706/1‑143 (MQ=255)
                                                                        gCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTCCTTTCGAgc                                                             >  2:493913/1‑142 (MQ=255)
                                                                         cGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTg                                                            <  1:52235/143‑1 (MQ=255)
                                                                                atatGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGCGTGACAATCt                                                     >  1:233061/1‑143 (MQ=255)
                                                                                  atGAACGGACTGGAATATATTGGCTAATGCCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                   <  1:403648/143‑1 (MQ=255)
                                                                                                   tatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTc                                  >  1:158548/1‑143 (MQ=255)
                                                                                                   tatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTc                                  >  1:436963/1‑143 (MQ=255)
                                                                                                     ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaa                                <  1:200090/143‑1 (MQ=255)
                                                                                                      tGGCTAATACCCAATTTTTCTGAATCTAAGATGCGCGCTCTGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAACTACACCCCACTTACACGAAATACTACGATCAACCGAAAATCTACCTATCTCTTTGATTTTCaac                               >  2:408603/1‑142 (MQ=255)
                                                                                                            aTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTACTCACTTACACGTCATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCCAAttcttc                         >  2:75156/1‑143 (MQ=255)
                                                                                                            aTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGCGCGCTTTTTTTGTCTCACTTTGATCTCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGCGATAATCTACCTATCTCTTTGATTTTCAAATTATTc                         <  2:137249/143‑1 (MQ=255)
                                                                                                            aTACCCAATTTTTCTGAATCTAAAAAGCGCCGGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAACACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                         <  1:271117/143‑1 (MQ=255)
                                                                                                                  aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCCAATTATTCGATGta                   >  1:177640/1‑143 (MQ=255)
                                                                                                                  aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGta                   >  2:520443/1‑143 (MQ=255)
                                                                                                                       ttCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGACACCGAACAATAATTACTCCTCACTTCCCCGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACaa              <  1:164261/143‑1 (MQ=255)
                                                                                                                         cTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGc            <  2:98788/143‑1 (MQ=255)
                                                                                                                          tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTCCCTATCTCTTTGATTTTCCAATTATTCGATGTATACAAGct           >  1:434332/1‑142 (MQ=255)
                                                                                                                          tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGcc           >  1:132350/1‑143 (MQ=255)
                                                                                                                                gaaaaaGAGCCTGCGGGCGCTTTTTTTGTCGCCCTTTGATACCGAACAATAATTACTCCTCCCTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatata     <  2:140634/142‑1 (MQ=255)
                                                                                                                                aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGCTTTTCAAATTATTCGATGTATACAAGCCtatata     >  2:379921/2‑143 (MQ=255)
                                                                                                                                aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtatata     >  2:326455/2‑143 (MQ=255)
                                                                                                                                 aaaaaGCGCCTGCGGGCGTTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAg    <  2:269916/143‑1 (MQ=255)
                                                                                                                                  aaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGc   >  1:160731/1‑143 (MQ=255)
                                                                                                                                   aaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg  >  1:296725/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                   
TGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCG  >  NC_000913/940612‑940885

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: