Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 441,092 G→A P4S (CCG→TCG)  ispA ← geranyl diphosphate/farnesyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913441,0920GA100.0% 61.3 / NA 22P4S (CCG→TCG) ispAgeranyl diphosphate/farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/13);  total (9/13)

ATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCGAACTCGTTCAGCGCCTCTTCCAGCGGCAGGTCGCCACTTTCCAGACGGGTTACAATCTGTTCCAGCTCGCT  >  NC_000913/440857‑441298
                                                                                                                                                                                                                                           |                                                                                                                                                                                                              
aTTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTTGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGGGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCa                                                                                                                                                                                                 <  1:290702/251‑1 (MQ=255)
                                                     gtTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAg                                                                                                                                                                                                            >  1:178646/1‑187 (MQ=255)
                                                            cTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCGTGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGCGTTGCTGCGAAACGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcg                                                                                                                                     >  1:145542/1‑251 (MQ=255)
                                                                  ccGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGTGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCAtt                                                                                                                                                                                                      <  1:241935/180‑1 (MQ=255)
                                                                                                          nnnnnntANNNNNTAANAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTc                                                                                                                                                                                            <  2:314689/144‑1 (MQ=255)
                                                                                                                  ccaccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCGTCGAGTGGCTGCGAAAAGTCCATGACTCATTGTCCGGTTTAAAATGGGGTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGcc                                                                                                          <  1:492615/224‑1 (MQ=255)
                                                                                                                           aatGCGCCATACTGCATGGTTTCGACCCCGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAAAGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGGAAAAGGGGTTAGAGAGGCGTCTTCATTGACAGACAGCAGCATTTGTACGCGCTGTTCGACTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCGAAc                                                                      >  1:32576/1‑251 (MQ=255)
                                                                                                                                   aTACTGCATGGTTTCGACCACGGGAGTGTGCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGGTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAc                                                                                                                                      <  1:255252/179‑1 (MQ=255)
                                                                                                                                                 cGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGc                                                                                     <  1:75970/214‑1 (MQ=255)
                                                                                                                                                    ccACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGACTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTAATCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCCTTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTg                                                                                                        >  1:207458/1‑192 (MQ=255)
                                                                                                                                                                   gAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCGAACGCGTTCAGCGCCTCTTCCAGCGGCAGGTCGCCACTTTcc                               <  1:368806/250‑1 (MQ=255)
                                                                                                                                                                   gAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTGGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGTTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCGAACTCGTTCAGCGCCTCTTCCAGCGGCAGGTCGCCACTTTcc                               <  1:376361/250‑1 (MQ=255)
                                                                                                                                                                         gCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTaga                                                                                                                                                                         >  1:389952/1‑106 (MQ=255)
                                                                                                                                                                             tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGttn                                                                                                                             <  1:92821/146‑2 (MQ=255)
                                                                                                                                                                                 gCGATTAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGCGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTACAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCACTGACGTGCCAGCTGCACGCCGCGTTCGAACTCGTTCAGCGCCCCTTCCAGCGGCAggacaccac                                    >  1:106397/1‑224 (MQ=255)
                                                                                                                                                                                 gCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCCGCTGCACGCCGCGTTCGAACTCGTTCAGCGCCTCTTCCAGCGGCAGGTCGCCAc                                    >  1:106402/1‑231 (MQ=255)
                                                                                                                                                                                       aaaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACGCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCAttn                                                                                                                                                          <  1:56318/107‑2 (MQ=255)
                                                                                                                                                                                               cAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcgc                                                                                                                                    <  1:70351/121‑1 (MQ=255)
                                                                                                                                                                                               cAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCAGGTGAAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCCGAATTTGTACGAGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCGAACTCGTTCAGCGCCTCTTCCAGCGGCAGGTCGCCACTTTCCAGACGGGTTACAATCTGTTCCAgctcgct  >  1:325400/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTcc                                                                                                                                                                                           >  1:30489/1‑38 (MQ=255)
                                                                                                                                                                                                                           gCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTcc                                                                                                                                                                                           <  2:30489/38‑1 (MQ=255)
                                                                                                                                                                                                                           gCTTCGAGTTGCTGCGAAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCTAACTCGTTCAGCGCCTCTGGCAGCGGCAGGTCGc                                       <  1:179749/186‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                              
ATTAATGAGTAAGCGTGGATACACTCAACGGCGGCAGCGGGTGCGTCCAGCGTGTTTGTGCTAACGCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGCACGCCGCGTTCGAACTCGTTCAGCGCCTCTTCCAGCGGCAGGTCGCCACTTTCCAGACGGGTTACAATCTGTTCCAGCTCGCT  >  NC_000913/440857‑441298

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