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%	107.1 / NA	26	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 . (7/19);  total (7/19)

TAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAC  >  NC_000913/932201‑932457                                                                                                                                     |                                                                                                                             tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAACCATTCGAGgaaga                                                                                                                              >  2:97930/1‑133 (MQ=255)       tCCCTCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                        <  2:85605/133‑1 (MQ=255)       tCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                        <  2:492636/133‑1 (MQ=255)         ccGCTTTATAAGCCGATTAAATGATGAATAAACGACCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                      >  2:428905/1‑131 (MQ=255)          cGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacag                                                                                                                     >  1:105322/1‑130 (MQ=255)                   aGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                           <  2:454335/133‑1 (MQ=255)                   aGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                           <  2:105322/133‑1 (MQ=255)                   aGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                           <  1:361933/133‑1 (MQ=255)                        aTTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTCTTGCTTTGACAATCCCCTGGTCTTTTGCAAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATgcg                                                                                                      <  2:41180/133‑1 (MQ=255)                                  gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            <  1:368217/133‑1 (MQ=255)                                  gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            >  2:368217/1‑133 (MQ=255)                                  gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            <  2:342288/133‑1 (MQ=255)                                  gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            <  2:304161/133‑1 (MQ=255)                                  gaATAAACGCCCCTGGTAGTGAGTATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                            <  1:190298/133‑1 (MQ=255)                                   aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatgc                                                                                           >  1:390226/1‑133 (MQ=255)                                         cGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                     <  2:271756/133‑1 (MQ=255)                                           ccccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAAt                                                                                   <  1:428905/133‑1 (MQ=255)                                             ccTGTTAAGGAATATCTGGCATGTTGTACTTAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                 <  1:163737/133‑1 (MQ=255)                                                        atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCg                                                                      <  1:40396/133‑1 (MQ=255)                                                         tatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGt                                                                     >  1:37772/1‑133 (MQ=255)                                                          atCTGGCATGTTGTACTAAAAATCGATGCTTTACTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTg                                                                    <  2:157202/133‑1 (MQ=255)                                                           tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATACCCATGCATGTAGATCCCATGTATACCGTGc                                                                   >  1:180667/1‑133 (MQ=255)                                                                         cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACg                                                     >  1:449061/1‑133 (MQ=255)                                                                                 cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGt                                             >  1:469807/1‑133 (MQ=255)                                                                                 cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGt                                             <  2:491543/133‑1 (MQ=255)                                                                                 cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATCTTTACCGGGCTAGTGAAATCTTCGTATGGCGt                                             >  2:177827/1‑133 (MQ=255)                                                                                    tGTTTTGCTTAGACATTCCCCTGGGGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                          <  2:180667/133‑1 (MQ=255)                                                                                                   aTCCCCTGTTGTTTTTCGTAAAAATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                           <  1:177827/133‑1 (MQ=255)                                                                                                                   cGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATGCCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCca           <  2:390226/133‑1 (MQ=255)                                                                                                                   cGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCca           <  2:449061/133‑1 (MQ=255)                                                                                                                            tCGAGGAAG‑AAAAAACCAGTATTCTTATATGCCCATAACCAAGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAc  <  2:190238/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                             TAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAAC  >  NC_000913/932201‑932457

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