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_000913	3,898,609	0	G	A	100.0%	144.3 / NA	45	*196* (TAG→TAA)	cbrC	UPF0167 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/23);  total (22/23)

TGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG  >  NC_000913/3898484‑3898743                                                                                                                              |                                                                                                                                       tgtgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  1:440925/139‑1 (MQ=255) tgtgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  1:363176/139‑1 (MQ=255)   tgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATga                                                                                                                         >  2:623825/1‑139 (MQ=255)   tgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGT‑‑GGTGATTTTTCGTAATTATTTAAATAATgaga                                                                                                                       >  1:1234049/1‑139 (MQ=255)                aTTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAg                                                                                                            <  1:623825/139‑1 (MQ=255)                    gAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               <  2:408681/132‑1 (MQ=255)                    gAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               >  1:408681/1‑132 (MQ=255)                     aaaTTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  2:245668/1‑114 (MQ=255)                     aaaTTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  1:245668/114‑1 (MQ=255)                         tCTGATCTAGCTAAATGCCTGCTAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTca                                                                                                   <  2:201743/139‑1 (MQ=255)                                       aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     >  1:268535/1‑139 (MQ=255)                                            tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           >  2:260187/1‑96 (MQ=255)                                            tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  1:260187/96‑1 (MQ=255)                                            tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCt                                                                                       <  2:795424/132‑1 (MQ=255)                                            tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCt                                                                                       >  1:795424/1‑132 (MQ=255)                                                        gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  2:123154/79‑1 (MQ=255)                                                        gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  1:123154/1‑79 (MQ=255)                                                        gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             <  2:1131184/114‑1 (MQ=255)                                                        gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             >  1:1131184/1‑114 (MQ=255)                                                        gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    >  1:1439269/1‑139 (MQ=255)                                                        gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACGGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    >  1:270704/1‑139 (MQ=255)                                                           aTTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  2:594973/81‑1 (MQ=255)                                                           aTTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           >  1:594973/1‑81 (MQ=255)                                                             tGTCAGGGTTATCTCTTGCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTCAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    >  2:864383/1‑134 (MQ=255)                                                             tGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    <  1:864383/134‑1 (MQ=255)                                                              gTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCt                                                              <  2:1439269/139‑1 (MQ=255)                                                                   ggTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  1:325022/1‑68 (MQ=255)                                                                   ggTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  2:325022/68‑1 (MQ=255)                                                                        tctcTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCt                                                                                >  1:517955/1‑111 (MQ=255)                                                                        tctcTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCt                                                                                <  2:517955/111‑1 (MQ=255)                                                                                    tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGc                                                <  2:290453/131‑1 (MQ=255)                                                                                    tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGc                                                >  1:290453/1‑131 (MQ=255)                                                                                           tGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGaa                                                                                                                      <  2:1083774/54‑1 (MQ=255)                                                                                           tGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGaa                                                                                                                      >  1:1083774/1‑54 (MQ=255)                                                                                              ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAc                              >  1:964645/1‑139 (MQ=255)                                                                                                     tGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               <  2:1378698/83‑1 (MQ=255)                                                                                                     tGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               >  1:1378698/1‑83 (MQ=255)                                                                                                      gagaCTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCg                      <  2:964645/139‑1 (MQ=255)                                                                                                             tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATg               >  2:352420/1‑139 (MQ=255)                                                                                                              ggggTGATTTTTCGTAATTATTTAGATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    >  2:1022683/1‑85 (MQ=255)                                                                                                              ggggTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc              <  2:268535/139‑1 (MQ=255)                                                                                                              ggggTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCATTATTTCTTAAGCTTAACCg                                                                    <  1:1022683/85‑1 (MQ=255)                                                                                                                   gATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    >  2:1404839/1‑80 (MQ=255)                                                                                                                   gATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    <  1:1404839/80‑1 (MQ=255)                                                                                                                          cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg  <  2:270704/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                       TGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG  >  NC_000913/3898484‑3898743

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: ‑