Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 345,649 G→A *82* (TAG→TAA yahM → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,6490GA100.0% 99.3 / NA 32*82* (TAG→TAAyahMuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/17);  total (15/17)

AAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACA  >  NC_000913/345516‑345774
                                                                                                                                     |                                                                                                                             
aaGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                          <  1:261849/139‑1 (MQ=255)
  gCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTaa                                                                                                                        <  1:151048/139‑1 (MQ=255)
                 agTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                         >  2:114701/1‑139 (MQ=255)
                      gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCg                                                                                                    >  1:696323/1‑139 (MQ=255)
                                    cTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTaa                                                                                      >  1:659241/1‑139 (MQ=255)
                                       aaGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATc                                                                                   <  2:659241/139‑1 (MQ=255)
                                                   aaCTTATCGGAAATAAGGATCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                    >  2:105647/1‑126 (MQ=255)
                                                   aaCTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                    <  1:105647/126‑1 (MQ=255)
                                                         tCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTaaag                                                                                                                             >  1:267033/1‑79 (MQ=255)
                                                                  aGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                     >  2:292608/1‑126 (MQ=255)
                                                                  aGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                     <  1:292608/126‑1 (MQ=255)
                                                                        aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCa                                                  >  2:265738/1‑139 (MQ=255)
                                                                                            tcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAAcgcg                              <  1:114701/139‑1 (MQ=255)
                                                                                            tcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAAcgcg                              <  1:265738/139‑1 (MQ=255)
                                                                                               ccGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCtt                                                                                           <  1:56676/75‑1 (MQ=255)
                                                                                               ccGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCtt                                                                                           >  2:56676/1‑75 (MQ=255)
                                                                                               ccGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                    >  2:828288/1‑82 (MQ=255)
                                                                                               ccGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                    <  1:828288/82‑1 (MQ=255)
                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTATGTGAAGAGTAAATTATCTGCTCACTGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                       <  1:514608/107‑1 (MQ=255)
                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                       >  2:514608/1‑107 (MQ=255)
                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCa                                   >  2:228443/1‑127 (MQ=255)
                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCa                                   <  1:228443/127‑1 (MQ=255)
                                                                                                     gCCGTGATTCGGCAAATCAGGATGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCTCTTCGCAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACg                     <  2:432468/139‑1 (MQ=255)
                                                                                                         tGATTCGGCATATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                         >  1:370244/1‑51 (MQ=255)
                                                                                                         tGATTCGGCATATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                         <  2:370244/51‑1 (MQ=255)
                                                                                                               ggCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCAcc                                                   >  2:372309/1‑99 (MQ=255)
                                                                                                               ggCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCAcc                                                   <  1:372309/99‑1 (MQ=255)
                                                                                                                    aTCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCacgacg      >  1:244657/1‑139 (MQ=255)
                                                                                                                        ggTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGcaca  <  2:244657/139‑1 (MQ=255)
                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACg                     <  1:3316/115‑1 (MQ=255)
                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACg                     >  2:3316/1‑115 (MQ=255)
                                                                                                                                 cGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGc                                    >  2:153288/1‑96 (MQ=255)
                                                                                                                                 cGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGc                                    <  1:153288/96‑1 (MQ=255)
                                                                                                                                     |                                                                                                                             
AAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACA  >  NC_000913/345516‑345774

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