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%	233.8 / NA	72	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 (35/37);  total (35/37)

TAATTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGA  >  NC_000913/444506‑444756                                                                                                                                        |                                                                                                                    tAATTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTggg                                                                                                                  >  1:1627467/1‑139 (MQ=255)              tcctccTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGc                                                                                                     <  2:1438316/139‑1 (MQ=255)                  cctcACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 <  2:1065258/139‑1 (MQ=255)                       cccGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                        >  1:941246/1‑127 (MQ=255)                       cccGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                        <  2:941246/127‑1 (MQ=255)                          gCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                         <  1:1573938/139‑1 (MQ=255)                            aaCCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                       <  2:1627467/139‑1 (MQ=255)                                            cAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGTTTCGCATCaaa                                                                                               >  2:1702367/1‑115 (MQ=255)                                            cAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                       <  2:1542638/139‑1 (MQ=255)                                            cAATAAGGCTATCCCTTAATTGTACATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGTTTATTTAACTTTGTAGGTATTATGTTTCGCATCaaa                                                                                               <  1:1702367/115‑1 (MQ=255)                                                       tCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGata                                                            >  2:68931/1‑139 (MQ=255)                                                         ccTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                          <  2:647909/139‑1 (MQ=255)                                                         ccTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                          <  1:68931/139‑1 (MQ=255)                                                                 gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 <  2:1619114/92‑1 (MQ=255)                                                                 gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 >  1:1619114/1‑92 (MQ=255)                                                                 gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                    >  2:730343/1‑121 (MQ=255)                                                                 gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                    <  1:730343/121‑1 (MQ=255)                                                                  tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                         <  2:198490/99‑1 (MQ=255)                                                                  tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                         >  1:198490/1‑99 (MQ=255)                                                                   gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCAt                                                                     >  2:177974/1‑118 (MQ=255)                                                                   gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCAt                                                                     <  1:177974/118‑1 (MQ=255)                                                                   gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAg                                                >  2:119719/1‑139 (MQ=255)                                                                        cTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                        <  2:1351644/78‑1 (MQ=255)                                                                        cTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                        >  1:1351644/1‑78 (MQ=255)                                                                         tGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATc                                          >  2:1659662/1‑139 (MQ=255)                                                                              cGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                      <  2:1419345/106‑1 (MQ=255)                                                                              cGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                      >  1:1419345/1‑106 (MQ=255)                                                                                   aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                       >  2:159917/1‑100 (MQ=255)                                                                                   aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                       <  1:159917/100‑1 (MQ=255)                                                                                     gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCg                                                                                                      >  1:1725503/1‑67 (MQ=255)                                                                                     gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCg                                                                                                      <  2:1725503/67‑1 (MQ=255)                                                                                      cacaATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATg                                                                                           <  2:253726/77‑1 (MQ=255)                                                                                      cacaATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATg                                                                                           >  1:253726/1‑77 (MQ=255)                                                                                           ttAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTgg                                                                                                                   <  2:1664249/48‑1 (MQ=255)                                                                                           ttAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTgg                                                                                                                   >  1:1664249/1‑48 (MQ=255)                                                                                           ttAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcg                        >  2:664347/1‑139 (MQ=255)                                                                                           ttAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcg                        <  1:664347/139‑1 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 >  1:72468/1‑56 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 <  2:72468/56‑1 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATg                                                                                           <  2:1523724/62‑1 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATg                                                                                           >  1:1523724/1‑62 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                       <  1:1498856/66‑1 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                       >  2:1498856/1‑66 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                    <  1:1485630/85‑1 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                    >  2:1485630/1‑85 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCttt                                                                 <  1:1564386/88‑1 (MQ=255)                                                                                                     cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCttt                                                                 >  2:1564386/1‑88 (MQ=255)                                                                                                             tGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGtctc                                                     <  2:759899/92‑1 (MQ=255)                                                                                                             tGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGtctc                                                     >  1:759899/1‑92 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                       >  1:1092165/1‑54 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                       <  2:1092165/54‑1 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                          <  2:1212862/83‑1 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                          <  1:46580/83‑1 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                          >  1:1212862/1‑83 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                          >  2:46580/1‑83 (MQ=255)                                                                                                                 ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGa  >  2:1397000/1‑139 (MQ=255)                                                                                                                      gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                            <  2:1341869/92‑1 (MQ=255)                                                                                                                      gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                            >  1:1341869/1‑92 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 >  2:1450515/1‑38 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                 <  1:1450515/38‑1 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAACTGAAGGGAGAGGAAAAGATGcc                                                                       >  1:1208317/1‑64 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGaa                                                                               <  2:1266433/56‑1 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGaa                                                                               >  1:1266433/1‑56 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGAt                                                                          >  1:734692/1‑61 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGAt                                                                          <  2:734692/61‑1 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                       <  2:1208317/64‑1 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGatat                                                           >  1:764158/1‑76 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGatat                                                           <  2:764158/76‑1 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGa                                                  <  2:368735/85‑1 (MQ=255)                                                                                                                       tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGa                                                  >  1:368735/1‑85 (MQ=255)                                                                                                                              tttAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAgg                                    >  2:909530/1‑92 (MQ=255)                                                                                                                              tttAACTTTGTGGGTAGTAAGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGAGGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAgg                                    <  1:909530/92‑1 (MQ=255)                                                                                                                                        |                                                                                                                    TAATTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGA  >  NC_000913/444506‑444756

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