Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	640,165	C→T	A138V (GCG→GTG)	ahpF →	alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	640,165	0	C	T	100.0%	170.8 / NA	49	A138V (GCG→GTG)	ahpF	alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/26);  total (23/26)

GCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCAGTGTTCGTAAACGGGAAAGAG  >  NC_000913/640034‑640292                                                                                                                                    |                                                                                                                                gCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGTCGTGGTGCAGGTGCTGAAc                                                                                                                          >  1:1106139/1‑139 (MQ=255) gCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAAc                                                                                                                          >  1:1059934/1‑139 (MQ=255)            gtggtCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAAc                                                                                                                          <  1:924066/128‑1 (MQ=255)            gtggtCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAAc                                                                                                                          >  2:924066/1‑128 (MQ=255)                  aTCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGaa                                                                                                         >  2:869308/1‑139 (MQ=255)                       tCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGTCGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCcgc                                                                                                    <  2:1106139/139‑1 (MQ=255)                       tCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCcgc                                                                                                    <  1:37150/139‑1 (MQ=255)                              aagCGCAGTCTCTGCTGGAGCAGATTCGCTATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCaa                                                                                             >  2:480829/1‑139 (MQ=255)                                        tctGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGca                                                                                          <  2:518518/132‑1 (MQ=255)                                        tctGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGca                                                                                          >  1:518518/1‑132 (MQ=255)                                            ctgGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTg                                                                               <  2:771867/139‑1 (MQ=255)                                                aGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCtgat                                                                                                                     >  1:982177/1‑97 (MQ=255)                                                aGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCtgat                                                                                                                     <  2:982177/97‑1 (MQ=255)                                                    gATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACggcggc                                                                       >  1:219901/1‑139 (MQ=255)                                                     aTTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCa                                                                      <  1:549887/139‑1 (MQ=255)                                                                        aTTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCAc                                                   >  1:435954/1‑139 (MQ=255)                                                                        aTTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCAc                                                   >  2:1088329/1‑139 (MQ=255)                                                                            tGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCgcgc                                                                                                  <  2:233854/88‑1 (MQ=255)                                                                            tGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCgcgc                                                                                                  >  1:233854/1‑88 (MQ=255)                                                                            tGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGAt                                               <  1:1088329/139‑1 (MQ=255)                                                                               aTTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGc                                            <  2:1059934/139‑1 (MQ=255)                                                                                ttCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCtg                                                                                                                       >  2:494193/1‑63 (MQ=255)                                                                                ttCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCtg                                                                                                                       <  1:494193/63‑1 (MQ=255)                                                                                           ttACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGc                                >  1:271973/1‑139 (MQ=255)                                                                                           ttACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCT‑AACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCg                               <  1:591148/139‑1 (MQ=255)                                                                                                 gCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCg                          >  2:665606/1‑139 (MQ=255)                                                                                                 gCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCg                          <  1:665606/139‑1 (MQ=255)                                                                                                  ctctctTGCCACAACTGCCCGGACGTGGTGCAGGTGCTg                                                                                                                             <  1:731275/39‑1 (MQ=255)                                                                                                  ctctctTGCCACAACTGCCCGGACGTGGTGCAGGTGCTg                                                                                                                             >  2:731275/1‑39 (MQ=255)                                                                                                  ctctctTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCCCCGATCGCAACGTGATCGGCGTTCCgg                         <  2:582576/139‑1 (MQ=255)                                                                                                    ctctTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCa                       <  1:842675/139‑1 (MQ=255)                                                                                                     tctTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCa                       >  1:1088481/1‑138 (MQ=255)                                                                                                     tctTGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCa                       <  2:1088481/138‑1 (MQ=255)                                                                                                        tGCCACAACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCCGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCAgtgt                   >  1:161731/1‑139 (MQ=255)                                                                                                             caACTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCAGTGTTCGTa              <  1:480829/139‑1 (MQ=255)                                                                                                               aCTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca                                                                                        <  2:740570/63‑1 (MQ=255)                                                                                                               aCTGCCCGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca                                                                                        >  1:740570/1‑63 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCgcgc                                                                                                  >  2:143673/1‑49 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCgcgc                                                                                                  <  1:143673/49‑1 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGca                                                                                          >  2:472757/1‑57 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGca                                                                                          <  1:472757/57‑1 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCa                                                                      <  2:1029053/77‑1 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCa                                                                      >  1:1029053/1‑77 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCa                                           <  1:972488/104‑1 (MQ=255)                                                                                                                   cccGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCa                                           >  2:972488/1‑104 (MQ=255)                                                                                                                     cGGACGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCAGTGTTCGTAAACGGGaa      <  1:919300/139‑1 (MQ=255)                                                                                                                         cGTGGTGCAGGTGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCAGTGTTCGTAAACGGGAAagag  <  2:271973/139‑1 (MQ=255)                                                                                                                              tgCAGGTGCTGAACCTGATGAGCGTACTGAACCcgc                                                                                                    <  2:1111865/36‑1 (MQ=255)                                                                                                                              tgCAGGTGCTGAACCTGATGAGCGTACTGAACCcgc                                                                                                    >  1:1111865/1‑36 (MQ=255)                                                                                                                                    |                                                                                                                                GCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCCATATTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACACTGCAATTGACGGCGGCACCTTCCAGAACGAAATCACCGATCGCAACGTGATGGGCGTTCCGGCAGTGTTCGTAAACGGGAAAGAG  >  NC_000913/640034‑640292

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