Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,124	C→T	W481* (TGG→TAG)	cpdB ←	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,124	0	C	T	100.0%	97.4 / NA	31	W481* (TGG→TAG)	cpdB	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/9);  total (22/9)

GCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCAC  >  NC_000913/4434881‑4435362                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               gCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTATTtgat                                                                                                                                                                                                                                         <  1:246783/251‑1 (MQ=255)                    cGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTg                                                                                                                                                                                                                      >  2:98998/1‑251 (MQ=255)                           aCCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGCAACATGGCGTTCGGATCAATCGGCTTGCCCTTAAAGGTCAGGTTCTTAATCCTTTCCGCATTCGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCCCACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTgg                                                                                                                                                                                                               >  1:187777/1‑251 (MQ=255)                                   aTTTGCCGCCGTAAGCGCGATAGTTATTGGGGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCGGTTGGGATCactc                                                                                                                                                                                                       >  2:37912/1‑248 (MQ=255)                                       gccgccGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTggg                                                                                                                                                                                                   >  1:273347/1‑250 (MQ=255)                                                   gcgATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTGCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCcg                                                                                                                                                                                       <  2:107378/251‑1 (MQ=255)                                                     gATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCcgcg                                                                                                                                                                                     >  1:303412/1‑251 (MQ=255)                                                        aggtATTGGAGGCAACAAGGAACATGGCGTTCGGCTCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCATACCATCAATACCATCACCGTTATAAGTGCGGAACCCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCACTCTGGTTAAACTGTCCCGCggcg                                                                                                                                                                                  >  2:109594/4‑249 (MQ=255)                                                                  ggCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTGCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTTGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCaccc                                                                                                                                                                        <  2:99116/251‑4 (MQ=255)                                                                         aggaaaaaggaggtCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGATGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCttt                                                                                                                                                                 >  1:185502/14‑251 (MQ=255)                                                                          ggAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAACCTGTCCCGCGGAGCACTCCAGCCACTCTTTc                                                                                                                                                                >  2:275529/1‑251 (MQ=255)                                                                             aCATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTc                                                                                                                                                                                                                         <  2:284454/191‑1 (MQ=255)                                                                             aCATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTc                                                                                                                                                                                                                         >  1:284454/1‑191 (MQ=255)                                                                               aTGGCGTTCGGATCAATCGGCTTGCCATTAACGGTCAGGTTCTTCATCCTTTCCGCATTGGCATTAATCATCTGGAACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACctc                                                                                                                                                           >  1:185102/1‑251 (MQ=255)                                                                                     ttCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTAcc                                                                                                                                                     >  2:304786/1‑251 (MQ=255)                                                                                                               ggTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTAt                                                                                                                           >  1:2314/1‑251 (MQ=255)                                                                                                                  cAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATCACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTgg                                                                                                                        >  2:244565/1‑251 (MQ=255)                                                                                                                       ttttAATCCTTTCCGCATTGGCATTAATCATCTGGCACTAGCCGGCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGGCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATa                                                                                                                    <  2:265811/248‑1 (MQ=255)                                                                                                                              ccTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAATCTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATaaa                                                                                                            <  1:235397/251‑1 (MQ=255)                                                                                                                                       aTTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcggcgc                                                                                                   >  1:106801/1‑250 (MQ=255)                                                                                                                                            cATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCAtt                                                                                               >  1:216589/1‑250 (MQ=255)                                                                                                                                                     tcTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACAACATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTTACCTCTTTACCGCTGGCTTTCACCACAATAAGCGTATTGGGATAGAGATAAAGATCGGTGGCATTACGGAAggcc                                                                                     >  1:11013/1‑249 (MQ=255)                                                                                                                                                                         cgggcgggCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTc                                                                 >  1:211192/1‑251 (MQ=255)                                                                                                                                                                                                  gATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCAtt                                            <  1:244565/247‑1 (MQ=255)                                                                                                                                                                                                     aaTTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTGGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGATTTCACCACAATCAGCGTATTGGGATAGAGATAACGATCGGCGGCATTACGGCAGGTCAACTGGCCTTTTTACACCTCCACATAGCTTGCCTGGTCATTTTTGCg                                      >  2:198696/1‑250 (MQ=255)                                                                                                                                                                                                       ttCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGacc                                   >  2:260390/1‑251 (MQ=255)                                                                                                                                                                                                            acCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcg                                                                                                       >  1:38208/1‑178 (MQ=255)                                                                                                                                                                                                            acCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATcg                                                                                                       <  2:38208/178‑1 (MQ=255)                                                                                                                                                                                                                      aCATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACgg                    >  1:128869/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGaa       >  2:258438/1‑249 (MQ=255)                                                                                                                                                                                                                                        cGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGGGGTGTTGGGATCATTCTGGTTAAACTGTCCCGCGGATCACTCCAGCCACTCTTTCACCTCTTTACCGCTTGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGTAAGGTCAACTGGCCTGTTTCCTCCTCCACATAGTTGGCCGGGTGTGTTGTGCGACCACCGTCTTTAAACGGTGCGGCAGCTGAAAGCAc  <  1:180636/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               GCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGTCCCGCGGAGCACTCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCAC  >  NC_000913/4434881‑4435362

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