Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,025 G→A *118* (TAG→TAA ybaA → DUF1428 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,0250GA100.0% 136.8 / NA 43*118* (TAG→TAAybaADUF1428 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/21);  total (22/21)

CCCTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATG  >  NC_000913/476890‑477153
                                                                                                                                       |                                                                                                                                
cccTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTa                                                                                                                               <  1:59420/139‑1 (MQ=255)
   ttCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCg                                                                                                                            >  1:1038696/1‑139 (MQ=255)
   ttCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCg                                                                                                                            >  2:531966/1‑139 (MQ=255)
             gTCCGCGACGCTGCTACTCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCtg                                                                                                                  >  2:989981/1‑139 (MQ=255)
                          ctAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGt                                                                                                     >  1:1471169/1‑139 (MQ=255)
                                       gatgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAAtt                                                                                        >  1:304461/1‑139 (MQ=255)
                                               aCCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                <  2:304461/139‑1 (MQ=255)
                                                    cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                           <  2:1471169/139‑1 (MQ=255)
                                                    cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                           <  2:658013/139‑1 (MQ=255)
                                                         gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                      <  1:531966/139‑1 (MQ=255)
                                                         gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                      <  2:1038696/139‑1 (MQ=255)
                                                                ttCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTa                                                               >  1:1603657/1‑139 (MQ=255)
                                                                      gAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAt                                                         <  2:1156700/139‑1 (MQ=255)
                                                                        gTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtct                                                       >  1:580776/1‑139 (MQ=255)
                                                                           cATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttct                                                    >  2:1464107/1‑139 (MQ=255)
                                                                           cATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttct                                                    >  2:1478323/1‑139 (MQ=255)
                                                                              gCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                      <  1:1198410/118‑1 (MQ=255)
                                                                              gCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                      >  2:1198410/1‑118 (MQ=255)
                                                                                   ttGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                >  2:513826/1‑103 (MQ=255)
                                                                                   ttGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                <  1:513826/103‑1 (MQ=255)
                                                                                     gACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATg                                          <  1:1464107/139‑1 (MQ=255)
                                                                                      aCGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGc                                         <  2:307622/139‑1 (MQ=255)
                                                                                                      cTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCTCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagccag                         >  1:708321/1‑139 (MQ=255)
                                                                                                            cGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                           >  1:626471/1‑83 (MQ=255)
                                                                                                            cGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                           <  2:626471/83‑1 (MQ=255)
                                                                                                              gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGc                                                                       >  1:242507/1‑85 (MQ=255)
                                                                                                              gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGc                                                                       <  2:242507/85‑1 (MQ=255)
                                                                                                              gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCgg                                               >  2:164339/1‑109 (MQ=255)
                                                                                                              gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATGCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCgg                                               <  1:164339/109‑1 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGccc                                                                                               >  1:1007042/1‑60 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGccc                                                                                               <  2:1007042/60‑1 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaa                                                                                           >  2:1188334/1‑64 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaa                                                                                           <  1:1188334/64‑1 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                      <  2:547420/85‑1 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                      >  1:547420/1‑85 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAg                                                                   >  2:580217/1‑88 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAg                                                                   <  1:580217/88‑1 (MQ=255)
                                                                                                               aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCg                >  2:561061/1‑139 (MQ=255)
                                                                                                                ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGccc                                                                     >  1:704658/1‑85 (MQ=255)
                                                                                                                ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGccc                                                                     <  2:704658/85‑1 (MQ=255)
                                                                                                                ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGt               >  2:498845/1‑139 (MQ=255)
                                                                                                                       cAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCAc        <  2:580776/139‑1 (MQ=255)
                                                                                                                             tcGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATg  <  1:561061/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                
CCCTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATG  >  NC_000913/476890‑477153

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: