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, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	932,333	0	A	.	100.0%	234.5 / NA	57	intergenic (‑283/‑262)	trxB/lrp	thioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (30/27);  total (30/27)

CTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932196‑932465                                                                                                                                          |                                                                                                                                     cTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaagaa                                                                                                                                     >  2:1297528/1‑139 (MQ=255)        aaTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:634507/139‑3 (MQ=255)                            gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtata                                                                                                                  <  2:553449/130‑1 (MQ=255)                            gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtata                                                                                                                  >  1:553449/1‑130 (MQ=255)                            gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  2:1260108/1‑139 (MQ=255)                            gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  2:730119/1‑139 (MQ=255)                            gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATaa                                                                                                         >  1:631491/1‑139 (MQ=255)                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtata                                                                                                                  <  2:1000756/128‑1 (MQ=255)                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtata                                                                                                                  >  1:1000756/1‑128 (MQ=255)                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       >  2:1137666/1‑139 (MQ=255)                              ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAACCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       >  2:388385/1‑139 (MQ=255)                                        aaTAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                             >  1:459127/1‑139 (MQ=255)                                               gCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   >  2:80199/1‑110 (MQ=255)                                               gCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   <  1:80199/110‑1 (MQ=255)                                                  ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCAt                                                                                                           >  1:789716/1‑115 (MQ=255)                                                  ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCAt                                                                                                           <  2:789716/115‑1 (MQ=255)                                                              tatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       <  1:730119/139‑1 (MQ=255)                                                              tatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       >  1:433765/1‑139 (MQ=255)                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       <  1:175973/137‑1 (MQ=255)                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                       >  2:175973/1‑137 (MQ=255)                                                                tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                     >  2:1336141/1‑139 (MQ=255)                                                                  tGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaaa                                                                   <  1:1297528/139‑1 (MQ=255)                                                                   ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAAt                                                                  <  1:1137666/139‑1 (MQ=255)                                                                      aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         >  1:1443584/1‑113 (MQ=255)                                                                      aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:1443584/113‑1 (MQ=255)                                                                        gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              >  2:1448759/1‑73 (MQ=255)                                                                        gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaaca                                                                                                                              <  1:1448759/75‑3 (MQ=255)                                                                         ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   >  1:545715/1‑84 (MQ=255)                                                                         ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   <  2:545715/84‑1 (MQ=255)                                                                         ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         <  2:1373594/110‑1 (MQ=255)                                                                         ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                         >  1:1373594/1‑110 (MQ=255)                                                                         ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                            >  1:67401/1‑139 (MQ=255)                                                                                   aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  2:489755/1‑139 (MQ=255)                                                                                   aaTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  <  1:1260108/139‑1 (MQ=255)                                                                                      cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                <  1:1058835/59‑1 (MQ=255)                                                                                      cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                >  2:1058835/1‑59 (MQ=255)                                                                                        aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            >  2:1379831/1‑57 (MQ=255)                                                                                        aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                            <  1:1379831/61‑5 (MQ=255)                                                                                                     aCAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATg                                >  1:505715/1‑139 (MQ=255)                                                                                                         tCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  <  1:1423675/117‑1 (MQ=255)                                                                                                         tCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGa                                                  >  2:1423675/1‑117 (MQ=255)                                                                                                          ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAAt                                                                                           <  2:25877/75‑1 (MQ=255)                                                                                                          ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAAt                                                                                           >  1:25877/1‑75 (MQ=255)                                                                                                           cccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATa                          <  2:631491/139‑1 (MQ=255)                                                                                                             cTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   >  2:95980/1‑48 (MQ=39)                                                                                                             cTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                   <  1:95980/48‑1 (MQ=39)                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       <  2:999252/58‑1 (MQ=255)                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAcc                                                                                                       >  1:999252/1‑58 (MQ=255)                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATa                                                                                          >  2:1144478/1‑71 (MQ=255)                                                                                                               ggTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATa                                                                                          <  1:1144478/71‑1 (MQ=255)                                                                                                                gtgtTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTg                                                    >  2:635544/1‑108 (MQ=255)                                                                                                                gtgtTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTg                                                    <  1:635544/108‑1 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    <  1:1075885/55‑1 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACcatg                                                                                                    >  2:1075885/1‑55 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 >  2:49524/1‑74 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                 <  1:49524/74‑1 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                                   <  1:373388/120‑1 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                                   >  2:373388/1‑120 (MQ=255)                                                                                                                     ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACaa                >  1:128658/1‑139 (MQ=255)                                                                                                                      tGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATgcg                                                                                                              <  1:742773/44‑1 (MQ=38)                                                                                                                      tGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATgcg                                                                                                              >  2:742773/1‑44 (MQ=38)                                                                                                                       gCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAgg              <  2:433765/139‑1 (MQ=255)                                                                                                                              cATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGt       <  1:388385/139‑1 (MQ=255)                                                                                                                                   gAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  2:67401/139‑1 (MQ=255)                                                                                                                                   gAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct  <  2:128658/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     CTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCT  >  NC_000913/932196‑932465

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