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% 130.9 / NA 41*196* (TAG→TAA) cbrCUPF0167 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/22);  total (19/22)

TGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCT  >  NC_000913/3898474‑3898735
                                                                                                                                       |                                                                                                                              
tGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTa                                                                                                                             <  2:852622/139‑1 (MQ=255)
     aagaTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        <  1:946945/139‑1 (MQ=255)
        aTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaata                                                                                                                     >  1:1018489/1‑139 (MQ=255)
          tgtgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                   <  2:649707/139‑1 (MQ=255)
                  tttCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        <  2:294396/126‑1 (MQ=255)
                  tttCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        >  1:294396/1‑126 (MQ=255)
                    tCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGc                                                                                                         <  1:1070634/139‑1 (MQ=255)
                              aaaTTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTaa                                                                                               <  2:1116529/139‑1 (MQ=255)
                                        cTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                     <  2:1213324/139‑1 (MQ=255)
                                               aaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                              <  2:1008940/139‑1 (MQ=255)
                                                     tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        <  1:1207317/91‑1 (MQ=255)
                                                     tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        >  2:1207317/1‑91 (MQ=255)
                                                     tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                       >  1:35000/1‑108 (MQ=255)
                                                     tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                       <  2:35000/108‑1 (MQ=255)
                                                      gCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        >  1:10692/1‑90 (MQ=255)
                                                      gCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        <  2:10692/90‑1 (MQ=255)
                                                      gCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTc                                                                         >  1:556593/1‑137 (MQ=255)
                                                      gCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTc                                                                         <  2:556593/137‑1 (MQ=255)
                                                                 gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGa                                                                                                     <  2:264274/98‑1 (MQ=255)
                                                                 gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGa                                                                                                     >  1:264274/1‑98 (MQ=255)
                                                                 gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                            <  2:244736/139‑1 (MQ=255)
                                                                       gTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCt                                                      <  1:385331/139‑1 (MQ=255)
                                                                                aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        >  1:888009/1‑64 (MQ=255)
                                                                                aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        <  2:888009/64‑1 (MQ=255)
                                                                                aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCg                                             >  1:2607/1‑139 (MQ=255)
                                                                                       ccGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCgg                                      >  2:362028/1‑139 (MQ=255)
                                                                                                       ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        <  1:143484/41‑1 (MQ=255)
                                                                                                       ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                        >  2:143484/1‑41 (MQ=255)
                                                                                                       ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAc                      >  2:1173329/1‑139 (MQ=255)
                                                                                                       ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAc                      >  1:917404/1‑139 (MQ=255)
                                                                                                          aaGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCAt                   <  2:173654/139‑1 (MQ=255)
                                                                                                                agaCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                  <  1:981062/54‑1 (MQ=255)
                                                                                                                agaCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                  >  2:981062/1‑54 (MQ=255)
                                                                                                                  aCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg           >  2:1145733/1‑139 (MQ=255)
                                                                                                                      tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATg       >  2:958619/1‑139 (MQ=255)
                                                                                                                          gTGATTTTCCGTAATTATTTAAATAATGAGAACAGGACGGAGCGTAATTCAAACATCAGGCGTTATTAGTGAAGCTTATCCAAGTGATAATACTGCAGCTGCGGCTTTATCAGAGATTACCAACAGTGAAGGagggaga   >  1:127977/1‑133 (MQ=255)
                                                                                                                           tGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCt  <  1:1145733/139‑1 (MQ=255)
                                                                                                                             aTTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTca                                                                                           <  2:1052126/48‑1 (MQ=255)
                                                                                                                             aTTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTca                                                                                           >  1:1052126/1‑48 (MQ=255)
                                                                                                                                   cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                              >  2:143460/1‑55 (MQ=255)
                                                                                                                                   cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                              <  1:143460/55‑1 (MQ=255)
                                                                                                                                       |                                                                                                                              
TGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCT  >  NC_000913/3898474‑3898735

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