Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,073,232	A→G	F260F (TTT→TTC)	rutA ←	pyrimidine oxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,073,232	0	A	G	100.0%	153.0 / NA	43	F260F (TTT→TTC)	rutA	pyrimidine oxygenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/22);  total (21/22)

TTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAG  >  NC_000913/1073098‑1073367                                                                                                                                       |                                                                                                                                        ttGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGaacaa                                                                                                                                    >  1:47386/1‑140 (MQ=255)  tGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGaacaac                                                                                                                                   <  1:1000305/140‑1 (MQ=255)       tCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATa                                                                                                                              <  1:749691/140‑1 (MQ=255)          gTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAacaaca                                                                                                                                  >  2:314453/1‑133 (MQ=255)          gTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAacaaca                                                                                                                                  <  1:314453/133‑1 (MQ=255)          gTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAaga                                                                                                                           >  1:727181/1‑140 (MQ=255)           tACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAagag                                                                                                                          >  1:810525/1‑140 (MQ=255)           tACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAagag                                                                                                                          >  2:860076/1‑140 (MQ=255)                     ggTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTcgc                                                                                                                >  1:354493/1‑140 (MQ=255)                     ggTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTcgc                                                                                                                >  1:1077605/1‑140 (MQ=255)                     ggTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTcgc                                                                                                                >  2:1065614/1‑140 (MQ=255)                     ggTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTcgc                                                                                                                >  1:204877/1‑140 (MQ=255)                       tATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGt                                                                                                              <  2:589026/140‑1 (MQ=255)                                   tttGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGctct                                                                                                  >  2:601755/1‑140 (MQ=255)                                   tttGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGctct                                                                                                  <  1:682080/140‑1 (MQ=255)                                       ttCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCgg                                                                                              >  1:289852/1‑140 (MQ=255)                                       ttCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCgg                                                                                              <  2:542364/140‑1 (MQ=255)                                       ttCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCgg                                                                                              <  2:486292/140‑1 (MQ=255)                                       ttCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCgg                                                                                              >  2:366172/1‑140 (MQ=255)                                        tCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGc                                                                                             <  1:601755/140‑1 (MQ=255)                                                ccAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCa                                                                                     >  2:597546/1‑140 (MQ=255)                                                   gCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATcc                                                                                  <  2:47386/140‑1 (MQ=255)                                                              ctTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGTGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTc                                                                       <  2:445541/140‑1 (MQ=255)                                                              ctTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTc                                                                       <  1:597546/140‑1 (MQ=255)                                                              ctTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTc                                                                       >  2:270114/1‑140 (MQ=255)                                                              ctTCATCCGCGCCCGCTTTGTAGTGTTACCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTc                                                                       <  2:453307/140‑1 (MQ=255)                                                                  aTCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCgggg                                                                   <  2:107253/140‑1 (MQ=255)                                                                    ccGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCg                                                                 <  1:270114/140‑1 (MQ=255)                                                                      gcgcCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaa                                                               >  1:198099/1‑140 (MQ=255)                                                                       cgcCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaaa                                                              >  1:165390/1‑140 (MQ=255)                                                                         cccGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGc                                                            >  1:434871/1‑140 (MQ=255)                                                                         cccGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGc                                                            >  1:944219/1‑140 (MQ=255)                                                                         cccGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGc                                                            >  1:748312/1‑140 (MQ=255)                                                                          ccGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGcc                                                           <  1:366172/140‑1 (MQ=255)                                                                           cGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCg                                                          <  2:198099/140‑1 (MQ=255)                                                                             cTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTg                                                        <  2:727181/140‑1 (MQ=255)                                                                                tGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGgt                                                     <  2:748312/140‑1 (MQ=255)                                                                                tGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGgt                                                     <  2:434871/140‑1 (MQ=255)                                                                                       ttCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTa                                              <  1:1065614/140‑1 (MQ=255)                                                                                                  cgcgAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTATGCCg                                   >  2:314995/1‑140 (MQ=255)                                                                                                     gAGCGGCATCGTCGGTTTCATCGGCAATCACCATGAACAACACATAAGAGCCATCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGaaa                                <  2:944219/140‑1 (MQ=255)                                                                                                                    tttCATCGGCAATCACCATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCg                 >  2:1071131/1‑140 (MQ=255)                                                                                                                                   ccATGAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgag  <  2:165390/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        TTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCGTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAG  >  NC_000913/1073098‑1073367

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