Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,073,080	T→A	D311V (GAT→GTT)	rutA ←	pyrimidine oxygenase

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

TTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCG  >  NC_000913/1072943‑1073208                                                                                                                                          |                                                                                                                                 ttCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACg                                                                                                                                >  2:874491/1‑140 (MQ=255)           aaaaTCGTCGAAGGTTAACAGCACGCCTTCGTCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGAc                                                                                                                      >  1:472190/1‑140 (MQ=255)             aaTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGa                                                                                                                    <  1:264713/140‑1 (MQ=255)              aTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGaa                                                                                                                   >  1:127865/1‑140 (MQ=255)                cgtcgAAGGTTAACAGCACGCCTTCGTCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACg                                                                                                                 <  2:472190/140‑1 (MQ=255)                 gtcgAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGt                                                                                                                <  2:765741/140‑1 (MQ=255)                  tcgAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGtt                                                                                                               <  2:200410/140‑1 (MQ=255)                                  cGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTg                                                                                                              <  2:766434/125‑1 (MQ=255)                                  cGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTg                                                                                                              >  1:766434/1‑125 (MQ=255)                                  cGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAg                                                                                               >  2:759120/1‑140 (MQ=255)                                                ggCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTg                                                                                 >  2:617459/1‑140 (MQ=255)                                                 gCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGa                                                                                <  1:360261/140‑1 (MQ=255)                                                                      aCATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCtt                                                           <  1:1014903/140‑1 (MQ=255)                                                                         tGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAAc                                                        <  2:127865/140‑1 (MQ=255)                                                                          gcgcgcGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACg                                                       >  2:304260/1‑140 (MQ=255)                                                                           cgcgcgACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGc                                                      >  1:253713/1‑140 (MQ=255)                                                                            gcgcgACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGcc                                                     >  1:636253/1‑140 (MQ=255)                                                                            gcgcgACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGcc                                                     >  1:361458/1‑140 (MQ=255)                                                                                gACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAAc                                                        <  2:484403/133‑1 (MQ=255)                                                                                gACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAAc                                                        >  1:484403/1‑133 (MQ=255)                                                                                        cGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCcgcg                                         >  1:497665/1‑140 (MQ=255)                                                                                         gTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCgcgc                                        <  1:371351/140‑1 (MQ=255)                                                                                          tAAGAACCGCCTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGcc                                       <  2:497665/140‑1 (MQ=255)                                                                                          tAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGcc                                       <  1:304260/140‑1 (MQ=255)                                                                                           aaGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGccc                                      <  2:636253/140‑1 (MQ=255)                                                                                           aaGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGccc                                      <  1:799827/140‑1 (MQ=255)                                                                                           aaGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGccc                                      <  2:361458/140‑1 (MQ=255)                                                                                                 cGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTg                                >  2:149990/1‑140 (MQ=255)                                                                                                 cGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTg                                >  1:948132/1‑140 (MQ=255)                                                                                                   aCTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTa                              <  2:948132/140‑1 (MQ=255)                                                                                                    cTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAg                             <  2:701427/140‑1 (MQ=255)                                                                                                        cGTCCCCATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGtt                         <  1:158594/140‑1 (MQ=255)                                                                                                            cccATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCa                     >  2:886949/1‑140 (MQ=255)                                                                                                              cATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAg                             <  1:887088/130‑1 (MQ=255)                                                                                                              cATATTGATGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAg                             >  2:887088/1‑130 (MQ=255)                                                                                                                     aTGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGcgcg            >  2:403656/1‑140 (MQ=255)                                                                                                                     aTGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGcgcg            >  2:363229/1‑140 (MQ=255)                                                                                                                     aTGTTTACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGcgcg            >  1:336895/1‑140 (MQ=255)                                                                                                                          tACCGCCGAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCgg       >  2:234733/1‑140 (MQ=255)                                                                                                                             cgccgAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGa           >  2:624684/1‑133 (MQ=255)                                                                                                                             cgccgAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGa           <  1:624684/133‑1 (MQ=255)                                                                                                                             cgccgAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCAt    >  2:252208/1‑140 (MQ=255)                                                                                                                               ccgAAGTGGGAACGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCAtcg  <  1:149990/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                 TTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTCTGACTTTGTTCGGTTAGCCAGCTTAACGCCTCTTCATCCGCGCCCGCTTTGTAGTGTTCCCATTTGGCGCGAGCGGCATCG  >  NC_000913/1072943‑1073208

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