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

AAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACA  >  NC_000913/345516‑345774
                                                                                                                                     |                                                                                                                             
aaGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                          <  2:187587/139‑1 (MQ=255)
            gggAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCt                                                                                                              <  1:510438/139‑1 (MQ=255)
                 agTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                         >  2:1489161/1‑139 (MQ=255)
                               ggAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCtt                                                                                           >  2:604511/1‑139 (MQ=255)
                                  aGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGt                                                                                        <  1:1468847/139‑1 (MQ=255)
                                       aaGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATc                                                                                   <  2:1303490/139‑1 (MQ=255)
                                           cTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGc                                                                               <  2:1200836/139‑1 (MQ=255)
                                            tCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCt                                                                                                              <  2:1023073/107‑1 (MQ=255)
                                            tCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCt                                                                                                              >  1:1023073/1‑107 (MQ=255)
                                             cAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                             <  1:1489161/139‑1 (MQ=255)
                                             cAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                             <  1:1182648/139‑1 (MQ=255)
                                              aaGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTa                                                                            >  2:535504/1‑139 (MQ=255)
                                                     cTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                     <  1:535504/139‑1 (MQ=255)
                                                      ttATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                    <  2:1356957/139‑1 (MQ=255)
                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGt                                                                                                   <  2:537337/106‑1 (MQ=255)
                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGt                                                                                                   >  1:537337/1‑106 (MQ=255)
                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATac                                                                  >  1:1224479/1‑139 (MQ=255)
                                                              aaTAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCaat                                                            <  2:948624/139‑1 (MQ=255)
                                                                tAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                               >  2:433498/1‑102 (MQ=255)
                                                                tAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                               <  1:433498/102‑1 (MQ=255)
                                                                         cGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCaa                                                 >  2:628484/1‑139 (MQ=255)
                                                                         cGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAACTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                    >  1:1248747/1‑104 (MQ=255)
                                                                         cGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAACTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                    <  2:1248747/104‑1 (MQ=255)
                                                                                   gggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTaaag                                                                                                                             >  2:1433492/1‑53 (MQ=255)
                                                                                    ggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                   >  2:663651/1‑62 (MQ=255)
                                                                                    ggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                   <  1:663651/62‑1 (MQ=255)
                                                                                      cTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                          <  1:427511/53‑1 (MQ=255)
                                                                                      cTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                          >  2:427511/1‑53 (MQ=255)
                                                                                      cTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGc                                    <  2:1224479/139‑1 (MQ=255)
                                                                                        atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGc                                                                                                        <  1:1151850/69‑1 (MQ=255)
                                                                                        atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGc                                                                                                        >  2:1151850/1‑69 (MQ=255)
                                                                                        atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                       <  1:40953/118‑1 (MQ=255)
                                                                                        atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                       >  2:40953/1‑118 (MQ=255)
                                                                                           atcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCa                                                  >  2:105736/1‑120 (MQ=255)
                                                                                           atcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCa                                                  <  1:105736/120‑1 (MQ=255)
                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAAccag                                          >  2:636319/1‑118 (MQ=255)
                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCg                                          <  1:636319/120‑1 (MQ=255)
                                                                                                                gCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                              <  2:832463/55‑1 (MQ=255)
                                                                                                                gCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                              >  1:832463/1‑55 (MQ=255)
                                                                                                                        ggTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGcaca  <  2:1448451/139‑1 (MQ=255)
                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTc                                                                                                  >  1:35004/1‑38 (MQ=255)
                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTc                                                                                                  <  2:35004/38‑1 (MQ=255)
                                                                                                                                     |                                                                                                                             
AAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACA  >  NC_000913/345516‑345774

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