Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 276,946 T→G 7.1% F78V (TTC→GTC)  mmuM → CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913276,9460TG 7.1% 199.9 / 6.0 56F78V (TTC→GTC) mmuMCP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base G (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCG  >  NC_000913/276807‑277093
                                                                                                                                           |                                                                                                                                                   
gtgtAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcg                                                                                                                                             <  1:709077/148‑1 (MQ=255)
   tAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa                                                                                                                                          <  4:998581/148‑1 (MQ=255)
      cTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAg                                                                                                                                         >  1:1029835/1‑146 (MQ=255)
      cTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAg                                                                                                                                         <  2:1029835/146‑1 (MQ=255)
         aGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTg                                                                                                                                    <  2:301041/148‑1 (MQ=255)
                    tgtgGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTg                                                                                                                                    <  3:787123/137‑1 (MQ=255)
                    ggtgGTCAGCCAAAGTGCGGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATGTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTg                                                                                                                                    >  4:787123/2‑137 (MQ=255)
                              cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcg                                                                                                                                                  <  1:1031643/113‑1 (MQ=255)
                              cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcg                                                                                                                                                  >  2:1031643/1‑113 (MQ=255)
                              cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa                                                                                                                                          <  2:821025/121‑1 (MQ=255)
                              cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa                                                                                                                                          >  1:821025/1‑121 (MQ=255)
                              cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcgc                                                                                                                          <  1:470104/137‑1 (MQ=255)
                              cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcgc                                                                                                                          >  2:470104/1‑137 (MQ=255)
                                       ggTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGc                                                                                                      >  1:23232/1‑148 (MQ=255)
                                        gTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCa                                                                                                     >  3:446218/1‑148 (MQ=255)
                                         tAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGa                                                                                                                                  >  2:1053379/1‑118 (MQ=255)
                                         tAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGa                                                                                                                                  <  1:1053379/118‑1 (MQ=255)
                                          aGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCaaa                                                                                                   >  4:489693/1‑148 (MQ=255)
                                             aaaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAAGCCCGCCGGCGTcgcctgata                                                                                                                                           <  1:298112/105‑6 (MQ=38)
                                             aaaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAAGCCCGCCGGCGTcgcctgata                                                                                                                                           >  2:298112/1‑100 (MQ=38)
                                              aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGa                                                                                                                               <  1:337483/116‑1 (MQ=255)
                                              aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGa                                                                                                                               >  2:337483/1‑116 (MQ=255)
                                              aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAAGCCCGCCGGCGTcgcct                                                                                                                                               >  3:700704/1‑99 (MQ=255)
                                              aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCAATCAGGCGAAGCCCGCCGGCGTcgcct                                                                                                                                               <  4:700704/100‑2 (MQ=255)
                                               aCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcgcg                                                                                                                                           <  3:756988/103‑1 (MQ=255)
                                               aCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcgcg                                                                                                                                           >  4:756988/1‑103 (MQ=255)
                                                 ccGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTgg                                                                                            <  2:23232/148‑1 (MQ=255)
                                                           tCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGcgcg                                                                                  >  1:600714/1‑148 (MQ=255)
                                                               cgAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGaa                                                                                                                   >  4:582228/1‑111 (MQ=255)
                                                               cgAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGaa                                                                                                                   <  3:582228/111‑1 (MQ=255)
                                                                        tCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa                                                                                                                                          >  2:670438/1‑79 (MQ=255)
                                                                        tCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa                                                                                                                                          <  1:670438/79‑1 (MQ=255)
                                                                        tCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGaa                                                                     >  3:995996/1‑148 (MQ=255)
                                                                           tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTgg                                                                                                                                   <  2:930843/83‑1 (MQ=255)
                                                                           tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTgg                                                                                                                                   >  1:930843/1‑83 (MQ=255)
                                                                           tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCTTTCGactgcg                                                                                                               >  1:591986/1‑98 (MQ=255)
                                                                           tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCTTTCGactgcg                                                                                                               <  2:591986/103‑6 (MQ=255)
                                                                             aTTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg                                                                                               <  3:719117/117‑1 (MQ=255)
                                                                             aTTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg                                                                                               >  4:719117/1‑117 (MQ=255)
                                                                              ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGAtt                                                                                                         <  4:75104/106‑1 (MQ=255)
                                                                              ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGAtt                                                                                                         >  3:75104/1‑106 (MQ=255)
                                                                              ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt                                                                                              >  3:768182/1‑117 (MQ=255)
                                                                              ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt                                                                                              <  4:768182/117‑1 (MQ=255)
                                                                              ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTAt                                                               <  4:995996/148‑1 (MQ=255)
                                                                               tactacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGa                                                                      >  2:664550/1‑140 (MQ=255)
                                                                               tactacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGa                                                                      <  1:664550/140‑1 (MQ=255)
                                                                                actacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg                                                                                               <  4:27400/114‑1 (MQ=255)
                                                                                actacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg                                                                                               >  3:27400/1‑114 (MQ=255)
                                                                                   acCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTaa                                                                               <  1:369079/127‑1 (MQ=255)
                                                                                   acCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTaa                                                                               >  2:369079/1‑127 (MQ=255)
                                                                                          gggggCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAAc                                                   <  2:600714/148‑1 (MQ=255)
                                                                                           ggggCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAAcc                                                  <  3:489693/148‑1 (MQ=255)
                                                                                                   aTGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCg                                          <  4:446218/148‑1 (MQ=255)
                                                                                                              cTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt                                                                                              <  1:403942/85‑1 (MQ=255)
                                                                                                              cTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt                                                                                              >  2:403942/1‑85 (MQ=255)
                                                                                                                                          cTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGgc   >  3:1030137/1‑148 (MQ=255)
                                                                                                                                           ttCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGgcg  <  4:1030137/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCGGGAACGCTTCTGGTGGCGGGATCCGTCGGGCCTTACGGCG  >  NC_000913/276807‑277093

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑