Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,340 T→C 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,3400TC100.0% 168.4 / NA 47T65A (ACT→GCT) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/17);  total (30/17)

CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAAT  >  NC_000913/1197244‑1197433
                                                                                                |                                                                                             
cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAAc                                                                                            >  2:508228/1‑100 (MQ=255)
   ttttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACt                                                                                           >  2:956551/1‑98 (MQ=255)
   ttttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtc                                                                                         >  1:825058/1‑100 (MQ=255)
          atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAt                                                                                  >  1:510920/1‑100 (MQ=255)
           tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAtt                                                                                 >  1:313873/1‑100 (MQ=255)
               gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                             >  2:1323894/1‑100 (MQ=255)
                tGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTg                                                                            <  1:193677/100‑1 (MQ=255)
                  ggACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAt                                                                          >  2:182154/1‑100 (MQ=255)
                   gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCTAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                         >  2:311675/1‑100 (MQ=255)
                       tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCACTTTCAAGTTGCCTTTTCTAGTCTAGAAAATCAACTTCATCAAATTTACTGATTAAtt                                                                     >  2:861675/1‑100 (MQ=255)
                        aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                     >  2:184539/1‑99 (MQ=255)
                         gTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtca                                                                   <  2:190257/100‑1 (MQ=255)
                           cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                 >  1:376840/1‑100 (MQ=255)
                             ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAg                                                               <  2:29451/100‑1 (MQ=255)
                               aatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGaa                                                             >  2:156008/1‑100 (MQ=255)
                               aatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGaa                                                             <  2:1169839/100‑1 (MQ=255)
                                 ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                           <  2:75874/100‑1 (MQ=255)
                                     aGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTa                                                       <  1:685250/100‑1 (MQ=255)
                                      gTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTCCTGATTAATTCATCAGataat                                                          >  2:1340911/1‑92 (MQ=255)
                                      gTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTaa                                                      >  1:58279/1‑100 (MQ=255)
                                       ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAg                                                     <  2:858934/100‑1 (MQ=255)
                                       ttttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAg                                                     <  1:1289078/100‑1 (MQ=255)
                                              tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                              >  1:486507/1‑100 (MQ=255)
                                              tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                              <  1:426437/100‑1 (MQ=255)
                                              tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                              <  2:1514333/100‑1 (MQ=255)
                                                ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                            >  1:39575/1‑100 (MQ=255)
                                                ctctAATAATTTGCTCTCTTTCAAGTTGCCTATTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                            >  1:39583/1‑100 (MQ=255)
                                                   taataaTTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                         <  1:1137624/100‑1 (MQ=255)
                                                      taatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCaaa                                      >  1:205902/1‑100 (MQ=255)
                                                      taatTTGCTCTATTTCAAGTTGCCTTTTCGAGTATAGAAAAGCAACTTAATCAAATTTACTGATTAATTCATCAGAATATCTACGTATTCCGCAGcta                                        >  2:674360/1‑96 (MQ=255)
                                                             ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                               >  2:165823/1‑100 (MQ=255)
                                                             ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                               >  2:227097/1‑100 (MQ=255)
                                                                 ctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                           >  1:478013/1‑100 (MQ=255)
                                                                 ctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                           >  1:567988/1‑100 (MQ=255)
                                                                       aGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAt                     >  2:1380474/1‑100 (MQ=255)
                                                                       aGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAt                     >  2:810791/1‑100 (MQ=255)
                                                                        gTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                    >  1:1046323/1‑100 (MQ=255)
                                                                            ccTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtc                <  2:514108/100‑1 (MQ=255)
                                                                             cTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtct               <  1:483294/100‑1 (MQ=255)
                                                                              ttttCGAGTCTAGAAAAGCAACCTCATCAAATTTACTCATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc              <  2:1065675/100‑1 (MQ=255)
                                                                                  cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg          <  2:255361/100‑1 (MQ=255)
                                                                                   gAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGc         >  2:1227949/1‑100 (MQ=255)
                                                                                    aGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt        <  2:62597/100‑1 (MQ=255)
                                                                                    aGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt        >  1:1163331/1‑100 (MQ=255)
                                                                                     gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCtt       >  2:801943/1‑100 (MQ=255)
                                                                                         aGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTaa   >  1:811477/1‑100 (MQ=255)
                                                                                          gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat  <  2:1229700/100‑1 (MQ=255)
                                                                                                |                                                                                             
CATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAAT  >  NC_000913/1197244‑1197433

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