Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,581,621 G→A A476A (GCC→GCT hsdM ← DNA methyltransferase M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,581,6210GA100.0% 86.9 / NA 28A476A (GCC→GCThsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/13);  total (15/13)

CATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCG  >  NC_000913/4581515‑4581753
                                                                                                          |                                                                                                                                    
cATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGga                                                                                                      <  1:751070/139‑1 (MQ=255)
                ttCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGt                                                                                                                         <  1:660964/104‑1 (MQ=255)
                ttCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGt                                                                                                                         >  2:660964/1‑104 (MQ=255)
                  cACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGTAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGa                                                                                    >  2:746759/1‑139 (MQ=255)
                          aGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGTAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                            <  1:746759/139‑1 (MQ=255)
                                      tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                               >  2:820106/1‑124 (MQ=255)
                                      tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                               <  1:820106/124‑1 (MQ=255)
                                       cAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTc                                                               <  2:560372/139‑1 (MQ=255)
                                              cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATc                                                                                       <  2:948771/108‑1 (MQ=255)
                                              cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATc                                                                                       >  1:948771/1‑108 (MQ=255)
                                              cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTg                                                        >  1:626331/1‑139 (MQ=255)
                                                            cGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCg                                          >  1:836014/1‑139 (MQ=255)
                                                              cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGGTGAACTTGCGCCAGCGGc                                        <  2:836014/139‑1 (MQ=255)
                                                                    gCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAg                                  >  2:243327/1‑139 (MQ=255)
                                                                      ttCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCa                                <  1:164352/139‑1 (MQ=255)
                                                                         tGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGa                             >  1:1108867/1‑139 (MQ=255)
                                                                             gCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                             <  2:1063310/39‑1 (MQ=38)
                                                                             gCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                             >  1:1063310/1‑39 (MQ=255)
                                                                                     aTCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGgtg                 >  2:297159/1‑139 (MQ=255)
                                                                                        cggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTg                                                  >  2:638620/1‑103 (MQ=255)
                                                                                        cggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTg                                                  <  1:638620/103‑1 (MQ=255)
                                                                                           ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCa                                                                                              >  2:1146128/1‑56 (MQ=255)
                                                                                           ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCa                                                                                              <  1:1146128/56‑1 (MQ=255)
                                                                                           ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTg           >  2:1288451/1‑139 (MQ=255)
                                                                                           ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTg           >  1:744603/1‑139 (MQ=255)
                                                                                              cggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                          <  2:421031/41‑1 (MQ=255)
                                                                                              cggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                          >  1:421031/1‑41 (MQ=255)
                                                                                                    gCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCg  <  2:787358/139‑1 (MQ=255)
                                                                                                          |                                                                                                                                    
CATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCG  >  NC_000913/4581515‑4581753

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