Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,684	+A	coding (400/501 nt)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,681	1	.	A	100.0%	167.0 / NA	54	Q133K (CAA→AAA)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (32/22);  total (32/22)

ACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGG  >  NC_000913/4435607‑4435750                                                                            |                                                                      aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a                                                                      >  1:945321/1‑76 (MQ=255)  cGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa                                                                     <  1:4460304/76‑1 (MQ=255)    aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑caaa                                                                   >  2:108125/1‑75 (MQ=255)     agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aaaac                                                                  >  2:606121/1‑74 (MQ=255)      gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aaaaca                                                                 >  1:1365402/1‑73 (MQ=255)          tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg                                                              >  2:1244316/1‑76 (MQ=255)          tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg                                                              <  1:4029467/76‑1 (MQ=255)          tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg                                                              >  1:2080573/1‑76 (MQ=255)          tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg                                                              >  1:1433928/1‑76 (MQ=255)          tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg                                                              >  1:2065924/1‑76 (MQ=255)            cccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGAt                                                            <  2:2732472/76‑1 (MQ=255)              cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTg                                                          <  1:2649501/76‑1 (MQ=255)               aTAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGt                                                         <  2:505816/76‑1 (MQ=255)                 aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTAt                                                       >  2:4505761/1‑76 (MQ=255)                  aaaaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATg                                                      >  2:768138/3‑76 (MQ=255)                  aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATg                                                      >  1:495869/1‑76 (MQ=255)                  aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATg                                                      >  2:2322147/1‑76 (MQ=255)                      tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTg                                                  >  2:2811093/1‑76 (MQ=255)                       aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTgg                                                 <  1:3930164/76‑1 (MQ=255)                          tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGat                                              <  1:1378896/76‑1 (MQ=255)                            aaTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGatat                                            <  2:2943588/76‑1 (MQ=255)                              tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTc                                          <  1:588064/76‑1 (MQ=255)                              tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTc                                          >  2:1144104/1‑76 (MQ=255)                              tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTc                                          <  2:2102833/76‑1 (MQ=255)                                taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTa                                        >  2:4233481/1‑76 (MQ=255)                                 aaCTCATTGGTTGATCTCGAGAACCGGGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTAc                                       <  2:2301741/76‑1 (MQ=255)                                  tcTCATTCATGGATTTCAAAAACCGTGAATTTCCATTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACa                                      >  2:1130753/2‑76 (MQ=255)                                  aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGCTATTCtccc                                      >  2:4872049/1‑73 (MQ=255)                                  aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACa                                      >  2:697389/1‑76 (MQ=255)                                  aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACa                                      >  2:4835745/1‑76 (MQ=255)                                    tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGa                                    <  2:1148296/76‑1 (MQ=255)                                    tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGa                                    >  1:2109867/1‑76 (MQ=255)                                    tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGa                                    >  1:4741905/1‑76 (MQ=255)                                    tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGa                                    <  1:3665420/76‑1 (MQ=255)                                    tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGa                                    >  1:3524018/1‑76 (MQ=255)                                    tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGa                                    >  1:2114726/1‑76 (MQ=255)                                     cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGaa                                   <  2:3910531/76‑1 (MQ=255)                                      attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGaaa                                  <  2:1400692/76‑1 (MQ=255)                                        tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGaaaaa                                >  2:4703546/1‑76 (MQ=255)                                           tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACa                             >  2:3904033/1‑76 (MQ=255)                                            ggATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAg                            <  1:1779844/76‑1 (MQ=255)                                              aTTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                          >  2:1305991/1‑76 (MQ=255)                                               tttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                         <  2:4853469/76‑1 (MQ=255)                                                ttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTg                        >  2:2824139/1‑76 (MQ=255)                                                 tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                       >  1:2345962/1‑76 (MQ=255)                                                   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                     >  1:3307008/1‑76 (MQ=255)                                                    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                    >  2:3094150/1‑76 (MQ=255)                                                     aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATa                   <  1:1770555/76‑1 (MQ=255)                                                      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                  <  2:1676496/76‑1 (MQ=255)                                                       aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                 >  1:3076218/1‑76 (MQ=255)                                                          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa              <  2:817766/76‑1 (MQ=255)                                                                    gTTAGAACAAAACAGCCGATTGTATGCCTGGGATATTCTACAGAAAACACACTCTGACATCCCGAAAATGATTaca    >  1:1874072/1‑74 (MQ=255)                                                                    gTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa    >  1:1868327/1‑76 (MQ=255)                                                                    gTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa    <  2:3344090/76‑1 (MQ=255)                                                                    gTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa    >  1:2368304/1‑76 (MQ=255)                                                                    gTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa    <  2:3309061/76‑1 (MQ=255)                                                                      tAGAACAAAACAGACGATTGTATGCCTGGGATGTTCTACAGAAAAAACAGTCTGTCATACCGAAAAGGAGTAGAgg  >  1:1145585/1‑76 (MQ=255)                                                                      tAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg  <  2:1210295/76‑1 (MQ=255)                                                                      tAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg  >  2:1208660/1‑76 (MQ=255)                                                                            |                                                                      ACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGG  >  NC_000913/4435607‑4435750

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