Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	444,641	C→T	intergenic (‑126/‑42)	panE ← / → yajQ	2‑dehydropantoate reductase, NADPH‑specific/phage Phi6 host factor, ATP/GTP binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	444,641	0	C	T	100.0%	142.2 / NA	45	intergenic (‑126/‑42)	panE/yajQ	2‑dehydropantoate reductase, NADPH‑specific/phage Phi6 host factor, ATP/GTP binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/20);  total (25/20)

ATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAC  >  NC_000913/444514‑444775                                                                                                                                |                                                                                                                                       aTTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGc                                                                                                                             >  1:646829/1‑139 (MQ=255)   tGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                                           <  2:438854/139‑1 (MQ=255)                            cGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGaa                                                                                                  <  1:51566/139‑1 (MQ=255)                                     aaTAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                         <  2:646829/139‑1 (MQ=255)                                               tCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGata                                                                               >  1:151534/1‑139 (MQ=255)                                               tCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGata                                                                               >  2:308394/1‑139 (MQ=255)                                                                           aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                    <  1:112420/74‑1 (MQ=255)                                                                           aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                    >  2:112420/1‑74 (MQ=255)                                                                           aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGc                                                   >  2:544250/1‑139 (MQ=255)                                                                             gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                            <  2:772674/80‑1 (MQ=255)                                                                             gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                            >  1:772674/1‑80 (MQ=255)                                                                             gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                               <  1:492628/125‑1 (MQ=255)                                                                             gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                               >  2:492628/1‑125 (MQ=255)                                                                                   ttAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATc                                                                                                                     >  2:701880/1‑65 (MQ=255)                                                                                   ttAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATc                                                                                                                     <  1:701880/65‑1 (MQ=255)                                                                                      aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGc                                                                                                                        >  2:76449/1‑59 (MQ=255)                                                                                      aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGc                                                                                                                        <  1:76449/59‑1 (MQ=255)                                                                                        ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                          <  2:100001/87‑1 (MQ=255)                                                                                        ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                          >  1:100001/1‑87 (MQ=255)                                                                                           tACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGAt                                   >  1:585269/1‑139 (MQ=255)                                                                                             cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                          <  2:90862/66‑1 (MQ=255)                                                                                             cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                          >  1:90862/1‑66 (MQ=255)                                                                                             cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGga                                                                                                        <  2:591714/68‑1 (MQ=255)                                                                                             cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGatat                                                                              >  1:284621/1‑94 (MQ=255)                                                                                             cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGatat                                                                              <  2:284621/94‑1 (MQ=255)                                                                                             cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCATCGCATCAAAAATGAAGGga                                                                                                        >  1:591714/1‑68 (MQ=255)                                                                                                        aGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                      <  1:544250/139‑1 (MQ=255)                                                                                                          gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGaa                    >  1:783698/1‑139 (MQ=255)                                                                                                          gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGaa                    >  2:404879/1‑139 (MQ=255)                                                                                                              gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAATGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGtt                                                                 <  2:700123/90‑1 (MQ=255)                                                                                                              gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGtt                                                                 >  1:700123/1‑90 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCaaa                                                                                                                  <  2:821828/40‑1 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCaaa                                                                                                                  >  1:821828/1‑40 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTc                                                                                   <  2:418286/71‑1 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTc                                                                                   >  1:418286/1‑71 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcg                                           <  2:368634/111‑1 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcg                                           >  1:368634/1‑111 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGa               >  1:763165/1‑139 (MQ=255)                                                                                                               tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGa               >  2:797689/1‑139 (MQ=255)                                                                                                                  gTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGatat                                                                              <  2:501298/73‑1 (MQ=255)                                                                                                                  gTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGatat                                                                              >  1:501298/1‑73 (MQ=255)                                                                                                                  gTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTc            <  1:308394/139‑1 (MQ=255)                                                                                                                    tATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTccc          >  2:3825/1‑139 (MQ=255)                                                                                                                         aaCTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGtttt     >  2:632288/1‑139 (MQ=255)                                                                                                                            tttGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAc  <  1:632288/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                       ATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAC  >  NC_000913/444514‑444775

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