Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,898,609 G→A *196* (TAG→TAA cbrC → UPF0167 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,898,6090GA100.0% 96.3 / NA 31*196* (TAG→TAAcbrCUPF0167 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/15);  total (16/15)

TGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCG  >  NC_000913/3898486‑3898723
                                                                                                                           |                                                                                                                  
tgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATga                                                                                                     >  2:476482/1‑139 (MQ=255)
     aTTTCTGTATTATAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACa                                                                                                <  1:327844/139‑1 (MQ=255)
      tttCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAg                                                                                               >  1:1039808/1‑139 (MQ=255)
          ggTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                           <  2:839874/139‑1 (MQ=255)
                        tGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTcaca                                                                             <  2:657566/139‑1 (MQ=255)
                                  aaaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAattatt                                                                                                               >  1:110806/1‑95 (MQ=255)
                                  aaaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAattatt                                                                                                               <  2:110806/95‑1 (MQ=255)
                                    aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGc                                                                                             >  2:771801/1‑111 (MQ=255)
                                    aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGc                                                                                             <  1:771801/111‑1 (MQ=255)
                                    aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                 >  2:1284141/1‑139 (MQ=255)
                                     tGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                            <  2:934614/95‑1 (MQ=255)
                                     tGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                            >  1:934614/1‑95 (MQ=255)
                                         tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                           >  1:92466/1‑108 (MQ=255)
                                         tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                           <  2:92466/108‑1 (MQ=255)
                                         tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                 <  1:539935/134‑1 (MQ=255)
                                         tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                 >  2:539935/1‑134 (MQ=255)
                                         tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCt                                                            >  1:243819/1‑139 (MQ=255)
                                                    ggtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTc                                                             <  1:660833/127‑1 (MQ=255)
                                                    ggtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTc                                                             >  2:660833/1‑127 (MQ=255)
                                                                    aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATcc                                                                        >  2:1274231/1‑100 (MQ=255)
                                                                    aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATcc                                                                        <  1:1274231/100‑1 (MQ=255)
                                                                              cTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCtt                       >  2:1136122/1‑139 (MQ=255)
                                                                                         gCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                       <  2:1307233/48‑1 (MQ=255)
                                                                                         gCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                       >  1:1307233/1‑48 (MQ=255)
                                                                                               aGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                      >  2:712196/1‑59 (MQ=255)
                                                                                               aGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                      <  1:712196/59‑1 (MQ=255)
                                                                                                   gagaCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCg  <  1:476482/139‑1 (MQ=255)
                                                                                                      aCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                           <  2:440479/47‑1 (MQ=255)
                                                                                                      aCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                           >  1:440479/1‑47 (MQ=255)
                                                                                                          tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                <  2:1206684/86‑1 (MQ=255)
                                                                                                          tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                >  1:1206684/1‑86 (MQ=255)
                                                                                                                       cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                  <  2:1944/55‑1 (MQ=255)
                                                                                                                       cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                  >  1:1944/1‑55 (MQ=255)
                                                                                                                           |                                                                                                                  
TGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCG  >  NC_000913/3898486‑3898723

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