Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	932,340	(A)8→7	intergenic (‑290/‑255)	trxB ← / → lrp	thioredoxin reductase/DNA‑binding transcriptional dual regulator Lrp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	932,333	0	A	.	100.0%	178.7 / NA	42	intergenic (‑283/‑262)	trxB/lrp	thioredoxin reductase/DNA‑binding transcriptional dual regulator Lrp
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (18/24);  total (18/24)

TGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932194‑932465                                                                                                                                            |                                                                                                                                     tGCTACTTAAATTTCCCGCTTTATAAGCCTATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaaga                                                                                                                                      >  1:865564/1‑140 (MQ=255)  gCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                                     >  1:118183/1‑140 (MQ=255)  gCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                                     >  1:122123/1‑140 (MQ=255)  gCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                                     <  2:707927/140‑1 (MQ=255)    tACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaa                                                                                                                                   >  2:712454/1‑140 (MQ=255)    tACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaa                                                                                                                                   <  2:617451/140‑1 (MQ=255)       ttAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  2:576921/140‑1 (MQ=255)         aaaTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:308906/140‑3 (MQ=255)           aTTTCCCGCTTTATAAGCCTATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            <  2:865564/140‑5 (MQ=255)           aTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            >  1:443437/1‑136 (MQ=255)                   cTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   <  2:118183/140‑1 (MQ=255)                     ttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                 <  2:122123/140‑1 (MQ=255)                         aaGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATgcgc                                                                                                             <  2:338367/140‑1 (MQ=255)                          aGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCa                                                                                                            <  1:712454/140‑1 (MQ=255)                           gCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCAt                                                                                                           <  1:750281/140‑1 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                        <  1:150466/140‑1 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                        >  2:484009/1‑140 (MQ=255)                                ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      >  2:77786/1‑140 (MQ=255)                                ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                      >  1:186648/1‑140 (MQ=255)                                  aaaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  1:112982/140‑1 (MQ=255)                                         gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                             >  1:787635/1‑140 (MQ=255)                                         gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                             >  2:1260/1‑140 (MQ=255)                                         gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                             <  2:559796/140‑1 (MQ=255)                                            tAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATa                                                                                          >  2:141942/1‑140 (MQ=255)                                             aaaCGCCCCTGTTAATGAATATCTGGCATGTCGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  1:704556/140‑1 (MQ=255)                                               aCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCa                                                                                       >  1:362463/1‑140 (MQ=255)                                               aCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCa                                                                                       >  2:469048/1‑140 (MQ=255)                                                    ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGttt                                                                                  >  1:457546/1‑140 (MQ=255)                                                               atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       <  1:141942/140‑1 (MQ=255)                                                                  tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaa                                                                    >  2:8628/1‑140 (MQ=255)                                                                    tGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       <  1:769430/135‑1 (MQ=255)                                                                    tGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       >  2:769430/1‑135 (MQ=255)                                                                        aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTAc                                                              >  2:669088/1‑140 (MQ=255)                                                                        aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTAc                                                              >  2:663284/1‑140 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:418471/140‑1 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  2:457546/140‑1 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            <  1:768142/140‑1 (MQ=255)                                                                           ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCACGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTa                                                           <  2:335594/140‑1 (MQ=255)                                                                               aCTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGc                                                       <  2:362463/140‑1 (MQ=255)                                                                                    aaaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  2:384891/1‑140 (MQ=255)                                                                                        cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGa                                              <  1:484009/140‑1 (MQ=255)                                                                                        cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGa                                              >  1:140812/1‑140 (MQ=255)                                                                                                          aTCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGa                            <  2:787635/140‑1 (MQ=255)                                                                                                            ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGcc                                          <  2:677237/124‑1 (MQ=255)                                                                                                            ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGcc                                          >  1:677237/1‑124 (MQ=255)                                                                                                            ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGAt                           <  1:746446/139‑1 (MQ=255)                                                                                                            ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGAt                           >  2:746446/1‑139 (MQ=255)                                                                                                             cccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAg                         <  1:444799/140‑1 (MQ=255)                                                                                                                 ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGc                     >  1:460137/1‑140 (MQ=255)                                                                                                                            aaaaCATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAc          <  1:8628/140‑1 (MQ=255)                                                                                                                               aCATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGt       <  1:384973/140‑1 (MQ=255)                                                                                                                                    cGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  2:186648/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                     TGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932194‑932465

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