Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,874,741 T→C T227A (ACA→GCA)  cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,874,7410TC100.0% 93.1 / NA 25T227A (ACA→GCA) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/10);  total (15/10)

CCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCA  >  NC_000913/2874611‑2874857
                                                                                                                                  |                                                                                                                    
ccGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGc                                                                                                              >  1:70487/1‑139 (MQ=255)
               ccGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATAc                                                                                               <  1:52901/139‑1 (MQ=255)
                    gTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCa                                                                                          <  2:284798/139‑1 (MQ=255)
                                      aTTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTg                                                                        <  2:37008/139‑1 (MQ=255)
                                           ggCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGa                                                                                                                  >  1:396678/1‑92 (MQ=255)
                                           ggCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGa                                                                                                                  <  2:396678/92‑1 (MQ=255)
                                               aTTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCg                                                                     <  1:509277/133‑1 (MQ=255)
                                               aTTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCg                                                                     >  2:509277/1‑133 (MQ=255)
                                                         tGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGa                                                                                                                  >  2:597824/1‑78 (MQ=255)
                                                         tGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGa                                                                                                                  <  1:597824/78‑1 (MQ=255)
                                                          gCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCagagag                                                    >  2:347455/1‑139 (MQ=255)
                                                           cACCTGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGCGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTATCGCCTTCCAGTGCAGAGAGc                                                   >  1:358229/1‑139 (MQ=255)
                                                              cGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGc                                                >  2:283699/1‑139 (MQ=255)
                                                                ggAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGcac                                              >  2:464276/1‑139 (MQ=255)
                                                                            gCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATc                                  <  1:464276/139‑1 (MQ=255)
                                                                                          acTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                    >  2:499581/1‑139 (MQ=255)
                                                                                          acTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCa                    >  2:25800/1‑139 (MQ=255)
                                                                                            tCGCTGGATCTCCACGGTTTCCAGCACTTCGAGTAGTGCCGGACCGCTGTACCACGGCATAATTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGc                  >  1:245043/1‑139 (MQ=255)
                                                                                              gCTGGATCTCCACGGTTTCCAGCACTTCGCGCAGTGCCGGACCGCAGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTcc                                                               >  1:355590/1‑92 (MQ=255)
                                                                                              gCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTcc                                                               <  2:355590/92‑1 (MQ=255)
                                                                                              gCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTcc                                                               <  2:453521/92‑1 (MQ=255)
                                                                                              gCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTcc                                                               >  1:453521/1‑92 (MQ=255)
                                                                                                tGGATCTCCACGGTTTCCAGCACTTTGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGcc              >  2:308594/1‑139 (MQ=255)
                                                                                                          cGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGAt                             <  1:80809/114‑1 (MQ=255)
                                                                                                          cGGTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGAt                             >  2:80809/1‑114 (MQ=255)
                                                                                                            gTTTCCAGCACTTCGAGCAGTGCCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCa  >  2:239936/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                    
CCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTTCGAGCAGTGTCGGACCGCTGTACCACGGCATACTTTCACTTTGCGATGCCACGTTGTCGCCTTCCAGTGCAGAGAGCGGCACAAAGCGGATATCCAGATTACCCGGCAGCTGCCCGGCAAAGGTCA  >  NC_000913/2874611‑2874857

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: