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% 134.1 / NA 42*82* (TAG→TAAyahMuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

CACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCAAA  >  NC_000913/345523‑345783
                                                                                                                              |                                                                                                                                      
cacaGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                            <  1:542772/139‑1 (MQ=255)
cacaGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                            <  2:360544/139‑1 (MQ=255)
          agTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                                  <  2:148711/139‑1 (MQ=255)
          agTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                                  >  1:213191/1‑139 (MQ=255)
               gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCg                                                                                                             >  2:437273/1‑139 (MQ=255)
               gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCg                                                                                                             >  1:232220/1‑139 (MQ=255)
               gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCg                                                                                                             >  2:65809/1‑139 (MQ=255)
                cTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGt                                                                                                            <  2:24392/139‑1 (MQ=255)
                               aaaGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             >  2:625752/1‑139 (MQ=255)
                                aaGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                            <  2:518481/107‑1 (MQ=255)
                                aaGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                            >  1:518481/1‑107 (MQ=255)
                                      cAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                                      <  2:88324/139‑1 (MQ=255)
                                      cAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                                      <  2:213191/139‑1 (MQ=255)
                                               ttATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                             <  1:188055/139‑1 (MQ=255)
                                                tATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATa                                                                            <  1:731021/139‑1 (MQ=255)
                                                 aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTg                                                                                                               >  1:428688/1‑103 (MQ=255)
                                                 aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTg                                                                                                               <  2:428688/103‑1 (MQ=255)
                                                 aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCg                                                                                                  <  1:155772/116‑1 (MQ=255)
                                                 aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCg                                                                                                  >  2:155772/1‑116 (MQ=255)
                                                   cGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATacac                                                                         >  1:111248/1‑139 (MQ=255)
                                                    ggAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACAcc                                                                        <  1:146811/139‑1 (MQ=255)
                                                                cACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCAcc                                                            >  1:41170/1‑139 (MQ=255)
                                                                 aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCCCCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                              >  2:176640/1‑120 (MQ=255)
                                                                 aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             >  1:84103/1‑105 (MQ=255)
                                                                 aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             <  2:84103/105‑1 (MQ=255)
                                                                 aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                              <  1:176640/120‑1 (MQ=255)
                                                                 aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGACAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCa                                                           >  2:252110/1‑139 (MQ=255)
                                                                       tACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATa                                                                            <  2:199714/116‑1 (MQ=255)
                                                                       tACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATa                                                                            >  1:199714/1‑116 (MQ=255)
                                                                             ggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             <  1:537886/93‑1 (MQ=255)
                                                                             ggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             >  2:537886/1‑93 (MQ=255)
                                                                              gCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             <  1:63416/92‑1 (MQ=255)
                                                                              gCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             >  2:63416/1‑92 (MQ=255)
                                                                                 atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                                >  1:558108/1‑118 (MQ=255)
                                                                                 atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                                <  2:558108/118‑1 (MQ=255)
                                                                                    atcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAAcgc                                        <  2:138436/139‑1 (MQ=255)
                                                                                        ccGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTg                                    >  1:352059/1‑139 (MQ=255)
                                                                                            tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                                >  1:739982/1‑139 (MQ=255)
                                                                                              gCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                                      <  1:413259/83‑1 (MQ=255)
                                                                                              gCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                                      >  2:413259/1‑83 (MQ=255)
                                                                                                                 ggTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGcaca           <  1:489560/139‑1 (MQ=255)
                                                                                                                          cGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCaaa  >  2:117616/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                      
CACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCAAA  >  NC_000913/345523‑345783

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