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_000913444,6410CT100.0% 155.2 / NA 49intergenic (‑126/‑42)panE/yajQ2‑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 (24/25);  total (24/25)

CCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAC  >  NC_000913/444523‑444775
                                                                                                                      |                                                                                                                                      
cctcACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGCCTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                    <  1:876492/139‑1 (MQ=255)
cctcACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                    <  1:1034569/139‑1 (MQ=255)
            ccACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGga                                                                                                        <  1:24800/139‑1 (MQ=255)
                       gtTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGga                                                                                                        <  1:713641/128‑1 (MQ=255)
                       gtTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGga                                                                                                        >  2:713641/1‑128 (MQ=255)
                           aaTAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                         <  1:1344807/139‑1 (MQ=255)
                                                 gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTgg                                                                                                                                      <  1:1370556/72‑1 (MQ=255)
                                                 gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTgg                                                                                                                                      >  2:1370556/1‑72 (MQ=255)
                                                 gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAg                                                                   >  1:192501/1‑139 (MQ=255)
                                                     gCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                               <  1:1392814/139‑1 (MQ=255)
                                                       tGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                         >  1:653392/1‑111 (MQ=255)
                                                       tGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                         <  2:653392/111‑1 (MQ=255)
                                                       tGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATc                                                             <  2:577116/139‑1 (MQ=255)
                                                              aCTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTGATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGa                                                      <  2:192501/139‑1 (MQ=255)
                                                                 aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGc                                                   >  1:1046958/1‑139 (MQ=255)
                                                                   gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                    <  2:525582/72‑1 (MQ=255)
                                                                   gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                    >  1:525582/1‑72 (MQ=255)
                                                                   gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                          <  2:1395344/98‑1 (MQ=255)
                                                                   gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                          >  1:1395344/1‑98 (MQ=255)
                                                                            aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCttt                                                                                    <  2:659439/95‑1 (MQ=255)
                                                                            aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCttt                                                                                    >  1:659439/1‑95 (MQ=255)
                                                                              ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGta                                                                                                                                   >  2:683119/1‑46 (MQ=255)
                                                                              ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGta                                                                                                                                   <  1:683119/46‑1 (MQ=255)
                                                                                 tACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                                           <  2:328071/51‑1 (MQ=255)
                                                                                 tACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCtt                                                                                                                           >  1:328071/1‑51 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGtatt                                                                                                                                 <  1:450313/43‑1 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGtatt                                                                                                                                 >  2:450313/1‑43 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                          >  1:177460/1‑82 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                          <  2:177460/82‑1 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCtt                                                                                     >  1:534308/1‑87 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCtt                                                                                     <  2:534308/87‑1 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGata                                                                               >  2:1078129/1‑93 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGata                                                                               <  1:1078129/93‑1 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                 >  2:1350026/1‑139 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                 >  1:236716/1‑139 (MQ=255)
                                                                                   cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                 >  2:724816/1‑139 (MQ=255)
                                                                                      ccTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAAcgc                              >  2:1038097/1‑139 (MQ=255)
                                                                                               ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCttt                                                                                    <  2:250973/76‑1 (MQ=255)
                                                                                               ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCttt                                                                                    >  1:250973/1‑76 (MQ=255)
                                                                                                gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGaa                    >  2:3807/1‑139 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                       >  2:1308586/1‑67 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                       <  1:1308586/67‑1 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                               <  2:516750/91‑1 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                               >  1:516750/1‑91 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCtt                                                           >  2:1113746/1‑95 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCtt                                                           <  1:1113746/95‑1 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGa                                                      >  2:1001285/1‑100 (MQ=255)
                                                                                                     tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGa                                                      <  1:1001285/100‑1 (MQ=255)
                                                                                                                  tttGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAc  <  1:1038097/139‑1 (MQ=255)
                                                                                                                      |                                                                                                                                      
CCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAC  >  NC_000913/444523‑444775

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