Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,436 C→T L47L (CTG→CTA thiF ← ThiS adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,4360CT100.0% 58.7 / NA 18L47L (CTG→CTAthiFThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/6);  total (12/6)

GTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCAC  >  NC_000913/4193306‑4193566
                                                                                                                                  |                                                                                                                                  
gTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAgc                                                                                                                           >  1:505414/1‑140 (MQ=255)
     tGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCgcagc                                                                                                                      >  2:44982/1‑140 (MQ=255)
        tGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAgg                                                                                                                   <  1:136025/140‑1 (MQ=255)
              aCCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTAcc                                                                                                             >  1:481899/1‑140 (MQ=255)
                         cGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCa                                                                                                  <  1:44982/140‑1 (MQ=255)
                              gatcgatATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCg                                                                                             <  2:505414/140‑1 (MQ=255)
                                           aGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCt                                                                                >  2:175030/1‑140 (MQ=255)
                                                  aaGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTc                                                                         >  1:412727/1‑140 (MQ=255)
                                                       gATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCa                                                                    >  1:27551/1‑140 (MQ=255)
                                                                  gCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTtgctgc                                                         <  2:481899/140‑1 (MQ=255)
                                                                                                 tGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATAtcg                                       >  1:267126/1‑127 (MQ=255)
                                                                                                 tGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATAtcg                                       <  2:267126/127‑1 (MQ=255)
                                                                                                       tACCAGCGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTggcggc                    >  2:296129/1‑140 (MQ=255)
                                                                                                             cGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAAc              >  2:558563/1‑140 (MQ=255)
                                                                                                             cGTCCCGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCACGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAAc              >  2:35478/1‑140 (MQ=255)
                                                                                                                  cGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATa         >  1:523990/1‑140 (MQ=255)
                                                                                                                  cGACGCCAGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATa         <  1:296129/140‑1 (MQ=255)
                                                                                                                         aGCGCCCGCTAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCAc  >  1:169941/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                  
GTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCGACGCCAGCGCCCGCCAGGTACAGCGCAGCAGGTGTACCCAGCCCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCAC  >  NC_000913/4193306‑4193566

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