Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,385,857 C→T *91* (TAG→TAA yhcO ← putative barnase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,385,8570CT100.0% 138.6 / NA 43*91* (TAG→TAAyhcOputative barnase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/24);  total (19/24)

TAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCAGG  >  NC_000913/3385759‑3385992
                                                                                                  |                                                                                                                                       
tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                               <  1:377711/109‑1 (MQ=255)
tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGaa                                                                                                                               >  2:377711/1‑109 (MQ=255)
 aGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTc                                                                                                >  1:1094914/1‑139 (MQ=255)
   cGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCct                                                                                              <  1:866280/139‑1 (MQ=255)
        acacCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGcc                                                                                         <  2:1372037/139‑1 (MQ=255)
                    cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                               >  2:57665/1‑105 (MQ=255)
                    cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                               <  1:57665/105‑1 (MQ=255)
                    cTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGc                                                                             <  2:1314445/139‑1 (MQ=255)
                       cGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAat                                                                          >  1:1231128/1‑139 (MQ=255)
                            gAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTaa                                                                     <  1:936178/139‑1 (MQ=255)
                             aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                               >  2:293196/1‑96 (MQ=255)
                             aaCTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATg                                                                                                               <  1:293196/96‑1 (MQ=255)
                               cTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGAc                                                                                                                                  >  2:317963/1‑75 (MQ=255)
                               cTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGAc                                                                                                                                  <  1:317963/75‑1 (MQ=255)
                                 gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACgcgc                                                                >  1:380709/1‑139 (MQ=255)
                                   aCAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACgcgc                                                                <  1:1339682/137‑1 (MQ=255)
                                   aCAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACgcgc                                                                >  2:1339682/1‑137 (MQ=255)
                                     aaaTAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCaaa                                                            <  2:523674/139‑1 (MQ=255)
                                           aCCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCg                                                      <  2:1094914/139‑1 (MQ=255)
                                            cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaa                                                                                                                        <  2:1239236/72‑1 (MQ=255)
                                            cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTaaa                                                                                                                        >  1:1239236/1‑72 (MQ=255)
                                            cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                     >  1:1108080/1‑139 (MQ=255)
                                            cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                                     >  1:254304/1‑139 (MQ=255)
                                                 aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGccc                                                                                                            <  1:206230/79‑1 (MQ=255)
                                                 aTCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGccc                                                                                                            >  2:206230/1‑79 (MQ=255)
                                                              acgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                                   <  2:509079/139‑1 (MQ=255)
                                                              acgaTATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGa                                   <  2:254304/139‑1 (MQ=255)
                                                                  tatTGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGaa                               <  1:1338857/139‑1 (MQ=255)
                                                                    ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                        >  1:410925/1‑64 (MQ=255)
                                                                    ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                                        <  2:410925/64‑1 (MQ=255)
                                                                    ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAAc                              >  1:687485/1‑138 (MQ=255)
                                                                    ttGCCCCCGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAAc                              <  2:687485/138‑1 (MQ=255)
                                                                           cGGTTCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCaa                      <  2:380709/139‑1 (MQ=255)
                                                                              ttCGGGGGCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAAtct                   >  1:281944/1‑139 (MQ=255)
                                                                                 ggtggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCaa                <  1:2656/136‑1 (MQ=255)
                                                                                 gggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaaa                                                                                 >  1:730836/1‑74 (MQ=255)
                                                                                 gggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaaa                                                                                 <  2:730836/74‑1 (MQ=255)
                                                                                  ggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctct                                         <  2:308012/113‑1 (MQ=255)
                                                                                  ggggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtctct                                         >  1:308012/1‑113 (MQ=255)
                                                                                   gggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                                >  1:599368/1‑73 (MQ=255)
                                                                                   gggCTTTTTTTTGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAAc                                                                                <  2:599368/73‑1 (MQ=255)
                                                                                        tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                                    >  2:32976/1‑48 (MQ=255)
                                                                                        tttttttGCGTTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                                    <  1:32976/48‑1 (MQ=255)
                                                                                               gcgtTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCAgg  >  1:1164849/1‑139 (MQ=255)
                                                                                                  |                                                                                                                                       
TAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAATCTCAAGTGGTAGCGGCAGG  >  NC_000913/3385759‑3385992

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