Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,751 T→C 10.0% 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.0% 118.0 / 4.0 50intergenic (+31/‑208)serS/dmsAserine‑‑tRNA ligase/dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (26/18);  new base C (3/2);  total (29/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 = 6.17e-01

TGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCG  >  NC_000913/940609‑940885
                                                                                                                                              |                                                                                                                                      
tGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCt                                                                                                                                        >  1:51912/1‑143 (MQ=255)
  ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGc                                                                                                                                      >  1:136995/1‑143 (MQ=255)
   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)
                                                                                                                                              |                                                                                                                                      
TGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCG  >  NC_000913/940609‑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: