Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,340 T→C 100% T65A (ACT→GCT)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,3400TC98.8% 271.5 / ‑5.2 81T65A (ACT→GCT) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (38/42);  minor base G (1/0);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATAT  >  NC_000913/1197244‑1197435
                                                                                                 |                                                                                                
cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCACTCTTACAAGTTGCATTTTCGAGTCTAGAAAAGCAAc                                                                                              >  2:1726760/1‑100 (MQ=255)
      tataATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCTTCTAATAATTTGCTCTCTTTCAAGTTAACTTTTCGAGTCTAGAAAAGCAACTtcatc                                                                                        >  2:455897/1‑100 (MQ=255)
       ataataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcatca                                                                                       >  1:1184862/1‑100 (MQ=255)
       ataataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcatca                                                                                       >  2:2493445/1‑100 (MQ=255)
        taataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCaa                                                                                      >  1:1438911/1‑100 (MQ=255)
        taataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCaa                                                                                      >  2:313142/1‑100 (MQ=255)
          atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAt                                                                                    >  2:2554303/1‑100 (MQ=255)
          atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAt                                                                                    >  1:434286/1‑100 (MQ=255)
          atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAt                                                                                    <  2:2236470/100‑1 (MQ=255)
          atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAt                                                                                    <  1:1510534/100‑1 (MQ=255)
           tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAtt                                                                                   <  2:1735211/100‑1 (MQ=255)
              cGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTAc                                                                                <  2:697438/100‑1 (MQ=255)
               gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTAc                                                                                <  2:412238/99‑1 (MQ=255)
               gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               <  1:973704/100‑1 (MQ=255)
                tcatacacAGTCAATACTTTTACTTTTAATTCCCCTAATATTGTACTCTCTTTCAAGTTACTTTTTCGAGTCTACAAAAGCACCTTTCACTAATTTAc                                                                                <  2:2308780/94‑1 (MQ=255)
                tGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTg                                                                              <  2:290862/100‑1 (MQ=255)
                   gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                           <  1:1124406/100‑1 (MQ=255)
                   gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                           >  2:1882905/1‑100 (MQ=255)
                   gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTATCTTTCAAGTTGCCTTTTCGAGTCTAGAACAGCAACTTCATCAAATTTACTGaat                                                                           >  2:843296/1‑98 (MQ=255)
                    aCATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCCCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                               <  2:230311/95‑1 (MQ=255)
                       tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGGAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                       <  1:1726772/100‑1 (MQ=255)
                       tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                       >  2:1731449/1‑100 (MQ=255)
                       tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                       >  1:2144996/1‑100 (MQ=255)
                        aGTCAATAATTTCAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtc                                                                      >  1:1686957/1‑100 (MQ=255)
                         gTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtca                                                                     >  1:1121218/1‑100 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTGTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   >  1:1053497/1‑100 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   >  1:309515/1‑100 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   <  1:356666/100‑1 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTCACTGATTAATtcatc                                                                   <  2:2304772/100‑1 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCAAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                   >  2:1645112/1‑100 (MQ=255)
                            aataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatca                                                                  <  1:1284961/100‑1 (MQ=255)
                             ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCATCTTCATCAAATTTACTGATTAATTCATCAg                                                                 <  2:1686945/100‑1 (MQ=255)
                               natTTTAGTTTTAATTCCTCTAATAATTTGCTATCTTTCAAGTTGCATTTTCTATTATATAAAAGCAACTTAATCAAATTTACTGATTAATTAATCAGaa                                                               >  2:2467665/2‑100 (MQ=255)
                                 ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                             <  1:758565/100‑1 (MQ=255)
                                 ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                             <  2:2165216/100‑1 (MQ=255)
                                   ttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATc                                                           <  1:2290222/100‑1 (MQ=255)
                                        tttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGt                                                      <  1:1722857/100‑1 (MQ=255)
                                              tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                                <  1:2065403/100‑1 (MQ=255)
                                               ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                               <  1:233481/100‑1 (MQ=255)
                                                ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                              >  1:635415/1‑100 (MQ=255)
                                                  ctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGc                                            >  2:2288046/1‑100 (MQ=255)
                                                   taataaTTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                           <  2:1762785/100‑1 (MQ=255)
                                                    aataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCa                                          >  1:928468/1‑100 (MQ=255)
                                                    aataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCa                                          >  1:792293/1‑100 (MQ=255)
                                                    aataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCa                                          >  2:61736/1‑100 (MQ=255)
                                                    aataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCa                                          >  2:780641/1‑100 (MQ=255)
                                                        atTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAAt                                      <  1:2319026/100‑1 (MQ=255)
                                                          ttGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCa                                    <  1:2267483/100‑1 (MQ=255)
                                                           tGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCaa                                   <  1:291816/100‑1 (MQ=255)
                                                             ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                                 <  1:608630/100‑1 (MQ=255)
                                                             ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCAGCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                                 <  2:1437091/100‑1 (MQ=255)
                                                              tctctTTCAAGTTGCCTTTTTTATTTTTTTAAAGCAACTTCATCAAATTTACTGATTAATTTATCAGAATATCTAAGTTTTCCGCTGCCAAAATc                                     >  2:1524809/1‑95 (MQ=255)
                                                              tctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAggg                                <  2:1268169/100‑1 (MQ=255)
                                                              tctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAggg                                >  2:2432971/1‑100 (MQ=255)
                                                              tctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgcg                                <  1:1731461/100‑3 (MQ=255)
                                                              tctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTGCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAggg                                <  2:2399925/100‑1 (MQ=255)
                                                                  tttCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTg                            <  2:366325/100‑1 (MQ=255)
                                                                  gatgaaGTTGCTTTTCTAGACTCGAAAAGGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTACGCAGCCAAAATCAAGGGCg                              >  2:2524248/5‑98 (MQ=21)
                                                                   ttCACGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATCCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                           <  1:756709/100‑1 (MQ=255)
                                                                   ttCAAGTTGCCTTTTCGAGTCTAGCAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                                 >  2:1127097/1‑94 (MQ=255)
                                                                   ttCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                           <  2:2141418/100‑1 (MQ=255)
                                                                   ttCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                           <  1:61737/100‑1 (MQ=255)
                                                                   ttCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                           <  1:1529677/100‑1 (MQ=255)
                                                                      aaGTTGCTTTTCTAGACTCGAAAAGGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCa                        >  1:438681/1‑100 (MQ=37)
                                                                           gCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAttttt                   >  1:1624774/1‑100 (MQ=255)
                                                                           gCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAttttt                   >  1:2104027/1‑100 (MQ=255)
                                                                            ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                  >  1:2032703/1‑100 (MQ=255)
                                                                             cTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCTGCAGCCAAAATCAAGGGCGTGGCCATTTTtct                 <  1:2427304/100‑1 (MQ=255)
                                                                              ttttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc                <  2:2287753/100‑1 (MQ=255)
                                                                                  cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg            >  2:162411/1‑100 (MQ=255)
                                                                                   gAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGc           <  2:2352224/100‑1 (MQ=255)
                                                                                      tCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTc        >  2:452523/1‑100 (MQ=255)
                                                                                        tAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTa      >  1:2066452/1‑100 (MQ=255)
                                                                                        tAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTa      >  1:2040096/1‑100 (MQ=255)
                                                                                        tAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATAAAGGGCGTGGCCATTTTTCTCAATGCTTCTa      >  2:647246/1‑100 (MQ=255)
                                                                                        tAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGCTTTCCGCAGCCAAAATCAAGGACGTGGCCATTTTTCTCAAtcgtgata      >  2:1582271/1‑93 (MQ=255)
                                                                                         agcaaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGTAGCCAAGATCAAGGGCGTGGCCAGTTTTCTCAATGGTTCTaa     <  2:464121/97‑1 (MQ=255)
                                                                                          gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat    <  2:406727/100‑1 (MQ=255)
                                                                                           aaaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata   <  1:228766/100‑1 (MQ=255)
                                                                                           aaaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata   <  2:1652084/100‑1 (MQ=255)
                                                                                            aaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat  >  2:873929/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATAT  >  NC_000913/1197244‑1197435

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