Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,144	(T)7→6	100%	intergenic (‑356/‑384)	dinQ ← / → arsR	UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,138	0	T	.	100.0%	293.0 / NA	64	intergenic (‑350/‑390)	dinQ/arsR	UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (39/25);  total (39/25)

CAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648039‑3648233                                                                                                    |                                                                                                cAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                 >  1:160642/1‑100 (MQ=255)  aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                >  1:2052066/1‑100 (MQ=255)  aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                <  1:1243424/100‑1 (MQ=255)    tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              >  1:1235095/1‑100 (MQ=255)    tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              <  1:1742595/100‑1 (MQ=255)    tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              >  2:437297/1‑100 (MQ=255)     attgTATTCGGATAGATACGTAAGAGGATGTCTAGCTCAAGGCATTCAAAGCCAACAACTCAGGGCTGAAAAGTAAAAAACCGACGCAAAGTCGGttctt                                                                                             <  2:1312750/97‑1 (MQ=2)     aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                             <  2:1971848/100‑1 (MQ=255)      tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:954700/100‑1 (MQ=255)      tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:2325919/100‑1 (MQ=255)      tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:630509/100‑1 (MQ=255)       ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  1:270622/100‑2 (MQ=255)        gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          <  2:1251665/100‑3 (MQ=255)        gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          <  1:1846063/100‑3 (MQ=255)        gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          <  2:1905312/100‑3 (MQ=255)         tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         <  1:367133/100‑4 (MQ=255)           ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATc                                                                                      <  2:170556/100‑1 (MQ=255)             aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    >  1:636680/1‑100 (MQ=255)             aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    <  1:1338683/100‑1 (MQ=255)              gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCgg                                                                                   >  1:827841/1‑100 (MQ=255)                   gATGCGTTAGATGGTGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCg                                                                              <  1:486285/100‑1 (MQ=255)                     tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGa                                                                            >  1:1598260/1‑100 (MQ=255)                      gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGAc                                                                           <  2:1833265/100‑1 (MQ=255)                        gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACaa                                                                         >  2:748483/1‑100 (MQ=255)                        gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACaa                                                                         >  2:237748/1‑100 (MQ=255)                           aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGc                                                                      >  1:2271134/1‑100 (MQ=255)                           aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGc                                                                      >  2:583839/1‑100 (MQ=255)                             atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCtt                                                                    <  1:250538/100‑1 (MQ=255)                             atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCtt                                                                    <  2:1269439/100‑1 (MQ=255)                              tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTa                                                                   <  2:1567418/100‑1 (MQ=255)                               gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                  <  1:453863/100‑1 (MQ=255)                                 tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                >  1:1678279/1‑100 (MQ=255)                                    cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                             >  2:297441/1‑100 (MQ=255)                                     tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                            <  1:1303552/100‑1 (MQ=255)                                      aTCTCAGTCTAGTCAGAGACACCGACTCAGGGCTGGAAAGTAAAAAACCGACGCAACGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                           <  1:931539/100‑1 (MQ=18)                                        ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         >  1:1425181/1‑100 (MQ=255)                                         tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                        >  1:897268/1‑100 (MQ=255)                                          cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑CTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       >  1:1856077/1‑100 (MQ=255)                                           aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                      <  1:583842/100‑1 (MQ=255)                                            cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                     >  1:1010141/1‑100 (MQ=255)                                             tCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                    >  1:1423673/1‑100 (MQ=255)                                              ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   >  2:1640478/1‑100 (MQ=255)                                                aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   >  1:177163/1‑98 (MQ=255)                                                    agagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                             >  1:97173/1‑100 (MQ=255)                                                       gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  1:297443/100‑1 (MQ=255)                                                        ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                         >  1:2211410/1‑100 (MQ=255)                                                        ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                         >  2:518542/1‑100 (MQ=255)                                                         caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        <  1:1808077/100‑1 (MQ=255)                                                         caccaACGCATTGCTGGAAAGTAAAAAACCGACGCAAATGCGG‑TTTTTTACATCCGGATTCGGACAAGTCTTAATATGACGATGACCCAGTGAAAGtata                                        <  2:224159/100‑1 (MQ=25)                                                          accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                       >  1:927927/1‑100 (MQ=255)                                                          accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                       <  2:1678267/100‑1 (MQ=255)                                                            caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATAGaa                                     >  2:1580450/1‑100 (MQ=255)                                                                 cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  2:1181272/1‑100 (MQ=255)                                                                   gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATACATCGTc                               >  1:1896407/1‑99 (MQ=255)                                                                   gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                              >  2:1080465/1‑100 (MQ=255)                                                                     gCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                            <  1:1449356/100‑1 (MQ=255)                                                                      cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCGGGATTAGGACAAGGCTTAATATGACGATGACGCAGGGAAAGGATATAAATAGTCACTg                           >  2:446213/1‑100 (MQ=25)                                                                      cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                           >  1:202651/1‑100 (MQ=255)                                                                       tGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                          <  1:437300/100‑1 (MQ=255)                                                                        ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGAAAAGGCTTAATATGACGATGACCAAGTGAAATTATATACATTGTAAATGc                          >  2:1805161/1‑99 (MQ=25)                                                                          aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       <  2:1023616/100‑1 (MQ=255)                                                                           aaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                      <  2:663109/100‑1 (MQ=255)                                                                             gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    >  2:1008332/1‑100 (MQ=255)                                                                             gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    >  1:1815125/1‑100 (MQ=255)                                                                              tAAAAAACCGACTTTGCGTCGG‑TTTTTTACATCCGGATTCGGACAAGGATTAATATGACGATGAGCCAGTGAAAGTATATAAATAGTAAGTGCGatatat                   >  2:2044431/1‑100 (MQ=21)                                                                              tAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                   >  1:459785/1‑100 (MQ=255)                                                                               aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  1:1248094/1‑100 (MQ=255)                                                                               aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  2:900633/1‑100 (MQ=255)                                                                               aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  1:1027403/1‑100 (MQ=255)                                                                               aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  2:1578646/1‑100 (MQ=255)                                                                                aaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                 >  1:2241996/1‑100 (MQ=255)                                                                                 aaaaCCGACGCAAAGCCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                <  2:1913885/100‑1 (MQ=255)                                                                                     ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            <  2:1855980/100‑1 (MQ=255)                                                                                     ccGACGCAAAGTCGG‑TTTTTTACATCCGGAATCGGAAAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTACGATATATACCGa              >  2:38825/1‑98 (MQ=255)                                                                                       gACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGTCTATGACCCAGTGAATGTATATAAATCGTCACTGCGATATATACCAAAGt           >  2:1026433/1‑99 (MQ=255)                                                                                       gACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt           <  2:1055923/99‑1 (MQ=255)                                                                                        aCGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc         <  1:1181277/100‑1 (MQ=255)                                                                                            aaaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  1:428208/100‑1 (MQ=255)                                                                                              aGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctc   >  1:5589/1‑100 (MQ=255)                                                                                               gTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc  <  1:1277970/100‑1 (MQ=255)                                                                                                    |                                                                                                CAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648039‑3648233

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