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

AGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG  >  NC_000913/3898488‑3898743
                                                                                                                         |                                                                                                                                      
agaATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATgaga                                                                                                                       <  2:462173/139‑1 (MQ=255)
                    tCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTca                                                                                                   <  1:745722/139‑1 (MQ=255)
                      tGATCTAGCTAAATGCCAGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTcaca                                                                                                 >  2:1045251/1‑139 (MQ=255)
                       gATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGTGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               >  1:1138213/1‑124 (MQ=255)
                       gATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               <  2:1138213/124‑1 (MQ=255)
                          cTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             <  2:596639/139‑1 (MQ=255)
                                 aaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                                      <  2:1290319/139‑1 (MQ=255)
                                  aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGCGACAGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGt                                                                                                         >  1:951309/1‑119 (MQ=255)
                                  aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     >  2:945136/1‑139 (MQ=255)
                                  aTGCCTGCAAAAGGGTAGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGt                                                                                                         <  2:951309/119‑1 (MQ=255)
                                       tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAattatt                                                                                                                                   <  2:714477/88‑1 (MQ=255)
                                       tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAattatt                                                                                                                                   >  1:714477/1‑88 (MQ=255)
                                       tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           >  1:1088556/1‑96 (MQ=255)
                                       tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  2:1088556/96‑1 (MQ=255)
                                       tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                                      <  1:5125/133‑1 (MQ=255)
                                       tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCtt                                                                                      >  2:5125/1‑133 (MQ=255)
                                               ggtggtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGa                                                                                                             >  1:1414213/1‑102 (MQ=255)
                                               ggtggtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGa                                                                                                             <  2:1414213/102‑1 (MQ=255)
                                                   gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  1:669287/79‑1 (MQ=255)
                                                   gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  2:669287/1‑79 (MQ=255)
                                                   gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    <  2:791255/139‑1 (MQ=255)
                                                                  aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaat                                                                                                                              <  1:53018/66‑1 (MQ=255)
                                                                  aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaat                                                                                                                              >  2:53018/1‑66 (MQ=255)
                                                                   tctcTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               <  2:473634/112‑1 (MQ=255)
                                                                   tctcTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               >  1:473634/1‑112 (MQ=255)
                                                                   tctcTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               <  2:415318/112‑1 (MQ=255)
                                                                   tctcTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               >  1:415318/1‑112 (MQ=255)
                                                                    ctctTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGc                                                                                                                 >  2:150818/1‑77 (MQ=255)
                                                                    ctctTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGc                                                                                                                 <  1:150818/77‑1 (MQ=255)
                                                                    ctctTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGc                                                   >  1:69375/1‑139 (MQ=255)
                                                                        tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             <  2:1302661/93‑1 (MQ=255)
                                                                        tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             >  1:1302661/1‑93 (MQ=255)
                                                                         ccGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCt                                                                                       <  1:47660/98‑1 (MQ=255)
                                                                         ccGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCt                                                                                       >  2:47660/1‑98 (MQ=255)
                                                                            cTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCtt                                           >  1:644009/1‑139 (MQ=255)
                                                                               tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               <  2:985657/68‑1 (MQ=255)
                                                                               tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               >  1:985657/1‑68 (MQ=255)
                                                                                 tcCACTGCGGCAAGCTGTGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             <  1:528807/84‑1 (MQ=255)
                                                                                 tcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             >  2:528807/1‑84 (MQ=255)
                                                                                  ccACTGCGGCAAGCTGAGACTGGGTGGTGAGTTTTCATAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTCAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAg                                     <  1:1621264/139‑1 (MQ=255)
                                                                                     cTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGa                                  <  2:69375/139‑1 (MQ=255)
                                                                                         ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           >  2:1255345/1‑46 (MQ=255)
                                                                                         ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  1:1255345/46‑1 (MQ=255)
                                                                                            aaGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCAt                           <  1:733284/139‑1 (MQ=255)
                                                                                               cTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAg                        <  2:644009/139‑1 (MQ=255)
                                                                                                        tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTa                                                                                                        >  1:772011/1‑50 (MQ=255)
                                                                                                        tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTa                                                                                                        <  2:772011/50‑1 (MQ=255)
                                                                                                              gATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTg         <  1:945136/139‑1 (MQ=255)
                                                                                                               aTTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAtt                                <  1:1445381/115‑1 (MQ=255)
                                                                                                               aTTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAtt                                >  2:1445381/1‑115 (MQ=255)
                                                                                                                tttttCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCa       >  1:1668258/1‑139 (MQ=255)
                                                                                                                tttttCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCa       >  2:829713/1‑139 (MQ=255)
                                                                                                                   ttCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca    <  2:1668258/139‑1 (MQ=255)
                                                                                                                     cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg  <  1:1045251/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                      
AGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG  >  NC_000913/3898488‑3898743

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